Aly 
, Pw op 
eats Ye 


ADVANCES... 
IN PHYSICS 


A QUARTERLY SUPPLEMENT 
of the 
PHILOSOPHICAL MAGAZINE 


EDITOR 
PROFESSOR B. H. FLOWERS, M.A., D.Sc., F.R.S. 


CONSULTANT EDITOR 
PROFESSOR N. F. MOTT, M.A., D.Sc., F.R.S. 


EDITORIAL BOARD 
SIR LAWRENCE BRAGG, 0.B.E., M.C., M.A., D.Sc., F.R.S. 
SIR GEORGE THOMSON, M.A., D.Sc., F.R.S. 


VOLUME 10 OCTOBER 1961 NUMBER 40 


PRICE per part 25s. plus postage 
PRICE per annum £4 15s. Od. post free 


PRINTED AND PUBLISHED BY TAYLOR & FRANCIS LID 
RED LION COURT, FLEET ST., LONDON, E.C.4 


REVISED CHEAPER EDITION FOR LIBRARIES AND SCHOOLS 


A History of Mathematics 


From antiquity to the early nineteenth century 


By J. F. SCOTT, B.A., D.Sc., Ph.D. 
Vice-Principal of St. Mary’s College, Strawberry Hill, Twickenham, Middlesex 


Author of The Scientific Work of René Descartes (1596-1650), 
Mathematical Work of John Wallis, D.D., F.R.S. (1616-1703), and other works 


Contents: Mathematics in Antiquity—Greek Mathematics—The Invention of 
Trigonometry—Decline of Alexandrian Science and the Revival in Europe—Mathematics 
in the Orient—Progress of Mathematics during the Renaissance—New Methods in 
Geometry—The Rise of Mechanics—The Invention of Decimal Fractions and of 
Logarithms—Newton and the Calculus—Taylor and Maclaurin, the Bernoullis and 
Euler, Related Advances—The Calculus of Variations, Probability, Projective Geometry, 
Non-Euclidean Geometry—Theory of Numbers—Lagrange, Legendre, Laplace, Gauss. 
This volume is intended primarily to help students who desire to have a knowledge of 
the development of the subject but who have too little leisure to consult works and 
documents. ‘The author has availed himself of the facilities afforded by the Royal 
Society and other learned Societies to reproduce extracts from manuscripts and many 
scarce works. 


Size 92” x 62”. 266 pp. Price 27s. 6d. plus postage and packing 2s. Od. 


Some Reviews of the First Edition 


“The invention of trigonometry, decimal fractions, logarithms and the calculus are each discussed 
clearly and concisely. ‘The book is easy to read for anybody who knows the elements of mathematics 
and, although not free from minor errors, can be strongly recommended.’’—British Book News, 
April 1958. 


“* Physicists will find that the development in applied mathematics are clearly set out, from ancient 
times, through that of Archimedes, to the mechanics of the sixteenth century when interest was 
revived. Significant advances made by Stevin, Galileo, Descartes, Huygens and others are stressed, 
and help the reader to appreciate what Newton achieved. There are useful appendices giving 
brief biographical notes on mathematical topics and terminology, followed by a bibliography.’’— 
Proceedings of The Physical Society, September 1958. 


“ 


...it has been written with clarity and balance, and the excellent printing helps to make it a 
pleasure to read.’”’—The Times Educational Supplement, 21 March 1958. 


“...his (Dr. Scott’s) wide knowledge of the material, his careful description of methods combine 
to provide an account which at times gives a sense of the excitement of discovery.’’—Nature, 
26 July 1958. 


“The printers and publishers are to be congratulated upon having produced such an attractive 
volume, .... We feel sure the book will be received with delight by all those interested in mathe- 
matical histories.’—BEAMA Journal, August 1958. 


“The work comes to life mainly because of his admirable use of the writings of mathematicians 
themselves, which vividly illustrates the great difficulties under which many of them laboured. 
This is not a book for the layman but both the student and anyone to whom figures are a fascination 
will find the subject clearly and pleasantly presented.’’—Technical Bookguide, March 1958. : 


Printed and Published by 


TAYLOR & FRANCIS LTD 
RED LION COURT, FLEET STREET, LONDON, E.C4 


CONTENTS 


Methods of Measuring Debye Temperatures and Comparison of Results 
for some Cubic Crystals. By F. H. Herssrer, National 
Physical Research Laboratory, Council for Scientific and Industrial 
Research, Pretoria, South Africa : j ‘ ; 


The Theory of the Interaction of Electrons with Lattice Vibrations in 
Metals. By G. V. Cuuster, Department of Mathematical ee 
The University, Birmingham 


Crystal Stability and the Theory of eee Part II. Piezo- 
electric Crystals. By W. Cocuray, Cer eee sak 
Cavendish Laboratory, bepaice 


313 


357 


401 


morse 


Methods of Measuring Debye Temperatures and Comparison of 
Results for some Cubic Crystals 


By F. H. Hersstein 


National Physical Research Laboratory, Council for Scientific and 
Industrial Research, Pretoria, South Africa 


CONTENTS PAGE 

§ 1. Iytropuctrion. 313 
§ 2. AssuMPTIONS OF THE DEBYE APPROXIMATION. 315 
§ 3. EXPERIMENTAL METHODS ror DETERMINING DEBYE TEMPERATURES. 316 
3.1. Specific Heat Measurements. 316 

3.2. Entropy Measurements. 317 

3.3. Measurements of Elastic Constants. 317 
3.3.1. Polycrystalline samples. 317 

3.3.2. Single-crystal samples. 319 

3.4. Diffraction Measurements. 319 


3.4.1. Measurement of Bragg intensities at a fixed temperature. 319 
3.4.2. Comparison of Bragg-intensity measurements at different 


temperatures. 321 

3.4.3. Measurement of temperature-diffuse scattering in powder 
patterns. 323 
3.4.4. Calculation of atomic vibration amplitudes. 324 
3.5. Measurements using the Méssbauer Effect. 324 
3.5.1. Introduction. 324 

3.5.2. Measurements of fraction of recoilless y-rays emitted or 
absorbed. 325 
3.5.3. Use of second-order Doppler shift of absorption maximum. 327 
3.5.4. Rayleigh scattering of Mossbauer radiation. 327 
3.6. Resistivity Measurements. 328 
3.7. Methods Based on Thermal Properties of Crystals. 328 
3.8. Other Methods. 329 
3.9. Applications to Alloys and Compounds. 329 

§ 4. THEORETICAL DIFFERENCES BETWEEN DeEByE TEMPERATURES 
DETERMINED BY DIFFERENT EXPERIMENTAL METHODS. 330 
§ 5. VARIATION oF DEBYE TEMPERATURE WITH TEMPERATURE. 335 
§ 6. Tae EXPERIMENTAL RESULTS. 336 
§ 7. Discussion oF RESULTS. 337 
7.1. Op Values at 0° and 300°x. 337 
7.2. Variation with Temperature. 342 
§ 8. ConcLUDING SUMMARY. 347 
ACKNOWLEDGMENTS. 348 
REFERENCES TO TABLES. 348 
REFERENCES. 351 


§ 1. INTRODUCTION 
Tur frequency spectrum of the coupled vibrations of the atoms in a crystal 
lattice has fascinated theoretical physicists since before the beginning of 
the century because of its intrinsic interest as a mathematical problem, and 
its importance in many branches of solid-state physics. Despite con- 
siderable effort, a mathematically exact calculation of the frequency 
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spectrum has not yet been achieved although a number of semi-empirical 
calculations, with varying degrees of approximation, have been made 
(Blackman 1955a for summary). The frequency spectra of a few elements 
have been determined experimentally by applying Born’s (1942-3) theory 
to the measured diffuse scattering of x-rays from Cu (Jacobsen 1955), _ 
a—-Fe (Curien 1952) and Al (Walker 1956 a). Similar determinations have 
been made for V by using the incoherent inelastic scattering of cold neutrons 
' (Stewart and Brockhouse 1958, Eisenhauer et al. 1958) and for Al by using 
the coherent inelastic scattering of cold neutrons (Carter et al. 1957, 
Brockhouse and Stewart 1958). Actual frequency spectra are so complex 
that approxiniate forms are still of considerable practical importance in 
calculating such quantities as specific heats, atomic vibration amplitudes 
in crystals, x-ray temperature factors and temperature diffuse scattering in 
elements and alloys and also for estimating thermal entropies. An 
important new application is to the study of the recoilless emission and 
absorption of nuclear gamma radiation in solids (the Mossbauer effect). 

The most widely used approximate frequency spectrum is the parabolic 
distribution first introduced by Debye (1912) and since known by his 
name. This frequency spectrum is usually described in terms of the 
characteristic or Debye temperaturet which is related to the maximum 
vibration frequency of the lattice v,, by Q0=hy,,/k (hk Planck’s constant, 
k Boltzmann’s constant). 

Current interest in Debye temperatures is in two directions. Firstly 
there is the desire to test the theoretical prediction that different methods 
of measuring Debye temperatures will give the same numerical values at 
temperatures approaching absolute zero. Experimental results obtained 
in the last few years confirm the theory (Alers and Neighbours 1959, 
§7.1). Secondly there is the need to use Debye temperatures (now 
referring to arbitrary temperatures of interest) to calculate many physical 
parameters of solids. There are a number of different methods of deter- 
mining the Debye temperature of a solid; these depend differently on the 
deviation of the true frequency spectrum from its assumed parabolic 
form and also differ greatly in convenience and ease of application. It is 
therefore important to determine whether the values of © obtained from 
various experiments can be used interchangeably and whether the differ- 
ences among the experimental values agree with theoretical predictions. In 
this paper the assumptions on which the Debye theory are based are first 
reviewed and the different experimental methods are then discussed. 
Values of Debye temperatures, determined by different techniques at 0°K, at 
room temperature and over a range of temperatures, have been taken from 
the literature for nineteen different cubic elements and these are critically 


/ 


+ In this paper the Debye temperature is denoted by ® when reference to a 
particular experimental method is not intended. Debye temperatures deter- 
mined by the various experimental methods are differentiated by subscripts: 
thus ©, refers to specific heat measurements, ©, to diffraction methods, and 
©, to resistivity measurements. 
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compared. A brief account of the principal results has been presented 
(Herbstein 1960). 

There have been a number of excellent reviews published recently 
dealing with various aspects of the vibrational spectra of crystals. One of 
the most comprehensive of these is the article by Blackman (1955a) in 
Handbuch der Physik, which covers both the Debye and Born—von Karman 
theories of specific heats and includes a detailed bibliography which will not 
be duplicated here. De Launay’s (1956) review covers similar ground, 
while Cole (1959) has presented a short and not too technical introduction 
to the subject, with emphasis on the use of diffraction methods for deriving 
the frequency spectrum. An excellent account of the effects of lattice 
vibrations on x-ray diffraction patterns of crystals has been given by James 
(1950), while more detailed treatments are due to Laval (1958) and to 
Slater (1958) who has reviewed the general theoretical background. 


§ 2. ASSUMPTIONS OF THE DEBYE APPROXIMATION 


In the Debye approximation it is assumed that the crystal lattice can be 
treated as an isotropic elastic continuum. It follows that the frequency 
distribution function is 

N(v) dv =41 G + a) Viidv vv, 

=( V>Vy 

where N(v) is the number of frequencies between v and v+ dy, C, and C;, are 
the velocities of the longitudinal and transverse elastic waves respectively 
and V is the volume of the crystal. It was assumed that C, and C, are 
invariant with respect to direction in the crystal and wavelength of the 
elastic wave (i.e. there is no dispersion). The atomic structure of crystals 
was then taken into account in the theory of specific heats by setting the 
same lower limit to the wavelengths of the longitudinal and transverse 
waves. Consequently there will be different upper limits to the frequency 
distributions of the transverse and longitudinal waves and an avarage 
maximum frequency v,, must be chosen. This must be such that the 
number of normal vibrations is 3, where JN is the number of atoms in the 
crystal. Thus 

le avidv=3N or v,=(9N/a)® Se at sures L) 

0 

where a=47V(O,-3+ 2C,*). 

As mentioned above the Debye temperature is then defined in terms of the 
cut-off frequency by 9 =hv,,/k. As the method of averaging the maximum 
frequencies for the different vibrational modes is not the same in specific- 

_heat theory as in the theory of the x-ray temperature factor, the respective 
Debye temperatures will differ slightly. This is discussed in more detail in 

4. 

The basic assumptions of the Debye approximation—simple harmonic 
vibrations, elastic isotropy, no dispersion—are by no means fully realized 
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in practice and thus lead to restrictions on the applicability of the theory. 
Anharmonicity of atomic vibrations and the existence of dispersion appear 
to affect only the details of the vibration spectrum, but the condition of 
elastic isotropy suggests that it is best fo restrict application of the Debye 
theory to cubic crystals (even among cubic crystals only W, and to a 
lesser extent Al, are approximately elastically isotropic). Elastic aniso- 
tropy is compensated for semi-empirically to a considerable extent as the 
frequency spectrum is in fact an average taken all over all directions in the 
crystal lattice. A further restriction is to crystals with only one atom in the 
primitive unit cell; this is because only the acoustic branch of the spectrum 
is included. The setting of a lower bound to the elastic wavelength in the 
lattice, independent of direction, involves the approximation of replacing 
the first Brillouin zone of the crystal by a sphere. This holds fairly 
accurately for body-centred and face-centred cubic crystals but less well 
for simple cubic crystals (Slater 1939). 

Reliable vibrational spectra have been determined for a number of 
elements (e.g. Al, Cu, V) and these are very different from the corresponding 
parabolic forms and, indeed, appreciably different from the spectra calcu- 
lated by lattice dynamics on the basis of various simplifying approximations. 
Despite the lack of realism in the basic assumptions of the Debye approxi- 
mation, and the resultant differences between the true vibrational spectra 
and the parabolic approximation, the effects are not very severe for such 
properties as specific heat, vibrational entropy and x-ray temperature 
factor where integration over the whole frequency spectrum is required. 
This is a natural consequence of the averaging process implicit in the 
integration. A further compensating factor is provided by the semi- 
empirical nature of the Debye approximation because the parabolic 
frequency distribution, and hence ©, is chosen to give the best fit between 
theory and the particular experimental result obtained. On the basis of 
these arguments, the Debye approximation would be expected to be much 
more successful in practice than the inappropriateness of its initial assump- 
tions might lead one to suppose. This is borne out both by past experience 
and by the results summarized here. 


§ 3. EXPERIMENTAL METHODS FOR DETERMINING DEBYE TEMPERATURES 
3.1. Specific Heat Measurements 
The classical method of determining the Debye temperature is from 
measurement of the specific heat. It is necessary to subtract the electronic 
specific heat from the measured values and to convert the measured C p to 
C, using the thermodynamic equation C,,—C,,=6?V7/y where B is the 
coefficient of volume expansion and y the isothermal compressibility. 
Debye’s equation for the heat capacity of 1 g-atom of solid is 
73 x preP 
aes Ep ee eT 
Op lg (1)? 2 
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where NV is Avogadro’s number, & Boltzman’s constant, 


Op 


h 
Serna and @p= 220 


wae 

The usual procedure is to calculate ©, from eqn. (2) using measured 
values of C,. Plots of @p against 7’ show that ©, is not constant as it 
would be if Debye’s assumptions were valid. The variations of @p with T 
can be explained in terms of the Born-von Karman theory of lattice 
vibrations. For purposes of comparison with the results of other experi- 
mental methods, the variation of ©, with 7’ makes a choice of Op difficult. 
Blackman (1955a, table 9, p. 364) has used 7’ = 1@, as a convenient standard 
temperature as “This temperature is aninegnae high to give fair agreement 
over most of the specific heat curve, and sufficiently low for errors in the 
(Cp—C,) correction to be small”’. The correction for the electronic specific 
heat will also be negligible at this temperature. Blackman’s usage has 
been followed in selecting the @p values given later in the tables, where they 
are nominally referred to a temperature of 300°K. Recent measurements 
of @p for many elements over a range of temperatures have been summarized 
by De Sorbo (1954) while ©, values determined at very low temperatures 
have been given by Alers and Neighbours (1959). 


3.2. Entropy Measurements 


@©p values can be derived from calorimetric entropy measurements 
using the equation (Lumsden 1952): 


3x2 34 ae 8 
a — eS Sh ee ee eee Se 3 
nee (4 Sine + 70 — 2240 36288 1689600 * ) (3) 


where x= @,/7. Lumsden has also reported values for a number of metals 
which are included in the tables of experimental results. The measured 
entropy values used in (3) must, of course, be corrected for electronic and 
other non-vibrational contributions. One use of eqn. (3) is for the calcu- 
lation of vibrational contributions to the entropy of alloys, using Debye 
temperatures measured in other experiments. 


3.3. Measurements of Elastic Constants 


3.3.1. Polycrystalline samples 


Polycrystalline samples of small grain size are effectively isotropic if free 
from preferred orientation so that Debye’s original derivation (1912) of © 
from the elastic constants of an isotropic continuum applies. Different 
expressions are obtained for ©, and @p, the Debye temperatures appropri- 
ate to diffraction and specific-heat measurements respectively, because 
different methods of averaging over the longitudinal and transverse 
vibrations are used (see §4). It is convenient to express these values in 
terms of Young’s modulus (Z) rather than the bulk modulus (K) originally 
used by Debye because the equations then depend on functions of Poisson’s 
ratio (c) which vary very slowly in the region of interest. The equations 
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given below were derived from those given by Debye (cf. James ( 1950)) by 
making the substitution K=H/3(1—2c) before averaging over the 
longitudinal and transverse values of @. We obtain 


ha/3 (3N\'8 EN 
Oe (ey pei) es 


k 4a 
with 3—50 ua 
ii(a)= Feel ’ . . . ° . (4a) 
h ON 1/3 Fil 
®>- F(Z) mils pire S27): eet) 
ae f,(o)= (92) (282 (8— 5a) +9°*(1 —20)} ee 


In these equations m is the atomic weight (for an alloy the mean atomic 
weight should be taken), p is the density and NV is Avogadro’s number. 
f,(c) and f,(c) are insensitive functions of o in the range 0-2<oa<0-5 (see 


Fig. 1 
‘420 ‘480 
f,(¢) 
1) -410 -470 
*400 ‘460 
0-2 O3 O-4 O'5 


Poissons ratio @ 


Variation of f\(o) (eqn. (4a)) and f,(o) (eqn. (5 a)) with o (Poisson’s ratio). 


fig. 1) in which all the values for metals cited by Koster (1954) lie. Conse- 
quently simpler approximate forms of (4) and (5) can be given; for 
example @,* 12 (Ha/A)"? for face-centred cubic materials, where-E is in 
kg/mm?, A is in atomic units (i.e. 20 = 12) and a is the lattice parameter in 
angstroms. 
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Values of the elastic constants can be measured either by static or dynamic 
methods. Useful summaries of measured values of elastic constants for 


polycrystalline metals and alloys have been given by Koster (1948 a, b, 
1954), 


3.3.2. Single-crystal samples 
The Debye temperature can be expressed as 


where V, is the volume per atom and (C) is the mean velocity of the elastic 
waves in the lattice and is defined by 


BYO) ey Oi 1928 


in which the C; are the three velocities of propagation of plane elastic waves 
in the crystal, averaged over all directions in the lattice. A fairly large 
number of methods of calculating (C) from the measured elastic constants 
have been devised and are summarized by Blackman (1955a). The 
present calculations of ©, from measured ¢, values were made using the 
methody of Quimby and Sutton (1953) and Sutton (1955). 

Elastic constants are usually measured by ultrasonic pulse-velocity 
methods and comprehensive reviews of theory and results have been given 
by Hearmon (1946) and Huntingdon (1958). 


3.4. Diffraction Measurements 


Much of the subject matter of this section is discussed in detail by James 


(1950, chapter V) where derivations of many of the equations used here will 
be found. 


3.4.1. Measurement of Bragg intensities at a given temperature 


For a monatomic cubic lattice the structure factor of the Bragg reflection 
hkl is given by 
K F? (hkl) =f? (hkl) exp (— 2M). heer (0) 


Equation (6) is applicable to diffraction measurements in general if f (Ak) is 
understood to denote the atomic scattering factor f,(hkl) for x-ray diffrac- 
tion, the atomic potential factor : 


e (Z—f,(hkl) 
MN ral (sin 8/A)? ) 
for electron diffraction or the nuclear scattering length 6b for neutron 
diffraction. F (hkl) is the structure factor derived from the corresponding 
measured intensity by standard methods in each case. 2 = 2B (sin?6/A”) 
is the Debye-Waller temperature factor, 0 is the Bragg angle, A is the wave- 
length of the radiation used and K is a normalizing factor depending both 


i ee a 
+ The calculations were done on the Stantec Zebra electronic computer in 
our laboratory, using a programme prepared by Miss 8. A. Gous. 
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on the experimental arrangement and the symmetry of the lattice. An 
accurate experimental value of K can be obtained relatively easily for 
neutron diffraction (Bacon 1955, p. 93) but this is much more difficult 
for x-ray or electron diffraction. 

From (6) it follows that 
sin? 0 


- (7) 


Thus a plot of In (f2/F2) against sin? 0/A2 will give both B, from the slope of 
the plot and In K, from the intercept at sin?@/A?=0. The value of ©, is 
obtained from the measured value of B by using the Debye-Waller formula 
(Debye 1914, Waller 1923) 


is Rae 
In (3) =Ink+2B 


12h? T (x) 
2B = —_——— 8 
mk On? ( 
where m is the mass of the atom, 7' is the absolute temperature of the 
sample, x= ©),/T7' and 
l(*udu | &@ 
¢) = ae ee a ee 
Way= 2 SE +5 (9) 


u(x) varies from 1-000 to 1-164 for 0<x< 2-5 and thus (8) can be solved for 
©, by a series of successive approximations. 

The great majority of values of © ,,derived from diffraction measurements 
have been obtained by x-ray diffraction. Intensity measurements in 
electron diffraction are complicated by theoretical and experimental 
difficulties so that this method has been little used. Only one independent 
value of @,, has so far been derived from neutron-diffraction intensity 
measurements (Henshaw 1958) but this method does have one definite 
advantage over others which is discussed below. 

There are a number of problems involved in the use of eqn. (6) which must 
be mentioned. The first of these is the necessity of correcting Bragg 
intensity measurements for a contribution from the temperature-diffuse 
scattering (TDS). The form of this correction, which must be applied to 
intensity measurements in all types of diffraction experiments, has only 
been published recently (Nilsson (1957), Chipman and Paskin (1959b)). 
In x-ray diffraction measurements, for example, the Bragg peaks are 
superimposed on a background of Compton scattering and temperature- 
diffuse scattering. The former is a smooth function in diffraction space 
but the latter has maxima at the same positions as the Bragg reflections. 
Thus the true intensity of the Bragg reflection must be obtained by sub- 
tracting the calculated TDS peak area from the measured total peak area, 
taking both over the same angular range. The importance of this correc- 
tion varies with the material investigated ; for single-crystal measurements 
on NaCl, Nilsson estimates that application of the correction reduces 
Renninger’s (1952) measured ©, from 319°K to 302°K. 

There is a further difficulty which applies to x-ray and electron diffraction 
but not to neutron diffraction. This is that the atomic scattering factors 
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are themselves theoretical quantities calculated from necessarily imper- 
fectly known electron distributions in atoms; in x-ray diffraction it is 
further necessary to correct the tabulated values for anomalous dispersion 
(James 1950, chapter IV, Parratt and Hempstead 1954, Dauben and 
Templeton 1955) and this correction is sometimes difficult to apply accu- 
rately. Neutron diffraction has the great advantage that the nuclear 
scattering length does not vary with scattering angle. Although the 
nuclear scattering lengths are empirical quantities, accurate values are 
available for most elements (Bacon 1955). 

Borie (1956) has proposed a method of determining both the atomic 
scattering factor and the Debye temperature from the x-ray powder 
pattern of an element. This method requires measurements of the 
integrated intensities of the Bragg peaks and the intensity of the diffuse 
background, using monochromatic radiation. The temperature diffuse 
contribution to the background is obtained by subtracting the calculated 
Compton scattering from the measured values. The TDS (for a powder 
pattern) can be written as 


Imps ("5") = f2(1 exp (-2)) 4 (=") ee 10) 


where G (sin 0/A) is a radial function in reciprocal space whose form has been 
given by Warren (1953), Herbstein and Averbach (1955) and Paskin (1959)T 
for a variety of situations. Thus f?[1—exp(—2/)] can be obtained as a 
function of (sin@/A) from the TDS while the Bragg intensities give 
f?exp (— 2M) as a function of (sin 6/A). The sum of the two curves is f? and 
their ratio is [exp(2M)—1]. This method has been applied with reason- 
able success to copper. One difficulty is the calculation of the Compton 
scattering, which amounts to about 30% of the background for copper, and 
which is known with an accuracy of about 10%. Chipman (1960) has esti- 
mated that the resultant uncertainty in ©), of copper is + 10°K. 


3.4.2. Comparison of Bragg-intensity measurements at different temperatures 
Measurements at different temperatures can be used as follows: 


; : se 
hin. Ge tai (itl) = Kj (=") exp| -B, (Fi ve 
. 1 


, os 
ab Tack, (ED) KG (= *) sae be B (Fe) | 
A /2 eA 
(It is assumed that any changes in Bragg angle due to thermal expansion 
have been allowed for in the derivation of structure factors from measured 
intensities.) If (7) and (8) are now used, ©), being treated as varying with 
temperature, we obtain 


gf, | _ eee ae W+ke+e 
2 | | 7 tae i ee 


x 3h? ato. xe tf eee Saeed) 


2mk |ag? Oy” 44” Oyn” 


+ See note added in proof. 
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where (a) is defined in eqn. (9). g is a correction factor, differing only 
slightly from unity, which allows for the change of atomic scattering factor 
with angle in x-ray and electron diffraction. For neutron diffraction g =1.. 
Thus the slope of the plot of In(gF',/F2) against (A? + k° +?) will give an 
equation relating @,, and © ,,; this can be solved directly if © yy 1s known 
from a separate experiment. Alternatively it has often been assumed that 
it is permissible to take ©, as constant over small temperature intervals. 
The smaller the temperature interval the better this assumption, but also 
the smaller the difference in F, and Ff: clearly very accurate intensity 
measurements are necessary for this method to work well. 


Fig. 2 


600 


400 


300 


100 200 300 400 500 600 700 800 900 1000 
Tack 


®, of aluminium versus temperature obtained from experimental intensity 
measurements. The various curves were obtained by choosing ©,, at 
300°K equal to 380°, 390°, 400°, 410° and 420°K respectively. This figure 
taken from Chipman (1960). 


When measurements are made at a number of temperatures it is usual 
to compute all the results with reference to a standard temperature 7, 
usually room temperature, and then plot (h?+k?+/?)In (gF,/F.) against 
T',. Some workers have then proceeded by assuming a particular form for 
the temperature-dependence of ©, such as that given by assuming that | 
the temperature-dependence of @,, is due to thermal expansion only (see 
§5). This procedure fails when its basic assumption is violated, as has been 
demonstrated for silver (Herbstein 1961). A more objective method of 
analysing the results has been proposed by Chipman (1960) who pointed out 
that the dependence of ©,, on 7’ is itself dependent on the value chosen for 
©. This is illustrated in fig. 2. The latitude introduced in this Way is 
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limited by the fact that only certain forms of the curve of @ m against 7’ are 
physically reasonable. Choice of the appropriate value of @y,, is compli- 
cated by the fact that these curves differ significantly in shape only at 
fairly low temperatures, and thus experimental data below 7’~ @,, are 
essential for this method to work satisfactorily. The form of the.curve of 
©, against 7’ has not been established theoretically but qualitative con- 
siderations suggest that it will have a small negative slope at fairly low 
temperatures and will rise at very low temperature in a manner somewhat 
similar to that found for plots of @, (specific heat) against 7. A careful 
analysis would be necessary if this method were to be applied to materials 
where an appreciable difference might be anticipated between the ©, 
values at two temperatures of measurement (e.g. measurements at 300°K 
and 80°K on Si or diamond). 

An interesting variation on the above method has been used by 
Brockhouse (1955) who measured the change with temperature of the 
intensity of elastically scattered incoherent neutrons from vanadium and 
determined the Debye temperature by applying the Debye-Waller 
expression in its appropriate form. (The cross section for incoherent 
scattering from vanadium is very much larger than that for coherent 
scattering.) The.variation of intensity with temperature in reflection and 
transmission electron-diffraction patterns has also been used to estimate 
vibration amplitudes, and hence Debye temperatures, in a number of 
metals (Coster and van Zouten 1939, Yamaguchi and Katsurai 1948, 
Menzel-Kopp and Menzel 1955, 1956). The method does not appear to 
give very accurate quantitative results but interesting qualitative results 
have been obtained (see § 7.2). 


3.4.3. Measurement of temperature-diffuse scattering (TDS) in powder - 
patterns t 


Measurement of the temperature-diffuse scattering from single crystals 
can yield sufficient information to permit derivation of the frequency 
spectrum of the crystal (Born 1942-3, Jacobsen 1955); if the sample is in 
the form of a powder, experimental measurements are much easier but the 
amount of information that can be obtained is correspondingly less. 
Chipman and Paskin (1959a) have shown that the temperature-diffuse 
scattering from a cubic powder can be written (for 7’ > ©): 


Ipps=N frexp (— 2M) { (2p) O.+ 2M) St eee (2p) Cyth 


(12) 
‘N is the number of atoms, f the atomic scattering factor, 2M the usual 
Debye-Waller factor, 2M a modified Debye-Waller factor, and the C;’s — 
functions that depend on the details of the calculation of the contribution of 
the Ith order phonon process to the TDS. The calculation of the C/s 
has been discussed in detail by Paskin (1958, 1959). The values of f? and 
Pai eee ne ener 6 eS ee 


+ See note added in proof. 
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2M can be taken from other results or derived from measurements on the 
Bragg reflections and then 2M;, obtained by comparison of eqn. (12) with 
experiment. For both copper and lead Chipman and Paskin (1959 a) 
find that 2M, <2M and suggest that this may be due to differences in 
methods of averaging over the velocities of the elastic waves (see § 4). 
At present the determination of Debye temperatures from TDS is 
associated with both experimental and theoretical difficulties and it 
seems that one of the other x-ray methods is to be preferred for this purpose. 

Debye temperatures have also been deduced from the change in energy of 
cold neutrons after being inelastically scattered by solids (Brockhouse 
(1955), Eglestaff (1959)). These values were actually obtained in pre- 
liminary experimental work on inelastic scattering of neutrons and it is 
clear that future efforts in this field will be devoted to determining the true 
form of the vibration spectrum without recourse to the Debye approxi- 
mation. 


3.4.4. Calculation of atomic vibration amplitudes 


The measured B factor (Debye-Waller temperature factor) is related to 

the mean square displacement of an atom from its average position by 
B= 8n'a?/3. 

This equation applies only to monatomic cubic crystals. Using eqn. 
(8) one can write (Lonsdale 1948) : 

4-364 x 10°7' 
= A2 

46,2”) 

where A is the atomic weight (17C = 12). Lonsdale (loc. cit.) has given values 
of (a)! for a number of elements, using tabulated values of @,,;. Lonsdale 
has particularly emphasized that even in cubic crystals the vibration 
amplitudes are not isotropic and has shown how the components in particu- 
lar directions due to particular elastic waves can be calculated from a 
knowledge of the elastic constants. Measurements of integrated diffracted 
intensity (or TDS from powder samples), to which attention has been 
restricted here, give only the average vibration amplitude without per- 
mitting distinctions among directions in the crystal lattice. 


a 


3.5. Measurements using the Mossbauer Effect 

3.5.1. Introduction 

There is often an appreciable probability that the y-ray emitted from a 
nucleus in a crystal lattice will have the exact energy of the nuclear transi- 
tion because the recoil energy of the emitting atom will be taken up by the 
lattice of the host crystallite. As a similar situation can occur on absorp- 
tion, resonance is obtained if emitting and absorbing nuclei are under 
identical conditions. Resonant absorption of recoilless y-rays was first 
demonstrated by Méssbauer (1958 a, b, 1959) using the 129Kev y-ray of 
Tr, while more recent work has shown that similar effects occur in a 
number of other nuclei, including *7Fe, "Sn and Au. The line width is 
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so narrow that resonance can be destroyed by the (first-order) Doppler 
shift resulting from a small uniform relative velocity of the order of milli- 
metres per second between source and absorber. Accurate measurements 
of line shapes and positions can be made in this way. Conditions are 
particularly favourable for *’Fe and an energy resolution of about 1 part in 
10” can be obtained for the 14-4Kev y-ray from this nuclide. This remark- 
_ able resolving power has led to the execution of a number of striking 
experiments, including measurements of the gravitational red shift 
(Pound and Rebka 1960 b) and other tests of the general theory of relativity 
(Sherwin 1960) and has also made possible measurement of the electrical 
and magnetic fields at nuclei in solids. A comprehensive review of theory 
and applications of the Mossbauer effect has been given by Cotton (1960) 
while less detailed accounts have been given by Burcham (1960) and Lustig 
(1961). For present purposes the important feature of the Méssbauer 
effect is that it provides a number of different ways of measuring the 
vibration amplitudes of the emitting (or absorbing) nuclei, which are 
discussed in more detail below. 


3.5.2. Measurements of fraction of recoilless y-rays emitted or absorbed 

For a nucleus bound in a crystal lattice, the probability of emission or 
absorption of a y-ray without energy transfer to or from the lattice is 
governed by the Debye-Waller factor familiar from diffraction experiments. 
Thus measurement of the fraction of y-rays emitted or absorbed without 
recoil will give the Debye temperature of the lattice. A simplified account 
of the theory has been given by Lipkin (1960) and more detailed derivations 
by Visscher (1960) and Singwi and Sjolander (1960). We shall give here a 
simplified summary of the results of these calculations, following the 
paper by Singwi and Sjélander. If the emitter and absorber are identical 
but at different temperatures 7’. and 7’, and moving with a uniform relative 
velocity v, then the self-absorption cross section for a crystal with a cubic 
Bravais lattice is given by 


da(S) = 


Ret © (2We+2Wa)™_. 
«agit 5 Cree aa) See 5 MUD 


m0," 


exp | —(2We+2Wa)] 


n=1 

In this expression s=(v/c) 2) where H, is energy of the y-ray emitted, a 
is the resonance absorption cross section, ’is the natural width of the 
excited state of the nucleus, ,,(s) are functions given in detail by Singwi and 
Sjdlander and which are not required for the present discussion. 

In the Debye approximation 

Ey 6T'e 

= a 
© 2 — hi Om" Pn) 
SS  — —— —— ee ne ce a a Leen 

+ It should be noted that the second-order Doppler shift of the lines due to 
temperature difference between source and absorber (see below) has been 
neglected in deriving this equation. 
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with a similar expression for 2W,. (x) has been defined previously in 
eqn. (9). The analogy between eqn. (14) and the Debye-Waller formula 
for 2M (see eqn. (8), remembering that 21/=2B (sin? @/A?)) is clear if one 
writes the single equation 
67 y(a) 

#032 ~ 
For x-ray diffraction H,, is the recoil energy given to the free atom by a 
photon of energy H(=/c/A) scattered through a total angle 26. In this 
case 


2M=E (15) 


E? 2h? sin? 6 
Lg ee Oe ed ez: 


For emission of y-rays the photon of energy H gives a recoil momentum H/c 
to a free atom of mass m and thus the kinetic energy of recoil is H?/2mc?. 
Equation (13) shows that the absorption cross section consists of a sharp 
resonance peak of width 21 centred abouts=0. For 7'’e=7' the first term 
is 
oa’ (0) =40,exp(—4W). 


The resonance self-absorption cross-section is thus diminished by a factor 
exp(—4W). Measurements of the self-absorption cross section will there- 
fore give the Debye temperature at the temperature of measurement. <A 
number of measurements of this type have been made on a—Fe using the 
14:4Kev y-ray from *’Fe. In the earlier experiments agreement between 
theory and experiment was rather poor (for example Pound and Rebka 
(1959) reported that only about 20% of the emitted y-rays were recoilless, 
whereas the calculated fraction is 72° for @ =355 °K). More recent work 
has increased the experimental value to 60°% (Cordey-Hayes et al. 1960, 
Clark and Robson 1960) although the fractions of recoilless emissions and 
absorptions are still less than would be predicted by the Debye model with 
@p=420°K. Rather similar discrepancies have been reported for mea- 
surements on Au, Pt, Ir and Os (Nagle et al. 1960). It has sometimes been 
suggested that these discrepancies are due to real differences between 
©, and ©, for the materials investigated. However, both theory and 
experiment show that the differences between ©, and @,, values are seldom 
greater than 5% and the discrepancies reported must be ascribed to other 
effects. These may be partly a result of experimental difficulties but it has 
also been suggested (Bussieré de Nerey, Langevin and Spighel 1960) that 
Lamb’s theory (1939) predicts too large values for the fraction of recoilless: 
emissions or absorptions when the lifetime of the excited state is larger than 
the time for de-excitation of the lattice. 

It is not essential that emitter and absorber should be at the same 
temperature. In one of his original experiments Méssbauer (1958 a) 
determined the variation of the absorption cross section of 1Ir when the 
absorber was held at 88°K and the source temperature varied between 
88°K and 380°K. Visscher (loc. cit.) has shown that these experimental 
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results can be explained if © =316°K, which is somewhat larger than the 
specific heat value of ©, = 285°K. 

The higher-order terms in eqn. (13) give the shape of the wings to the 
sharp resonance absorption line. If 2W<1 then all terms in the sum 
except the first are negligible. In this case, as Visscher has pointed out, 
the shape of the wings is related in a simple way to the energy spectrum of 
the crystal vibrations. The experimental situation is, however, rather 
difficult and it seems probable that slow-neutron scattering provides an 
easier path to the vibration spectrum. 


3.5.3. Use of second-order Doppler shift of absorption maximum 


An alternative method of deriving Debye temperatures depends on the 
temperature-dependence of the position of the absorption maximum. 
The resonance absorption peak will shift because of the change in the 
energy of the lattice due to the change in nuclear mass consequent on the 
emission or absorption of a y-ray photon. The magnitude of this shift has 
been given by Pound and Rebka (1960a) and by Josephson (1960) as 


av __ Cp 
or = 2mC? 


where C, is the specific heat of the lattice at constant pressure. For 
a-iron the fractional frequency shift will be 2-2 x 10-15/°K at 300 °K and 
0-73 x 10-45/°K at 80 °K. Pound and Rebka (1960a) measured Av/y as a 
function of the temperature difference between source and absorber and 
found that their results could be explained if ©)p=430°K. It is worth 
noting that it is @p which is measured in this experiment, and not ©, as in 
the other Mossbauer effect experiments. 


3.5.4. Rayleigh scattering of Méssbauer radiation 


The methods described above can only be applied to solids containing 
nuclei which emit suitable y-rays. However, the method can be extended 
to other solids by using the high energy discrimination of the Méssbauer 
effect to measure the fraction of recoilless y-rays which have been Rayleigh 
scattered (i.e. without phonon exchange) by a solid at a given temperature 
(Tzara and Barloutard 1960). This fraction is given by the Debye-Waller 
factor in the form of eqn. (14). Tzara and Barloutard, using the 23-8kev 
radiation from 14°Sn and a !°Sn absorber, have determined Debye tem- — 
peratures of Pt, Al, graphite and paraffin wax in this way and obtained fair 
agreement with the expected values. 

Determination of Debye temperatures by various methods employing 
the Mossbauer effect may in many instances be simpler than the older 
techniques. For example Singwi and Sjélander (loc. cit.) have suggested 
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that a change in the absorption of Méssbauer radiation at the super- 
conducting transition temperature could provide experimental evidence 
for a change in the phonon spectrum at this temperature. Specific heat 
measurements for Nb have suggested that ©, ~ 243°k in the supercon- 
ducting phase and ~231°K in the normal phase (Boorse et al. 1960). 
Earlier neutron diffraction intensity measurements (Wilkinson et al. 1954) 
were able to show that the change in @p could not exceed 20°K. The 
Mossbauer method would, of course, only be applicable to superconductors 
with suitable isotopes. In the case of disordered alloys, use of the 
Mossbauer effect should give information about the vibration amplitudes 
of the emitting (or absorbing) atoms which cannot be obtained by diffraction 
measurements (see §3.9). 

At the time of writing (February, 1961) most experiments using the 
Mossbauer effect for determining Debye temperatures have been essentially 
exploratory and no accurate values of ©; have yet been reported by this 
technique. But the subject is still in its infancy and one must expect many 
striking developments in the future. 


3.6. Resistivity Measurements 


Using the Bloch (1928, 1929) theory of the temperature dependence of 
the electrical resistivity of metals, Griineisen (1933) related the resistivity 
to a characteristic temperature Oy. In general no close agreement is to be 
expected between the values of ©) and @p as, according to Bloch’s theory, 
Op depends on the longitudinal lattice vibrations only. This was pointed 
out by Blackman (1951) who, however, noted that there is actually a 
remarkable correspondence between ©, and © p values for many metals. 
Blackman states: “It is true that the assumptions of the Bloch theory do 
not apply strictly to most of these cases; why the deviations from these 
_ assumptions should lead to Og ~ @p is a matter which is in great need of 
clarification.” Because of the lack of a sound theoretical basis for the Ox 
values, they will not be considered further here. 


3.7. Methods Based on Thermal Properties of Crystals 


Because of the dependence of the thermal properties of crystals on the 
vibration amplitudes of the atoms, it should be possible to derive © values 
from such measured quantities as self-diffusion coefficients, thermal- 
expansion and conductivity coefficients and melting points. This field 
has been extensively reviewed by Cartz (1958b). For example, if the 
melting point is used, then Lindeman’s theory of the melting of a solid can 


be used to give 
e (ertace, [i 
hy | att | 


Here 7’; is the melting point in °K, A is the atomic weight, r is half the 
distance between neighbouring aroms, and L, is a constant determined 
empirically from @p values for F.C.C., B.C.C. and NaCl-type crystals. 
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Cartz (loc. cit.) used room temperatures values of Oy to obtain L, and found 
it to vary from the mean value by at most about 20%. Neshpor (1959) 
has used the thermal characteristics (melting point, Debye temperature, . 
coefficient of linear expansion) to check the various relationships proposed 
between these quantities. He finds that semi-empirical functions fit 
the data better than those derived theoretically. For example it was found 
that 

12-45 T 473 (Cip23)12 
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(C=specific heat in cal/g, p=density, M=atomic or molecular weight) 
fits the data better than the Lindemann formula. No great accuracy can 
be expected from © values derived from the thermal properties of crystals 
but they are useful when other methods cannot be applied. 


® = 


3.8. Other Methods 


Some other methods of deriving © values have been discussed by 
Blackman (1955a). For ionic crystals infra-red data (‘Reststrahlen ’) 
and compressibilities can be used but these will not be discussed further 
as the emphasis in this paper is on metals. . 


3.9. Applications to Alloys and Compounds 


Discussion of experimental methods and results in this paper (§§ 6 and 7) 
has been restricted to cubic elements because inadequacies of the Debye 
approximation should affect these substances least. The application of the 
methods described earlier to disordered alloys and ordered compounds 
(or alloys) should be considered briefly. These applications are of interest 
because diffraction methods in particular have been widely used in the 
Soviet Union to determine the effects of various alloying additions on the 
Debye temperatures of metals; the mechanical properties of the alloys 
have then been related to the changes in the Debye temperatures and to the 
static deformations in the alloys (see, for example, Kurdyumov et al., 1955). 

The specific-heat and elastic-constants methods of determining ©) 
can be extended to cubic alloys and compounds without difficulties beyond 
those encountered with the elements, except that the assumption of a 
monatomic lattice no longer holds for ordered compounds. Debye 
temperatures can also be calculated formally for non-cubic materials from 
specific heats or elastic constants but the value of the resulting Op’s is 
limited because the assumptions of the Debye approximation are so far from 
being fulfilled. 

There are some differences when diffraction methods are applied for here 
one must distinguish between the vibration amplitudes of the individual 
components. For example, for a disordered binary alloy eqn. (6) must be 
rewritten (see e.g. Murakami 1953). 


in2 @ sin? @ 
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where m4, 1p are the atomic fractions of A and B respectively and By and 
B, their respective temperature factors. The values of B, and By could be. 
obtained by a least squares solution of the observational eqns. (16) but in 
practice F(hkl) and the atomic scattering factors are seldom known 
accurately enough for this to be done and the simplifying assumption is 
made that B,=B,=8. A further complicating factor for disordered 
alloys is the possible existence of a quasi-temperature factor due to static 
displacements of the atomic centres from the sites of the average lattice 
(Herbstein et al. 1956, Iveronovna et al. 1957, Borie 1957). If this is so, 
then use of atomic scattering factors to determine B will give values which 
are too large and hence values of ©, which are too small. Intensity 
measurements made at two temperatures are necessary for the elimination 
of this effect. For ordered crystals of stoichiometric composition it is 
always possible to separate the contribution to (Akl) from each kind of 
atom and hence obtain a separate value of B; for each atom. The Debye 
approximation is not strictly applicable to ordered crystals because the 
primitive unit cells are not monatomic but it is possible to derive a value 
for ©, if a suitable mean value for B can be found. Lonsdale (1948) has 


suggested that R 


where W,, Wz are the weight fractions of A, B in the ordered crystal. 
The experimental results summarized by Lonsdale (loc. cit.)are in reasonable 
agreement with this suggestion. 

Use of the Mossbauer effect promises at least partial alleviation of the 
difficulties encountered when diffraction methods are applied to disordered 
alloys. The recoilless emission (or absorption) of y-rays by an atom de- 
pends on the vibration spectrum of the atom in its actual environment in 
the host lattice. Thus Méssbauer measurements should make it possible 
to determine the vibration amplitude of the component emitting or absorb- 
ing the recoilless y-ray and thus provide information which cannot at 
present be obtained from diffraction measurements. 


§ 4. THEORETICAL DIFFERENCES BETWEEN DEBYE TEMPERATURES 
DETERMINED BY DIFFERENT EXPERIMENTAL METHODS 


In terms of general frequency distributions the formulae for the Debye— 
Waller factor (Blackman 1937) and the specific heat can be written: 


4h? cd ¥m iL 1e 1 ‘ 
Ba N. eae 
3NmkT > ih i”) E (5 as] 5) | id OE) 
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where p=hv/kT' and N,,(v) is the frequency distribution function for elastic 
waves of polarization j. Equations (17) and (18) have been written in the 
form shown in order to emphasize that a summation over the three types of 
waves with different polarization vectors must be carried out. The 
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integrands in (17) and (18) both consist of frequency distributions multi- 
plied by weighting factors, denoted by W, (hv[kT) and W, (hv/kT) respec- 
tively. These weighting factors are shown in figs. 3 and 4 as functions of 
frequency and temperature. It is at once apparent that different parts 
of the frequency distribution are quite differently weighted}, the 


Fig. 3 
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Variation of W, (hv/kT’) with v and J (see eqn. (17)). The curves of 
W, (hv[/kT)->co as v0 but the integrand of (17) is finite at v=0. 
The position of ©,, for Al (Debye approximation) is shown. Thus at 
300°x for example, the low-frequency part of the frequency spectrum 
(below about 30 x 10" c/s) will make a much larger relative contribution 

to B than the high frequency portion. 
i eee ee Se 
+ James (1950, pp. 223-5, especially fig. 75) has previously pointed out that 
the different frequency ranges of the elastic spectrum make very different 
contributions to B. James’ discussion is in terms of the Debye’ frequency 


spectrum. 
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low-frequency region being relatively much more important for the x-ray 
temperature factor than for the specific heat. Thus the effects of devia- 
tions from a parabolic spectrum on B factors will be different from the 
effects on specific heats. As the actual frequency spectrum deviates least 
from a parabolic spectrum at low frequencies, the Debye approximation 
will fit the x-ray temperature factors better than the specific heats ; this 
should be most marked at low temperatures but unfortunately no experi- 
mental data are available to test this prediction. As only the very lowest 


‘Fig. 4 
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Variation of Wy (hv/kT) with v and 7 (see eqn. (18)). This shows that at 
300°K for Al all regions of the frequency spectrum make approximately 
equal relative contributions to C,. 


lattice frequencies are excited in measurements of elastic constants, the 
Debye approximation will apply accurately to the @, values derived in 
this way, and deviations from the Debye approximation will not influence 
the shape of the @p (elastic constants) versus 7’ curve. 

Parabolic forms for NV,(v) must be inserted in (17) and (18) when the 
Debye approximation is made, i.e. N,(v) =(9N/v,,;°) v2 where v,,; will be 
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different for the transverse and longitudinal modes of vibration. This 
substitution gives 


ST Tien). eter 
Ba eel age or | (19) 
D 2D 
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where ©, and ©, are Debye temperatures appropriate to longitudinal and 
transverse vibrations respectively 
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s(x) has been defined in eqn. (9) and 
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is the familiar Debye integral from specific-heat theory. Except at low 
temperatures the difference between %(x,) and u(x,) will be small; if an 
average value of ©,, be chosen such that 
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then (19) will become identical with (8). Similar considerations apply to 
(20) ; if an average value of @,, be chosen such that 


BO Ree Ot 00ers (99) 


then (20) will become identical with (2). Zener and Bilinsky (1936) were 
the first to point out that the average Debye temperatures in (21) and (22) 
were defined differently and that their ratio could be expressed as a function 
of Poisson’s ratio, using the expressions for ©, and @p derived from 
elastic constants (§3.3.1.). Their values of © ,/@p are reproduced in 
table 1. In this paper all numerical values of Debye temperature refer to 
@p, unless specifically stated otherwise; where necessary, values originally 
referring to ©, have been converted to Op, using table 1. The differences 
between @,, and @, never exceed 7% and are generally about 3-4%. 


Table 1. ©,,/@p as a function of Poisson’s ratio o (from Zener and 
Bilinsky 1936) 


o 0-25 0-3 0-35 0-4 0-45 0-5 


Oy/Op | 1-025 1-027 1-032 1-040 1-051 1-070 


Differences in the methods of averaging over the velocities of the elastic 
waves also lead to differences between the values of 2/, the ordinary 
Debye-Waller factor, and 2p, the corresponding factor for temperature- 
diffuse x-ray scattering (see §3.4.3.), with the effect that ©, derived from 
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Debye temperatures (Op) of f.c.c. metals 
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Elastic constants 


Young’s 


Diffraction measurements 
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Ses Temp.| Entropy | Specific heat 
Ag 0 = 926-0?) 
300 207 | 2168.4 
AL | 0 | —_ | 427-7) 
300 379M | 396%) 
380 
Au | 0 | — | 164-7 
300 168” 185(3, 4) 
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} TDS under Bragg peaks taken into account. 
{ Extrapolated value. 


§ Value of 413-437 °K derived from Méssbauer effect measurements 2), 
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Bragg intensity measurements would be expected to be lower than that 
derived from measurements of TDS (Chipman and Paskin 1959 a). Experi- 
mental values for Cu and Pb by the two methods (table 2) were in agreement 
with this expectation but an approximate quantitative treatment suggested 
that the velocity effects were too small to account for the differences actually 
found}. Uncertainty as to the theoretical significance of @,, from TDS 
measurements is an added reason for preferring the other diffraction 
methods of determining ©,,. 

One further consequence of the deviation of the true vibration spectrum 
from the parabolic form must be mentioned here. The limiting form of the 
true spectrum at low frequencies is parabolic (Blackman 1955a, p. 340); 
the value of @p derived from measurements of elastic constants will be 
that appropriate to this part of the spectrum. However the x-ray value 
of Op (i.e. after the correction for the difference between ©, and @, has 
been applied) will be less than the elastic-constant value if the real frequency 
spectrum rises above its initial parabolic shape. This is because the para- 
bola that fits the true spectrum best over an appreciable part of the total 
frequency range will be steeper, and thus have a smaller value of v,,, then 
the parabola chosen from the limiting shape of N(v) at the very lowest 
values of v. The rule @p(x-ray) <.@p (elastic constants) is a very general 
one as most real frequency spectra have ‘humps’ in the lower frequency 
range (Blackman 1955b). It is difficult to make a similar prediction for 
@p (specific heat)eas the higher frequencies play a much more important 
role in determining C,, than B. 


§ 5. VARIATION OF DEBYE TEMPERATURE WITH TEMPERATURE 


The formal derivations of eqns. (16) and (17) do not require that either 
the frequency spectrum or the maximum frequency should be invariant 
with temperature (see Born and Huang 1954, chapter IT). As @p depends 
explicitly on the volume (e.g. Fowler 1955), it is to be expected that it will 
in turn vary with temperature. Using this volume-dependence of 0 
Paskin (1957) has shown that the temperature-dependence of © is given by 


OT) aol 
sere = | 7a | Spoah ae ieee TCLS) 


where 1(7') is the linear dimension of the crystal at temperature T and 
(7) that at 7’), and y is Griineisen’s constant (= VB/C,x where Bis the 
coefficient of volume expansion and x the compressibility). This formula 
is closely equivalent to an earlier formula due to Zener and Bilinsky (1936) 
which, however, contains less accessible quantities than Paskin’s formula. 
A rather similar formula was given earlier by Owen and Williams (1947), 
who based their derivation on the Einstein form of the frequency spectrum. 
The predictions of Paskin’s formula are compared with the available 


experimental data in $7.2. 


+ See note added in proof. 
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§ 6. THE EXPERIMENTAL RESULTS 


Values of Debye temperatures at 0°K and 300°K determined by various 
experimental methods are given in tables 2 (face-centred cubic metals), 
3 (body-centred cubic metals), and 4 (diamond-cubic crystals). 

The aim in preparing these tables has been to present a comprehensive 
list of values for those substances where more than one experimental method 
has been used. In general only values from the most recent references are 
given, and some attempt has been made to select the most reliable results 
where a number of determinations by the same method are available. 
Discussions of comparable results are usually given in the references cited 
and have not been duplicated here. 

The results listed are all in terms of ©, (see §4); where the value listed 
in the original paper was in terms of © y,, it has been recalculated here using 
table 1 and values of Poisson’s ratio given by Koster (1954). In many 


Table 3. Debye temperatures ©, of b.c.c. metals 


Elastic constants Diffraction measurements 
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instances the value of @, given in the tables has had to be calculated here 
from experimental data such as Young’s modulus or single-crystal elastic 
constants; the reference given is then to the experimental data used. 
A number of the values in the tables differ appreciably from those given by 
Furukawa and Douglas (1957) in the Handbook of the American Institute of 
Physics. : 

The necessity for correcting x-ray determinations of Debye temperature 
for peaking of temperature-diffuse scattering under the Bragg peaks has 
been mentioned earlier (§3.4.1.); this can only be done when the experi- 
mental conditions are known and only those values in the tables marked with 
an asterisk have been corrected in this way. This correction is not more 
than a few per cent for @,, greater than 300°K, and will always reduce ©, 
below the uncorrected value. 

For a number of substances (Al, Pb, Cu, Ag, and Au) it is possible to 
extend these comparisons of different ©, values to cover a wide range of 
temperatures and these results are given graphically in figs. 5, 6, 7, 8 and 9. 


Table 4. Debye temperatures ©, of Diamond-type cubic crystals 


Elastic constants Diffraction measurements 

Temp.| Entropy | Specific heat Tague 
f-curve different 

temperatures 


Young’s 
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§ 7. Discussion OF RESULTS 
7.1. Op Values at 0° and 300°K 
It is convenient to discuss the results obtained from related types of 
experimental measurements first and then to compare the tabulated values 


by structure type. 
The @, values derived from the elastic constants of single crystals, or 
Young’s modulus (Z) of polycrystalline samples, can only be satisfactorily 
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compared at 300 °K because of the lack of experimental values of EF at very 
low temperatures (curves of H against temperature have been extrapolated 
to 0 °K by Koster (1948 a) but it has been felt that these values could not be 
considered as true experimental determinations and so comparisons have 
not been made at 0 °K). In general there is agreement between the two sets 
of values to within 1-2%, with rather larger deviations for Ni (5%) and 
Mo (4%). Differences are found for Pb (20%), Si (22%) and Ge (25%) and 
either the ¢,; or H values must be wrong. For Si and Ge there does not 
seem to be any possibility of error in the c,; values because of agreement 
between different authors (McSkimin 1953, Fine 1955) and it seems most 
likely that the values given for Young’s modulus are unreliable. The 
various values given for Li are also very scattered. 

The principal x-ray methods which have been used in practice are 
(a) evaluation of @,, from intensities measured at a given temperature by 
using a theoretical f-curve, (b) ©, from comparison of intensities measured 
at two or more temperatures. Both methods have been applied to Cu, Ni, 
Pb, Al, o-Fe, and Si and agreement is good for all these elements except 
Si. Too few determinations of ©, from temperature diffuse scattering 
alone, or from a combination of TDS and Bragg intensities (Borie’s (1956) 
method), have been made for the capabilities of these methods to be 
reliably assessed but one would expect ©, from TDS alone to be less 
accurate than the other methods because of the relative weakness of the 
TDS, the need for correction for Compton scattering and the approximations 
in the theory. Borie’s method is very attractive in principle but the value 
given for Cu is not very accurate because the TDS had to be treated on the 
basis of the approximate theory available at the time the work was doney. 
The different results for silicon illustrate two important points in the 
application of the x-ray methods: the value of 758°K was obtained from a 
comparison of intensities at two temperatures, but this comparison was 
restricted to one reflection and hence it will be very susceptible to experi- 
mental errors; this value cannot therefore be considered very reliable. 
Secondly, two different values were obtained by fitting different f-curves to 
the same set of very accurate intensity measurements (Géttlicher et al. 
1959, Gottlicher and Wo6lfel 1959). The f-curves used in the first of these 
references was that due to Viervoll and Ogrim (1949), where electron 
exchange was not taken into account, while in the second paper the more 
recent calculations of Freeman (1959), which include the effect of exchange, 
were used. 

The determination of values of ©, from specific heats are beset by two 
main difficulties. At very low temperatures (7'~0°xK) the correction for 
the electronic specific heat is important and this fact, added to the con- 
siderable experimental difficulties of the measurements, makes many 
reported values unreliable. At higher temperatures (7'’>@,) the specific 
heat tends towards the equipartition value of 6cal/gatom and is not 


} This remark is incorreet—see note added in proof. 


Methods of Measuring Debye Temperatures 339 


sensitively dependent on the value of @p. There are also theoretical 
difficulties in interpreting specific heats at relatively high temperatures and 
so this method is only really useful for 7< ~4@,. Although values given 
in the tables for ©) from specific heats refer nominally to 300°x, in fact the 
values were taken from curves of @p against 7’ at 7’=1@, (see $3.1). 
Many of these values were taken from the curves given by De Sorbo (1954). 
There are often fairly large differences between ©, values at 0°K and those 
at higher temperatures. This is in qualitative agreement with theoretical 
calculations based on the Born—von Karman lattice model, which show 
that @p will be greatest at 0°xK. 

One important practical difference between specific-heat and elastic- 
constant values of © on the one hand, and x-ray values, on the other, is 
that the former are determined at small temperature intervals and thus a 
detailed curve of @p against temperature can be obtained ; the shape of this 
curve at low temperatures (for specific-heat values of ©p) is particularly 
important for the comparison of experiment with theory. The x-ray 
values, at least in the past, have been determined at rather wide tempera- 
ture intervals and further only a few measurements} have so far been 
reported at temperatures below 80°K. 

Strictly speaking, the ©p values calculated from the vibrational entropy 
at 298°K (Sy9g) do not belong in the tables, as they do not refer to a particular 
temperature of measurement. This is because the entropy is obtained as 
an integral over the whole temperature range 


298 (G@ 
geen Fat). 
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It will be noticed that the ©, (entropy) values are usually lower than the 
values obtained from other experimental methods; this is a consequence 
of the dip in the curve of ©, against 7 at low temperatures which is caused 
by the discontinuous, lattice nature of the solid. ©) values from other 
methods have not been included in the tables because they do not refer | 
to definite temperatures or because the theoretical position is not clear. 

It would be expected that results from the solidified rare gases would 
be most useful for comparison with theory but experimental measurements 
are difficult and the only element for which two roughly comparable 
results are available is Ne. Here specific-heat measurements (Dobbs and 
Jones 1957) give © = 64°K at T'~ 60 °K, while neutron-diffraction intensity 
measurements of Bragg peaks at 4°K give @y=73°K (Henshaw 1958). 
The Bragg intensities were not corrected for peaking of TDS. It would, 
of course, be expected that the Debye temperature near absolute zero would 
be somewhat larger than at higher temperatures. 


Par ea A ee ee 

+ Of these, that for Al at 5 °K (Walker 1956 b) was made as a check only and is 
not very accurate (private communication from Professor Walker). Renninger 
(1952) has made intensity measurements on NaCl at 20 °K but these results were 
not corrected for the peaking of the TDS beneath the Bragg peaks (Nilsson 
1957). 
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A careful comparison has recently been made of elastic-constant, (cy) and 
specific-heat values of @p at 0°K (Alers and Neighbours 1959) and it has 
been shown that excellent agreement between the two values is obtained 
for Cu, Ag and Au. This is also true for Al (Phillips 1959) and Pd (Rayne 
1960). The specific-heat values for Ni and Th are less accurate than those 
for Cu, Ag, Au, Al and Pd but it seems that agreement for these metals is 
also likely. Good agreement is also obtained for diamond, Si and Ge. 
All this is in agreement with theory. The agreement for Liis worse but this 
is probably a consequence of the occurrence of a martensitic transformation 
(Barrett 1956) at low temperatures. Because of this transformation 
neither specific heats nor elastic constants of the body-centred cubic form 
can be measured below about 60°K (Martin 1960 b, Nash and Smith 1959). 
Extrapolation of the elastic constants to 0°K is not difficult but this cannot 
be done satisfactorily for the specific heat. Thus a comparison of the two 
@y values at 0°K is of doubtful validity. No other suitable experimental 
values have been determined for the substances included in the tables. 

The most extensive sets of comparable experimental data available at 
300°K are those for Cu, Al and Pb and the agreement is generally good. 
The elastic-constant values of ©p do, however, lie consistently higher than 
the x-ray results by a few per cent. As mentioned earlier (§ 4) the relation 
©, (elastic)/@) (x-ray) >1 is of a very general nature and is due to the 
existence of peaks in the vibrational spectrum at low frequencies. The only 
exceptions to this rule noticed in the tables are Cr and Li; in both instances 
confirmation of the experimental measurements is desirable before a differ- 
ence in behaviour can he considered as established. The agreement between 
elastic-constant and some x-ray values of @, for the other face-centred 
cubic metals is not good: for example for Ni the difference is about 25% 
(using the values given by Iveronova et al. 1951, and Ilyina and Kritskaya 
1955 a) and for Pt about 30%. In some instances this may be due to 
systematic errors in the experimental results. The x-ray methods seem 
to be particularly liable to experimental errors and one hesitates to accept 
an X-ray value for @p which has not been independently confirmed. Thisis 
illustrated by Ni for which three comparable diffraction measurements of 
©, have been reported. In the most recent of these experiments 
(Zhuravlev and Katsnel’son 1959) accurate intensity measurements were 
made at 215°C, 20°C and —175°c with Mo Ka radiation and a value of Oy 
derived which is in much better agreement with the elastic-constant value 
than the results reported in the earlier work. It is interesting to note that 
for Ni @p (elastic constant) > ©, (x-ray), in agreement with the results 
for other elements. 

The results for the three diamond-cubic elements (diamond, Si and Ge) 
provide an interesting pattern of discrepancies and incidentally show how 
much the accepted experimental values of ©, for these elements have 
changed over the past few years. For diamond the earlier measurements 
of specific heats seemed to show ©, ~ 1900°K at 0°K (Pitzer 1938, De Sorbo 
1953) but more recent measurements by Desnoyers and Morrison (1958) 
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and Burk and Friedberg (1958) agree, within combined experimental 
error, to give Op (0°K)=2232°K. Single-crystal elastic constants have 
only been measured at room temperature but the three sets of values 
reported, and the Debye temperatures calculated from them (see table 4) 
differ considerably. ©,, calculated from the latest set of values (McSkimin 
and Bond 1957) is in excellent agreement with the recent specific-heat 
valuesj. A value of ©, (=1990+80°K) has also been obtained from a 
comparison of x-ray diffraction intensities at 125 °K and 300 °K by Post 
and Miksic (Post 1960), and this is in good agreement with what would be 
expected from the specific-heat results. However the values of Oy 
deduced from measurements of intensities at room temperature, using 
available atomic scattering factors, do not agree well with those given 
above. The reason for this discrepancy is not clear as two sets of very 
accurate intensity measurements have been made on diamond (Brill et al. 
1939, Gottlicher and Wolfel 1959) and the results are concordant. Further, 
the atomic scattering factors calculated by various authors (Hoerni and 
Ibers 1954, Berghuis ef al. 1955, McWeeny 1951, Freeman 1959) are also in 
good agreement. _In their analysis of the experimental results Géttlicher 
and Wolfel, and Carpenter (1960) (using the measurements by Brill et al.), 
used intensities in the regions 


0-7< = <1:3A-1 and a >1A4 

respectively so that uncertainties due to the valence-electron distributions 
were avoided; the small difference in their results is probably due to the 
slightly different experimental data used. Correction for peaking of TDS 
under the Bragg peaks would be expected to be small and one may conclude 
that a value of ©,,=1520+60°K (where approximate 95% confidence 
levels are given) is well established from analysis of x-ray intensities at 
room temperature with available atomic scattering factors. 

Rather similar effects are found for Si. The specific-heat and elastic- 
constants measurements give @, values which agree excellently at 0°K and 
300°x. Accurate measurements of x-ray intensities have been made at 
300°K but the difficulty of obtaining an unambiguous x-ray value for Op of 
Si using atomic scattering factors has been mentioned earlier. IfFreeman’s 
(1959) f-curve is accepted because electron exchange was taken into account, 
then there is a discrepancy of about 20% between the two sets of Op values 
at 300°k. Presumably the differences found among the various © values 
for diamond and silicon can be accounted for in terms of the real frequency 
distributions of these elements, but no detailed treatment has been made. 
It should also be remembered that the diamond lattice is not monatomic. 
No x-ray measurements have yet been reported for Ge. 

The position with regard to a—Fe is as follows. Two different x-ray 
determinations, one using the f-curve (Herbstein and Smuts 1960) and 


Aner ees we ee eee 
+ Because of the high @p of diamond, the elastic constants will hardly change 
between 300 and 0°K. 
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the other using intensities measured at 300°K and 80°K (Ilyina and Kritskaya 
1955 b) are in satisfactory agreement and both sets of measurements were 
made under suitable experimental conditions. Although the single- 
crystal elastic constants (see Hearmon (1946) for summary) were measured 
more than 20 years ago, a number of independent measurements are in 
good agreement; in addition the values of Op calculated from Young’s 
modulus (Koster 1954) and the ¢, agree well. Thus «—Fe appears to be 
another example where @p (elastic constants) > Op (X-ray). 

A rough picture-of the extent of the agreement among the various values 
of @, reported may be obtained from the fact that for about two-thirds of 
the examples given in the tables, the deviation of the extreme value from 
the mean is greater than 10%. Satisfactory agreement (deviation of 
extreme from mean less than 5°) has been reported only for Cu and Al. 

Fig. 5 
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Variation of @) with 7' for Al. The references to the values given in the figures 
are: Sutton 1953, Koster 1948 a, Paskin theory 1957; JBW, James, 
Brindley and Wood 1929; BWW, Bensch, Witte and Wéolfel 1955: 
Chipman 1960; OW, Owen and Williams 1947, 


7.2. Variation with Temperature 


The most extensive results for the variation of @, with temperature 
have been given for Al, where values of single-crystal elastic constants 
(Sutton 1953), Young’s modulus (Koster 1948 a), specific-heats (De Sorbo 
1954) and x-ray intensities (James et al., 1929, Owen and Williams 1947 
and Chipman 1960) between 80°K and about 900°K are available. The @ 
values calculated from these various experimental measurements are ee 
in fig. 5, together with the curve given by application of Paskin’s (1957) 
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method. The two sets of elastic constants give @p values in good agreement 
with each other up to about 700°K while the x-ray results of Owen and 
Williams and Chipman are also in reasonably good agreement, with 
Chipman’s results scattering considerably less than those of Owen and 
Williams. At first sight this agreement is rather surprising as Chipman 
corrected his results for TDS peaking under the Bragg peaks whereas 
Owen and Williams’ work was done long before the importance of this 
correction was appreciated. It is possible that their method of intensity 
measurement was such that the TDS was largely corrected for uninten- 
tionally but sufficient details were not given in their paper to judge if this 
was indeed so. 

It should be noted that the @, values from elastic constants are obtained 
in absolute terms but that the x-ray values determined from the variation 
of Bragg intensities with temperature will vary according to the value 
assumed for © at some standard temperature, e.g. 300°K. Chipman’s 
values were calculated for 0, =382°K at 300°K, a value which is somewhat 
lower than the corresponding values obtained using the atomic scattering 
factor curve for Al. It is therefore possible that the x-ray curve should be 
shifted upwards by about 13°K, approximately parallel to its present 
position. The x-ray and elastic-constant values of @p show the same 
trend with temperature up to about 700°K but above this temperature the 
values from Young’s modulus fall off much more sharply than the x-ray 
values. It would be desirable to check the variation of Young’s modulus 
with temperature above 700°K. The fact that the temperature variation 
of ©, could not be accounted for in terms of volume changes alone was 
noted by Owen and Williams and the contrast with their results for Cu and 
Auemphasized. There is a marked difference in slope between the experi- 
mental values and the curve obtained from Paskin’s formula. 

Fairly extensive data are also available for lead (fig. 6). The absolute 
difference between the various curves is quite large. The discrepancy 
between the single-crystal and polycrystal elastic-constant ©p values has 
already been noted. The x-ray results (here fixed along the @p axis by the 
good agreement between the two x-ray results at 300°K) are closer to the 
Young’s modulus values. Volume changes alone are not sufficient to 
account for the temperature variations shown. 

Rather less detailed results are available for Cu (fig. 7). The x-ray value 
at 300°K has been determined independently by two groups of workers 
(Coyle and Gale 1954+, Chipman and Paskin 1959 a) and the agreement is 
good. This value can be used in turn to obtain values at 80°K from the 
measurements of Chipman and Paskin, and above room temperature from 
the measurements of Owen and Williams (1947). The agreement between 


elastic-constant and x-ray values of @, is not very good. In particular 
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+ The original data have been analysed by the author, using newer values of 
the atomic scattering factors and taking iato account peaking of TDS beneath 
Bragg peaks. The new value of 4, is oaly very slightly dferent from that 
given by Coyle and Gale. 
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Variation of @p with 7 for Pb. The references are given in the caption to 
fig. 5 except for Waldorf 1960. 
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Variation of ©, with 7' for Cu. The references to the values given in the figure 
are: OG, Overton and Gaffney 1955; CP, Chipman and Paskin 1959 a; 
CG, Coyle and Gale 1954 (also see caption to fig. 5). 
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the elastic-constant values show a much stronger temperature dependence 
than the x-ray values, and one may tentatively suggest that this is because 
the x-ray values of Owen and Williams were not corrected for the TDS 
under the Bragg peaks. Newx-ray measurements above room temperature 
would be desirable. 

Silver has been investigated with x-rays by four groups of workers 
(Andriessen 1935, Boskovits e¢ al. 1958, Spreadborough and Christian 
1959, Haworth 1960) and the results reported vary considerably among 
themselves and also differ appreciably from the @, values derived from 
measurement of elastic constants. The x-ray results have been re-inter- 
preted (Herbstein 1961) using Chipman’s method (1960) and the agreement 
with the other values (fig. 8) is good. Volume changes alone are not 
sufficient to account for the temperature variation of @>p. 
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Variation of @y with 7 for Ag. The references to the values given in the figure 
are: NA, Neighbours and Alers 1958; H, Herbstein 1961 (also see 
caption to fig. 5). 


X-ray results for Au have only been obtained by Owen and Williams 
(fig. 9) and the temperature variation of @p can be explained in terms of the 
effect of volume changes. Again the agreement with the values of @p 
deduced from the temperature variation of Young’s modulus is not very 

ood. 
B Other experimental results on the temperature variation of Op have been 
reported for lithium (Alexopoulos 1932, Pankow 1936), sodium (Dawton 


2C2 
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1937, 1938) and platinum. ‘The results for the first of these elements show 
fair agreement between the various methods, while for sodium the agreement 
is worse. For platinum there is marked disagreement between the 
x-ray results (Alexopoulos and Euthymiou 1954) and those calculated 
from Young’s modulus. The examples cited in the tables indicate the 
quality of the agreement for these three metals and detailed diagrams have 
not been given. 
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Variation of @p with 7’ for Au. The references to the values given in the figure 
are in the captions to figs. 5 and 8. 


Electron-diffraction measurements of vibration amplitudes over a range 
of temperatures have been reported for Au (Coster and van Zouten 1939, 
Menzel-Kopp and Menzel 1955), Ni (Yamaguchi and Katsurai 1948), 
Cu (Menzel-Kopp and Menzel 1955, 1956) and Ag (Menzel-Kopp and 
Menzel 1956). The Debye temperatures derived for Au (Coster and van 
Zouten) agree reasonably well with the specific-heat values. Menzel-Kopp 
and Menzel made measurements of vibration amplitudes in Cu and Ag by 
both reflection and transmission methods. The accuracy of these measure- 
ments is not high enough to warrant comparison with the values shown in 
figs. 7 and 8 but the results are nevertheless of considerable interest because 
they show that the vibration amplitudes in the surface of the metal 
(measured by reflection electron diffraction) are considerably larger than 
those in the bulk of the metal (measured by transmission electron diffrac- 
tion). ‘The approximation of simple harmonic vibrations is fairly accurate 
for the bulk of the metal, even near the melting point, but breaks down 
badly in the surface layers. Rather similar results were reported for Ni. 
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§ 8. CoNCLUSIONS 


Of the methods available for determining Debye temperatures, that 
based on single-crystal elastic constants appears to give the most reliable 
results as independent measurements are generally, but not always (cf. 
diamond) in good agreement. Use of Young’s modulus is nearly as 
satisfactory except where difficulties due to large grain size or preferred 
orientation make interpretation hazardous. Elastic-constant measure- 
ments can be conveniently made over the whole temperature range from 
liquid -helium temperature up to the melting point and at small temperature 
intervals. They do suffer from the serious disadvantage that only the low- 
frequency region of the vibration spectrum is involved in the measurement 
and consequently the curve of ©, against 7’ will not provide the information 
needed for testing frequency spectra derived by more elaborate theoretical 
or experimental techniques. This information can be obtained from 
specific-heat measurements, but these will not give accurate @, values at 
temperatures much above $0). The x-ray diffraction methods have been 
used fairly widely but difficulties of interpretation remain; for measure- 
ments at a given temperature a theoretical atomic scattering factor must be 
used. to derive ©,, while measurements over a range of temperatures do 
not give unambiguous results. The first of these difficulties does not 
exist for neutron diffraction. 

The careful comparison made by Alers and Neighbours (1959) shows that 
good agreement exists between elastic-constant and specific-heat values 
of ©» at 0°K as required by theory. When this comparison is extended to 
all available methods over a wide range of temperatures one finds that 
satisfactory agreement exists only for Al. It does seem that many of the 
discrepancies encountered may be due to errors of measurement and 
experimental verification of these discrepancies is desirable before attempts 
are made to find theoretical reasons for their existence. Nevertheless 
some large discrepancies between the results of different methods do appear 
to be established, e.g. for diamond, and there also appears to be definite 
evidence, shown most strikingly for Cu and Pb, that 


Op, (elastic constants) > Op (X-ray). , 


This is in accordance with theoretical expectations (Blackman 1955b). 

For Ag, Al and Pb, the temperature dependence of @, derived from 
x-ray and elastic-constant measurements is in substantial agreement and 
cannot be explained by assuming that the changes in @p are due to volume 
changes only. However, for the other elements for which results are 
available (Au, Cu), quite different temperature dependences are shown by 
Op (x-ray) and Op (Young’s modulus) ; the former are in agreement with 
the theory of the pure volume-dependence of @p while the latter are not. 
The situation is complicated by the fact that the x-ray measurements were 
not corrected for the peaking of TDS under the Bragg peaks and it is 
difficult to estimate the influence of this on the results (see discussion for Al). 
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At low temperatures no results are available for the variation of 
@, (x-ray) with temperature, nor has there been any theoretical discussion 
of this point. When the influence of the weighting factors W, (hv/kT) and 
W, (hv/kT’) (see figs. 3 and 4) is taken into account, it would be expected that 
@, (x-ray) would be less affected than Op (specific-heat) by deviations of 
the actual vibration spectrum from the Debye form. 

The present survey of experimental results has been restricted to elements 
because the best agreement among the results of the various experimental 
methods would be expected for the elements. The discrepancies found 
have been such that it has not been felt worth while to extend the compari- 
son to alloys, both disordered and ordered, or to salts, because it would be 
difficult at this stage to decide which effects were due to experimental 
errors, which to deviations from the Debye approximation and which to the 
effects of inhomogeneities or long-range order in the materials. Neverthe- 
less it should be remembered that many of the important potential applica- 
tions of Debye temperatures are to alloy systems for the calculation of such 
quantities as thermal entropies and x-ray temperature-diffuse scattering. 

Despite the 49 years that have elapsed since Debye first proposed his 
approximate form of the lattice vibration spectrum, its usefulness has not 
been exhausted, nor indeed is our present knowledge of the Debye tempera- 
ture, approximate parameter though it may be, adequate for many of the 
applications required. Many interesting problems remain to be solved 
and many useful measurements must still be made. 
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@y values from specific heats, although referred nominally to 300°K, are 
actually taken from graphs of ©, against 7’ at 7 ~3@p. The only exception 
to this rule is diamond (table 4) where @p (300) does refer to 300 °K. 


Note added in proof.—Since this paper was submitted (February, 1961) a 
number of relevant papers have appeared and are discussed below under the 
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appropriate section headings. The literature available in Pretoria to the 
end of July, 1961 has been covered. 

3.3.1. A new summary of elastic constants of the elements (for poly- 
crystalline samples) has been given by Koster and Franz (1961). The 
values for Young’s modulus are generally the same as those quoted in 
Koster’s earlier paper (1954). The only noteworthy change in ©p value 
is for germanium (see table 4). The new value of ©p from Young’s modulus 
at 300°K is 313°K, thus reducing but not eliminating the difference between 
Op(#) and 05(c¢;;). 

3.4.3. Borie (1961) has shown that the values given by Paskin (1958, 
1959) for C, and C, are not correct and that 2M p= 2M when the Debye 
model is used. Borie’s accurate expression for the temperature-diffuse 
scattering from a monatomic cubic powder is 
Ips = Nf? [1—exp (— 2M)]+ Nf2exp (—2M)(C,—1)(2M+M?)... (12a) 
where the symbols have the same significance as before. Oj, is the function 
of sin@/A calculated by Warren (1953) for 7'>© and by Herbstein and 
Averbach (1955) for all 7’. . 

Another method of determining Debye temperatures from diffuse 
X-ray scattering, first proposed by Cartz (1954), has recently been dis- 
cussed by Canut and Amorés (1961). The intensity of diffuse scattering at 
a particular point in reciprocal space depends on the temperature of the 
crystal; in general the intensity of /th order diffuse scattering first increases 
as the temperature is raised and then falls. Thus there is a measurable 
inversion temperature 7',,, given by 


fe Sm Ody . . . . . . . (12 b) 


d,,, in this expression corresponds to the reciprocal-lattice position at 
which the measurement was made. Canut and Amords have used some of 
Cartz’s experimental measurements of 7',,, for lead to derive @~ 90°K. 
This method, which is' valid for 7’> © only, should be useful in obtaining 
approximate values of © for crystals, but the exact significance of the 
value obtained and its relationship to the other types of @ discussed in this 
review, will depend to a considerable extent on the exact conditions of 
measurement. 

3.5.1. Three further reviews of Méssbauer effect have appeared 
(Belozerskii and Nemilov 1961, Shapiro 1961, Méssbauer 1961) and the 
first two cover the thermal-vibration aspect in some detail. Some of the 
precautions necessary when comparing values of © from Méssbauer 

.experiments with those obtained by other methods are discussed in these 
papers. ‘Two articles have appeared on the Méssbauer effect in ™9Sn. 
In the first of these (Boyle et al. 1960) the frequency shift due to a tempera- 
ture difference between source and absorber was measured and could be 
accounted for in terms of the theories of Josephson (1960) and Pound and 
Rebka (1960a). In the second paper (Boyle et al. 1961) the intensity of the 
recoilless resonance absorption in tin was measured from 120°K to the 
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melting point. At lower temperatures @ ~ 142°K but at high temperatures 
the results indicated considerable anharmonicity of the lattice vibrations. 

3.9. A review has appeared on “The Physics of the Strength of 
Crystalline Materials” (Garber and Gindin 1960), which includes a dis- 
cussion of Soviet work on thermal-vibration amplitudes in alloys as 
measures of the strength of interatomic bonding. 

6. New measurements have been reported of the single-crystal elastic 
constants of w-iron (Rayne and Chandrasekhar 1961) over the temperature 
range 4-2-200°K. ‘The new values agree well with the earlier values for 
200°K (Hearmon 1946). The new ¢,; values give @, = 467°K (at 300°K) com- 
pared to 464°K (see table 3).. A value of @p =477 + 2°K at 0°K was calcu- 
lated from the elastic. constants by Rayne and Chandrasekhar. Specific 
heat measurements at very low temperatures (Cheng et al. 1960) had been 
interpreted to give @p =445°K at 0°K, but Rayne and Chandrasekhar point 
out that the spin wave contribution to the specific heat was neglected in that 
analysis. A reasonable value for the spin wave contribution was obtained 
on the basis of 0p =477°K. 

The Debye-Waller factor for copper has been determined from single- 
crystal x-ray intensity measurements over the temperature range 
4-2-500°K (Flinn et al. 1961). The results could be interpreted satisfac- 
torily in terms of the frequency spectrum calculated from the central-force 
model with first and second-neighbour interactions. This paper has been 
seen in abstract only and it has not been possible to compare the new results 
with those summarized in fig. 7. 
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§ 1. INTRODUCTION 


THE purpose of this article is to present a review of the theory of the 
interactions of electron and lattice vibrations (phonons) in metals. The 
only recent reviews in this field are those of Bardeen (1951, 1956). These 
were largely confined to the special problems that arise in the theory of 
superconductivity. Our aim, however, will be to present a general review 
of the field and we shall discuss the problems that arise in both the 
equilibrium and non-equilibrium properties of a metal. In general we 
shall not be concerned with the detailed results that emerge when any of , 
these are applied to particular metals. This field has been reviewed in 
detail by Ziman (1960, 1961). Rather, we shall be concerned with the 
question of the validity of the methods and the general results that follow 
from them. We must emphasize that the attitude we take is that 
fundamentally the so-called electron—phonon problem is one particular 
aspect of a many-body problem. The many-body problem being that of 
a number of Fermions (electrons) interacting with the ions that constitute 
the lattice. If we take this point of view, then not only do we have a 
many-body problem to deal with but the basic interaction between the 
particles will be a long-range Coulomb interaction. As is well known, 
the treatment of such interactions is much more difficult than the 
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short-range interactions that are of importance in electronically saturated 
systems. We shall see that the Coulomb forces are vital in understanding 
both the motion of the ions in the lattice and the effective coupling between 
the moving ions and the electrons. 

In §2 we shall give a qualitative discussion of the physical and 
theoretical problems in the field. We do this because much of what 
follows in §3 is rather technical. This technical discussion is necessary 
because, if we wish to discuss the validity of the methods that have been 
used then we must present a reasonably detailed outline of them. Our 
aim in §2 will be to avoid the technicalities of the problem and to give a 
clear physical picture of what the problems are and to describe in 
qualitative terms the methods that are discussed in §3. It is hoped that 
the reader who is not interested in the details of the methods will be able 
to get a reasonable idea of what they involve by reading this section. 

Section 3 is devoted to a fairly detailed exposition of some of the methods 
that have been used to discuss the problem of electron—phonon coupling 
in metals. We have been deliberately selective and have chosen for 
discussion only those methods which seem to us to be important. In fact 
most of the methods are both in principle and application rather closely 
related to one another. For this reason the major part of §3 is devoted 
to discussing the method used by Nakajima (1953). A much shorter 
discussion is given of the Born—-Oppenheimer (Chester and Houghton 
1959) and Hartree—Fock (Toya 1952) methods. The collective coordinate 
approach (Bardeen and Pines 1955) is only briefly mentioned as a clear 
description of this method is available in the review article by Bardeen 
(1956). 

In §4 we discuss the application of the methods reviewed in §3 to the 
calculation of thermodynamic functions and transport coefficients. As we 
have already emphasized, we shall be only interested in the general 
qualitative features that emerge. We conclude this section by pointing 
out that for temperatures at or above the Debye temperature of the metal 
it should be possible to approach the electron-phonon problem by 
expanding the quantities of interest as a power series in 6/7’, where @ is 
the Debye temperature of the metal. 

We finally discuss in § 5 the difficulties that arise when we try to apply 
any of the methods to the problem of superconductivity. Our main aim 
will be to give a critical discussion of the validity of the standard methods 
when applied to this problem. We shall briefly discuss the work of 
Bogolyubov (Bogolyubov e¢ al. 1959) in this field. We conclude this 
section by pointing out that there are possibly several quite different 
schemes of approximation available for this particular problem. None 
of these has as yet been explored. 

We now give a brief historical survey of the subject. The first funda- 
mental paper is really that of Bardeen (1937). In this paper he set up the 
Hartree equations for the metal allowing for the motion of the unscreened 
ions. Solving these to first order in the ion displacements he showed that 
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the electrons were able to screen the bare Coulomb interaction with the 
vibrating ions. This self-consistent field method was extended by Toya 
(1952) to calculate the normal mode frequencies of a metal. The next 
important step was taken by Frohlich (1950, 1951, 1952). He pointed 
out (a) that electron-phonon interactions were likely to be the origin of 
superconductivity—this conjecture was confirmed at the same time by 
the discovery of the isotope effect (Maxwell 1950a,b), and (b) that the 
electron—phonon interaction could lead to a phonon self-energy term 
which would ‘renormalize’ the lattice frequencies. Fréhlich’s work was 
based on a phenomenological screened interaction and ignored the 
long-range Coulomb effects. A similar approach was employed by 
Nakajima (1953). He started from an unscreened interaction and 
included the Coulomb interaction between the electrons. He found 
(a) that the Coulomb interactions ‘screened out’ the bare interaction; 
the effective interaction was identical with that calculated by Bardeen 
(1937), (6) that the longitudinal lattice frequencies were drastically 
renormalized, and (c) that an electron-electron scattering term— 
corresponding to the exchange of virtual phonons—appeared in the 
Hamiltonian. This interaction is usually assumed to be responsible for 
superconductivity. This paper is probably the most important in the 
- whole field and we devote a considerable portion of $3 to a discussion of 
it. ‘The methods developed by Bohn and Pines (1952) for an electron gas 
were extended to take account of the lattice degrees of freedom by Bardeen 
and Pines (1955). This paper contains a detailed account of the method 
and we shall not discuss it in this review. 

It has recently been proposed (Chester and Houghton 1959) that the 
parameter (m/J/)"* should be considered as the basic small parameter in 
electron—phonon coupling problems. This parameter is of course always 
assumed to be small but it is often assumed that the effective 
electron—phonon.coupling constant is also small. This is a much more 
dubious assumption. We shall discuss the consequences of this 
-assumption in §3.3. 

Some recent work by Bogolyubov et al. (1959) connected with super- 
conductivity will be discussed in §5. 


§ 2. QUALITATIVE DISCUSSION 


The conventional physical picture of a metal assumes that the electrons 
and ions are essentially decoupled from one another. For example, in 
calculating the specific heat of a metal it is usually assumed that it can 
be written as the sum of a lattice and electronic part and that any 
contribution from their interaction is negligible. Again, in calculating the 
electrical conductivity it is normally assumed that the electrons interact 
rather weakly with the vibrating lattice and that the scattering by the 
vibrating lattice can be calculated by perturbation theory. We should say 
at once that there is a considerable body of experimental evidence which 


P.M.S. 2D 


360 G. V. Chester on the 


supports this simple picture. On the other hand it is equally certain that 
it is quite wrong to think in terms of a model of unscreened electrons 
interacting with unscreened ions in the lattice. Such a model would lead, 
amongst other things, to an absurdly large electrical resistance due to the 
strong scattering of the electrons by the long-range Coulomb interactions 
with the ions. Moreover the frequency spectrum of the longitudinal lattice 
vibrations would not tend to zero, in the long-wavelength limit, but would 
tend towards the plasma frequency of the ions. This result is of course 
quite contrary to what we know experimentally about the frequency 
spectrum of a metal. What we expect to happen is that the electrons 
manage to screen out both the Coulomb interactions between themselves 
and the Coulomb interactions with the vibrating ions. It is just this 
latter screening problem which is one of the central problems in the theory 
of the interactions of electrons with lattice vibrations. 

On physical grounds we expect that the electrons would try to follow the 
motion of the ions in such a way as to keep the system locally electrically 
neutral. For example, when a longitudinal mode propagates in the lattice 
causing local compressions and rarefactions in the positive ion density we 
expect that the electrons will try to respond to the motion of the ions in 
such a way as to screen out the local charge fluctuations. It is in fact 
very likely that the electrons succeed in doing this. We say ‘likely’ » 
because there is certainly no convincing theoretical proof that this is so. 
On the other hand it would be very hard to see why. the electrical and 
thermal resistances of a metal would not be very much greater than they 
are if this screening did not take place. All the theoretical methods we 
shall describe in the next section will in fact confirm this statement. If we 
accept this conclusion for the moment then we would be led to the following 
tentative ‘picture’ of a metal. As the ions vibrate and the positive 
charges are thereby displaced, the electrons respond ‘easily’ to the forces 
generated and flow in such a way as to keep the metal locally neutral. 
Now it is important to realize that the maximum lattice frequency is of 
the order of 10!3 sec and is quite small compared with a typical atomic 
frequency. This latter frequency is of the order of #,/h 10! sec! (#, is 
of the order of a few electron volts). If the electrons are able to respond 
in times of the order of atomic times then they will effectively be 
following the motion of the lattice instantaneously at all frequencies 
of vibration. In other words the electronic motion will be essentially 
adiabatic; the wave functions of the electrons adjusting instantaneously 
to the motion of the ions. The methods we shall describe in §3 assume 
that this approximation is valid, except perhaps for a small group of 
electrons with energies very close to the Fermi energy. We should like 
to present a simple semi-quantitative argument that confirms this point 
of view. 

What we want to do is to obtain an estimate of the importance of the 
non-adiabatic terms in the equations of motion. We do this by writing 
the complete Schrédinger equation in a way which clearly isolates the 
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non-adiabatie terms from the rest. We are then led to an equation of 
the form 


Ge 

Aaa Sot RE ee 

Xa p2 E,— HX 
where 6x,,, is the lowest order non-adiabatic correction to the wave 
function and the matrix elements; C,,,,, are defined by, eqns. (2.5), (2.6) 
and (2.7). The reader who is not interested in the manipulations that 
lead to this expression for $y,,,, can turn at once to the discussion that 
follows eqn. (2.10). Let us write the exact Hamiltonian of the metal as 


A = H+ Hion, 


where He; is the complete Hamiltonian of the electrons including their 
interactions with the ions; Hjon is the sum of the kinetic and potential 
energies of the interacting ions. The Schrédinger equation for the 
electrons with the ions considered to be in fixed positions is then 


Hex(R)y,,(r|R) = £,,%5,,(r|R), eR ied Pas 


where we have written Hei(R) to indicate that the electronic Hamiltonian 
is parametrically dependent on the positions R of the ions. Likewise 
#,(r|R) indicates that the eigenfunctions also depend on the ion coordinates 
R. The exact Schrodinger equation for the entire system is 


He DAeae Oe eae sete claoe ce, 4 (0:5) 
We can expand ©,, in terms of the complete set db, 
®,,= LY Xmn(R)pn(r]R). = ewes «ne -2) (2-8) 


If we substitute this expansion into eqn. (2.2) and write the equation in 
matrix form we find that the y,,,, must satisfy (Born and Huang 1954) 


(Tion + Vion Et DanlXnn > Con Xmn = 9 mXmn> : (2.4) 


where . io 
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Equation (2.4) is really a matrix operator equation of the form 
H Ym =EmXms Scere ie, tere Scare) 


where X,, stands for the column vector [VpbeNmae 9s Xian «sould ns 
an operator with matrix elements 


A nn = (Lion + Vion t+ En + Ban)Onw + Cnn . oe (2.9) 
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Now if we set C,,,,=0 then eqn. (2.8) is a diagonal matrix equation and 
each ¥,, that is a solution of it has just one non-zero component. 
Consequently the wave function ©®,,=Y,,, and is entirely adiabatic. 
Thus the C,,,,, are the terms in (2.4) that lead to non-adiabatic motion and 
we wish to estimate their importance. If we treat them as small then to 
lowest order we may write 


0 
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where dy,,, is the non-adiabatic distortion of the wave function y,,,., 

and y,,,“ is a solution of the equation with C,,,,=0. 

Clearly 8y,,,, Will be small as long as the matrix element C,,,, is small 
compared with a typical energy denominator. We do not wish to enter 
into the details of the estimation of C,,,,, this is given in Appendix 1, 
we find that its order of magnitude is 6/1, where 0 is the average velocity 
of the ions, and / is the distance through which an ion must be moved to 
produce an appreciable change in the electronic wave functions ¢%,. The 
important point to notice is that the energy denominators #, — EZ, will, 
for most transitions, be of the order of magnitude of the Fermi energy 7. 
This is because the matrix elements C,,,, only connect states that differ 
by single electron excitations and the Pauli exclusion principle forbids 
all transitions except those from an occupied to an unoccupied state. 
The great majority of such states have an energy separation of the order 
of magnitude of the Fermi energy; consequently #,—E£,,~7. This 
argument obviously breaks down for transitions from just below to just 
above the last occupied state. However, the number of such states is. 
very small. For if we are a distance 6 away in energy from the Fermi 
level then the number of states is equal to the density of states at the 
Fermi level multiplied by 6; i.e. of order (1/7)8. Thus the product of 
the number of states and the energy denominator is of order of magnitude 


and is again of order 1/n. Let us suppose for the moment that this special 
group of states might also be treated in the same way as those which lie 
well below the Fermi level. Then 


—Hy In 
Now 6~ 10° cm/sec, 7~5ev and / cannot be much less than 10-8em. 
With these values we find 


ho 
= ¥2x 10-4, 


If this, estimate is at all reasonable then the non-adiabatie correction to 
the wave function is extremely small. The essential points in our argument 


_ Interaction of Electrons with Lattice Vibrations in Metals 363 


are (a) that the Pauli exclusion principle makes transitions with small 
energy denominators very unlikely and (b) that the velocity of the ions 
is quite small, The estimate we have obtained for the ratio C,,,,/E,, —E». 
is very similar to that conventionally used for the validity of the adiabatic 
_ approximation in electronically saturated systems with discrete electronic 
levels. In fact (Peierls 1955) the usual criterion differs from ours in that 
718 replaced by the minimum distance between the discrete levels. Since 
this distance is usually of the order of electron volts the two criteria are 
numerically identical. We may put our conclusions in a more picturesque 
language by saying that the Pauli principle ‘freezes out’ most of the 
non-adiabatic transitions and forces the electrons to respond continuously 
to the motion of the ions. We must now return to the problem of the 
small group of electrons which lie very close to the Fermi level. The 
argument we gave to show that these did not contribute abnormally to 
8Xmn can only be regarded as suggestive. Moreover it is precisely this 
group of electrons that dominates many physical phenomena. For 
example, the specific heat and electronic paramagnetism are due entirely 
to the response of these electrons to thermal and magnetic perturbations. 
It is therefore of the utmost importance that they be treated very carefully 
if we are to have any confidence in our results. It seems likely that if a 
property of a metal depends on the entire distribution of electrons in the 
Fermi sphere then the non-adiabatic terms in the equations of motion 
can be safely neglected. Examples of such properties are the cohesive 
energy and normal mode frequencies... On the other hand properties such 
as the specific heat and spin paramagnetism could be appreciably modified 
by the non-adiabatic terms. As far as we are aware none of the present 
methods of dealing with the electron-phonon interaction problem gives 
an accurate and convincing treatment of this special group of electrons. 
We can estimate how close to the Fermi surface we have to be for the 
adiabatic approximation to break down. If we ask how big C,,,/L,— Ey, 
has to be for Syn ~ Xmn” then we find that 


4 


Pr 
|Z, ~Eyl~ - < 10-4 ergs ~ hw max, 


where wmax is the maximum lattice frequency. We may therefore define 
a ‘danger zone’ of electrons that are within iw max of the Fermi energy and 
any property that depends significantly on these electrons must be 
calculated taking into account the non-adiabatic terms in the equations 
of motion. In particular it is very doubtful if any of the current methods 
gives a sufficiently reliable treatment of them to provide an adequate 
basis for a theory of superconductivity. 

We should now like to give a brief physical description of the methods 
we shall discuss in detail in §3. We shall describe three methods, that 
due to Nakajima, the Born—Oppenheimer method and the Hartree-Fock 
self-consistent field method. The Hartree-Fock self-consistent field 
method is essentially no more than a particular way of handling the 
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series of equations generated by the Born—Oppenheimer method and we 
need say very little about it at this state. 

The method due to Nakajima takes into account the long-range Coulomb 
forces between the electrons and ions. He treats the electrons as moving 
nearly (but not quite!) adiabatically. They are therefore able to screen 
out the long-range Coulomb interactions of the ions with each other. 
The longitudinal mode of vibration of the ions corresponds to the 
propagation of an ordinary sound wave rather than to the propagation 
of a plasma wave. The electrons also screen out the Coulomb interaction 
of the vibrating ions with the conduction electrons; the treatment leads 
to the usual effective screened matrix element that is used in the theory 
of electrical conductivity. The other result derived by Nakajima is that 
the Hamiltonian can be written with a new electron—electron interaction 
term in it. This term arises from the absorption and emission of virtual 
phonons by the electrons. These results are obtained by applying a 
canonical transformation to the Hamiltonian so that.the unscreened 
interaction between the electrons and vibrating ions is removed and 
replaced by a screened interaction. At the same time the frequency of 
the longitudinal mode is ‘renormalized’ so that it becomes well-behaved 
for long-wavelengths. The important technical point is of course that it 
is quite impossible to perform the canonical transformation exactly and 
we have to be content with an approximate expression for the transformed 
Hamiltonian (see eqn. (3.2)). In §3 we shall reach the conclusion that 
the transformed Hamiltonian is very likely to be an accurate approximation 
for most purposes, but that it may be quite unsound for electrons in the 
‘danger zone’. The transformed Hamiltonian contains energy denomi- 
nators that are quite small in this region and even vanish for real (energy . 
conserving) processes in which an electron absorbs or emits a single phonon. 
This latter defect can be remedied by modifying Nakajima’s method to 
take special account of these real transitions (§3.2). This modification 
does not however remove the difficulty for other transitions in the 
‘danger zone’. We have already mentioned that the method does not 
quite treat the electrons as moving adiabatically. The transformed 
Hamiltonian contains terms which represent both adiabatic and non- 
adiabatic motion. The latter terms give rise to a large modification of the 
electronic specific heat. This modification is at present an unobservable 
effect because it has exactly the same form as the usual electronic specific 
heat and is therefore indistinguishable from it experimentally. 

The basic ideas behind the application of the Born—Oppenheimer 
method (§3.3) to metals are firstly that we expect the motion of the 
electrons to be largely adiabatic and secondly we wish to make use 
of the fact that the parameter «= (m/M/)"4 is always less than 1/10. The 
importance of this parameter is that the Born—Oppenheimer method 
generates a series expansion for all the physical properties of a metal in 
powers of x; in fact in most cases in powers of x? rather than «. This 
type of expansion immediately opens up the possibility of avoiding 
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expansions in terms of coupling constants which are hardly ever as small 
as i and never as small as «?. If we did not believe that the motion of the 
ions is essentially adiabatic then there would be no point in applying 
this method ; it does not yield any non-adiabatic corrections to the energy 
levels until terms of order «* and the wave functions are adiabatic to 
order x* (Born and Huang 1954). From what we have said previously 
it should be clear that this method is also likely to be unsound for electrons 
in the ‘danger zone’. The advantage of the method over Nakajima’s 
lies in that it provides a perfectly systematic approach to the problem. 
Moreover the expansion parameter x? is really very small and so we may 
expect fairly rapid numerical convergence of the expansions. Any simple 
estimate of the convergence always confirms this prediction. Since the 
method is essentially an adiabtic approximation it produces very similar 
results to Nakajima’s. Only the important difference lies in the fact that 
Nakajima’s method contains non-adiabatic terms in low-order approxi- 
mations. The large modification to the specific heat therefore does not 
appear in this method—it can only be obtained by performing a partial 
summation on the higher-order terms in the series. Both this method and 
Nakajima’s suffer from the technical defect that to carry out the approxi- 
mations consistently one really has to solve exactly the problem of the 
conduction electrons moving in the static crystal lattice and interacting 
via the long-range Coulomb potential. The Hartree-Fock method 
provides a possible way of approximating to the solution of this basic 
problem. In fact at each stage of the Born—Oppenheimer method we 
can solve the appropriate equation by the Hartree-Fock method. We 
outline this procedure in § 3.4. 

We hope that the discussion in this section has provided a reasonable 
physical picture of the problems involved in the electron-phonon problem. 
In the next section we shall give a more detailed mathematical discussion 
of the methods. 


§ 3. Basic METHODS AND APPROXIMATIONS 
3.1. The Basic Hamiltonian 


We begin this section with a discussion of the Hamiltonian of the system. 
The purpose of this discussion is to show just what approximations are 
usually made in setting up the Hamiltonian, for, the electron-phonon 
problem in metals. 

We start from the Hamiltonian 


Ay=T eat VatT uct Vauct Vint Suen Jame 16s Pal) 


Here we have assumed that the usual non-relativistic approximation is 
valid—all spin-orbit and spin-spin interactions have been neglected. 
We believe that these cause no real difficulties and could easily be included 
in the Hamiltonian. The operators in eqn. (3.1) have the following 
meanings; 7'e) is the kinetic energy operator for the electrons, Tnuc the 
kinetic energy operator for the nuclei, Ve1 and Vnue are the Coulomb 
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interaction energies of the electrons and nuclei, and Vint is the Coulomb 
interaction energy between electrons and nuclei. It is to be emphasized 
that all three interaction energies are at this stage assumed to be unscreened 
Coulomb interactions. The Hamiltonian (3.1) is quite general and contains 
no specific reference as to whether we are dealing with a metal or not. 
To overcome this difficulty a more explicit model is usually introduced. 
The Hamiltonian H, is replaced by an approximate one H, which represents 
a collection of ions, rather than bare nuclei, interacting with a sufficient 
number of electrons so that the whole system is electrically neutral. 
This implies firstly that the number of electronic degrees of freedom in 
To; and V1 is reduced so that only a small fraction of the total number of 
electrons contributes to these terms. Secondly that the interaction 
energies Vion and Vint have now to be modified because at small distances 
of approach the electronic clouds around the nuclei will modify the 
Coulomb interactions. Both Vion and Vint will, however, be Coulomb 
potentials for large distances of separation. With these explanatory 
remarks in mind we can write: 


A,=T e+ Veit Vion + Vion+ Vint. : . (3.2) 


There is little reason to doubt that this approximation to H, is sufficiently 
accurate for almost all purposes. 

It is next assumed that there is an equilibrium configuration for the 
ions and both V,,, and Vint are expanded in a power series in the 
displacements of the ions from these equilibrium configurations. More 
precisely if #(R,...Ry) is any function of the coordinates R,...Ry of 
the ions then we write 


F(R, cee Ry) = F(R,°+ KU,, sees Ry°+ KUy) 
= FO(R,°...Ry°)+FO4 FO... 


where F® is linear in the displacement «xu,...«uy, and F™ is of the 
nth order. We have introduced the conventional parameter « into this 
expansion ; we shall see shortly that our assumptions are only consistent 
if «=(m/M)"/4, where m is the electronic mass and the ionic mass. 
Expanding in this way we have, 


1 » Hi”, i, Se ne eae . (3.3) 
where ‘ 
Hy (0) — = Tei + Veit Vion aS fa in {0 . . . : (3. (3.44) 
A= V ion? ) 4 Vine™, 5 Se oe oes (3. 4b) 
A = Vion + Vint® + Vion he eee (3.4¢) 
Fe Vint™ + Vion™, (pay + whe wagk Ae) 


If «=(m/M)"4 then Tion is of the same degree in k as Vion® and Vint®, 
and when this is so the simplest equation to determine the motions 
of the ions is the familiar harmonic equation; it is therefore consistent 
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with the assumption that the ions have an equilibrium configuration. 
The fact that «=(m/M)/4<10— for all metals suggests that this 
expansion scheme may be expected to ‘converge’ quite quickly. 
However, it must be borne in mind that some of the effects we are 
interested in are extremely small—for example, the difference in energy 
between the normal and superconducting states of a metal is about 
one part in 10°—and we should be careful in assuming that all the properties 
of a metal can be calculated from the first few terms in this expansion of H,. 
In the methods that we shall review in this section the series (3.3) is 
truncated at the third term, we call this approximate Hamiltonian H,: 


HoH pH Hele 2 es RB) 


where for simplicity we have dropped the suffix 2 on the terms on the 
right-hand side of this equation. The Hamiltonian H, is essentially 
identical with that used by Bardeen (1956), Nakajima (1953), Toya (1952), 
Bardeen and Pines (1955) and Chester and Houghten (1959). The only 
difference is that for many purposes the term Vint” in H® can be neglected 
and this has been done by all the authors referred to above. We wish to 
emphasize that H, is a reasonable and consistent approximation to the 
exact Hamiltonian H,. It is not however at all clear how justified one is 
in computing it from quantities which depend on higher powers of « than 
the second. This is simply because we have no guarantee that the results 
of these calculations would not be seriously modified if higher order terms 
in H were included in the calculations. 

It is customary to introduce a number of more technical simplifications 
in H;. By this we mean simplifications that are made in the detailed 
structure of the terms in H,—particularly simplifications connected with 
the band structure of the one-electron states. These are usually necessary 
before any but the most laborious calculations can be made. We first 
remark that it is convenient to re-define the three interaction terms 
Vion, Ver and Vint in H so that each of them is by itself finite. The point 
here is that each of them is Coulomb-like at large distances of separation 
_and consequently each of them contains an infinite electrostatic energy. 
The three infinite terms of course exactly cancel one another but it 
simplifies our considerations if we simply subtract from each interaction 
the infinite self-energy; the sum of three subtraction terms is of course 
zero and H, is unaltered. With this subtraction procedure each interaction 
term is itself finite. The simplifications that we shall make to H, are 
most easily formulated in terms of a representation for the operators. 
We therefore introduce a second quantization scheme for the electron 
degrees of freedom, using the one-electron Bloch states of the electron 
moving in the static crystal potential Vint} as the basic states. 


+ An alternative set of states would be the solution of the Hartree-Fock 
equation for an electron moving in the lattice. This set of states however 
probably has the defect that it partially includes the effect of the long-range 
Coulomb interactions and unless careful account if taken of the remaining part 
of the Coulomb interactions it could lead to absurd results. 
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The lattice degrees of freedom are most easily handled in terms of the 
standard normal coordinates 58 in terms of which we have 


xu, =(NM)> mp2 fi (f, s)q¢, ,exp (0R,;°. f), ie asan 


where f is the usual wave vector, s the polarization index and V,(f,s) 
the set of polarization vectors. 
The Bloch wave functions satisfy the equation 


(ter +O int), = (A), ety ity Steet Bele 


where te) is the kinetic energy operator for a single electron and Vint the 
static crystal potential for a single electron. 
We now write out H®, H® and H® explicitly : 


H= > a*a,e(k) + > Vee tet tes One, oe a cae ae 
k kk’f 

AD = > Vy? etess* es ee Oe eee eg SING 
kf 


H® = 4 » (De,2+Q¢ 2G 2?) + 2 Ve ester. « -  (3.8¢) 
»8 kft’ 
The sums over f and k vectors are to be taken in the extended zone scheme 
to allow for Umklapp processes (Peierls 1955) and a summation over 
spin indices is implied in the k summation. The three matrix elements 
Vines Viv and V2), are given by the equations: 


Vier = Shanes rs )bi(ri)p(r;) dr, dr, ves Ha gay 
Vie=— (NM) ae : Pa , 8). V,vint(r — R,°) exp (— if. Ry) dade 
(3 9.b) 


Vite = (NM) [cea LE Vy f,s).V,)(V;(f, 8’). V,)vint(r — R,°) 
J < 
x exp (—i(f+f’). Ry) bade : (3.9¢) 


The frequencies Q,;, that appear in H® are the normal mode frequencies 
of the Hamiltonian, 


Tion + Vion®. 


We have omitted the term Vion from H® because this term vanishes 
identically for any periodic lattice (Peierls 1955). 

Our first simplification is to assume a monovalent metal with one atom 
per unit cell. It is a simple matter to extend the formalism to more 
complicated situations. This means that there are only three values for 
the polarization index, s. 
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We assume further that the term 
p (Pes +24 674¢5") 

in H® can be split into a sum over two transverse modes and one 
longitudinal. This assumption is of course strictly valid only for very 
small wave numbers for it is only for these wave numbers that the three 
modes of vibration have either purely longitudinal or a purely transverse 
character. However, as many of our most important considerations will 
refer to these small wave numbers this is a not too drastic approximation. 

The frequencies Q,, of the transverse modes tend to zero as [f| tends to 
zero. However, the frequency of the longitudinal mode tends to a 


limiting value 
4ane*\ 12 
(dour ( = ) = OF. 


where Q,, is the plasma frequency of the ions in the lattice and n is the 
density of ions. The physical reason for this is that the frequency Q, 
has been calculated without taking into account the screening of the ions 
due to the conduction electrons. Consequently the long-range Coulomb 
forces between the ions dominate the small wave number limit of 
(26): ong: 

A rather more important approximation is that we shall assume that 
in Vint® coupling occurs only between the longitudinal mode of vibration 
and the electrons. This approximation is discussed in detail in Wilson 
(1953) and Ziman (1960). It is only correct for highly isotropic metals. 
Nevertheless we shall introduce it to facilitate our calculations—we believe 
that the qualitative results we shall find will remain valid when coupling 
to transverse modes is introduced. 

Our final simplification is to assume that all the oe elements of 
Ver, VY and VQ, depend only on the difference between the wave 
vectors of the initial and final states. 

We can now use these simplifications to rewrite eqn. (3.8) as 


H= > ad, €(k) iy V,< > Ug f *Qiest™ AN, . . . . ° (3.12) 
k k,k 
H®= DS V—& 2 Oy Oy 60 f; . ° . . . ° ° ° (3.13) 
f k 
H®= § > (Deer Oe a8 5°) Vets Dae ocaageVe 
ff" 
cant NE 
a EDA (p+ Q-a¢”). sweet (3. 14) 
It follows from the long-range nature of the Coulomb interactions that, 
V ordre? Vt fae eee eee (LD) 
[V | ™ Q,|f|-1 Se Ree ee Mee sacaine (3.16) 
and j 
Og (47€2n/M)42 eda aera eo. kd.) 


for small |f|, Bardeen and Pines (1955). 
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3.2. The Nakajima Transformation 


In this section we shall describe the canonical transformation introduced 
by Nakajima (1953) and Bardeen (1956). This transformation is a 
generalization of the one used by Fréhlich (1952). The basic idea is to 
transform the Hamiltonian so that the electron and ion degrees of freedom 
are decoupled to high order of accuracy. The ‘largest’ coupling term 
in H, is the term Vint, this is first order in the ionic displacements. More- 
over the matrix element in Vint is unscreened (see eqn. (3.16)) and would 
lead to entirely wrong temperature dependence for the electrical 
conductivity. We have already pointed out that the frequency of the 
longitudinal mode of the ions does not vanish for zero wave numbers ; 
this would lead to an exponential term in the lattice specific heat at low 
temperatures. The aim of Nakajima’s transformation is to (a) remove 
the first-order unscreened coupling term Vint” and (b) to replace the 
‘bare’ ion frequency 2, by a screened frequency w,. It has become 
customary to refer to the Q,’s as the unrenormalized frequencies and the 
w,’8 as the renormalized frequencies. It is of course essential that the 
w,~|f|, as |f|>0, in order that the lattice specific heat have the correct 
temperature dependence at low temperatures. 

We expect on physical grounds (see §2) that the response of the 
conduction electrons to the ionic motion screens the moving ions and 
hence renormalizes the frequency of the longitudinal mode. This implies 
that any mathematical transformation of H, which achieves this renormali- 
zation must handle the Coulomb interaction between the electrons to a 
reasonable approximation. The Nakajima transformation does in fact 
do this and it is for this reason that we regard it as an important method 
in the electron-phonon interaction problem. Before describing the 
_ transformation we should like to comment on the fact that the method 
‘transforms away’ the unscreened interaction Vint. We have already 
pointed out that in some respects this interaction term has quite the 
wrong structure and would lead to absurd results for the electrical 
conductivity. However, it is the only term in H,; that can lead to real 
(energy conserving) one-phonon absorption and emission processes, and 
it is just these processes that we believe give rise to the electrical resistance 
of an ideally pure metal. Consequently if we transform this term away 
we are apparently left with no means of accounting for the resistance in 
terms of real one-phonon processes. An immediate objection to this 
argument is that if we remove Vint” by a canonical transformation then 
a very similar interaction is likely to reappear in the higher order terms 
in the transformed Hamiltonian. We shall see however that the new 
terms which do appear do not have finite matrix elements on the energy 
shell for one-phonon processes ; consequently they cannot lead to a finite 
resistance. This difficulty with Nakajima’s method does not seem to 
have been discussed in detail before. We show that a simple way round 
the difficulty is to perform two transformations, one of which is confined 
entirely to the region of the energy shell for one-phonon processes. ‘The 
transformation can in fact be chosen so that all the higher order terms are 
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finite on the energy shell and the effective electron—phonon interaction is 
properly screened. 

Nakajima starts from the Hamiltonian H 3 With the term Vjn¢? omitted. 
This term does not contribute to the normal mode frequencies and merely 


leads to a two-phonon scattering term. Now H, can be written in the 
form 


A,=H®+ V9 + Hp + (App —Hpn™)+ Htrans. . (3.18) 
Here 
Ay™= > (pP+orq oe Ma Pee (eo) 
and 
Hp = ¥(pP+Qege). 6... « (3.20) 


f 


Both these terms in the Hamiltonian refer only to the longitudinal mode; 
Htrans is the Hamiltonian for the transverse modes. We now transform 
H, by means of a canonical transformation S* 


H,>H,' =exp (—iS*)H,exp(+i8). ress 4 
We have chosen our units so that #=1. To second order in Sy, H;' is given 
by 
ei ni [[9*,[S",H]].  . . (8.22) 
We choose S, so that 
—i[S¥, H+ Hyp] + Vine =0 Moy es) 


and the renormalized frequencies w, so that, 


5 Mao 5) lat pn — Bn) == 05-3 a (3.24) 


where the suffix d.p. means that the part of the operator proportional to 
9;9;* is to be taken. When these two equations are satisfied 


: : ”) 
fale = H°+ Hpn” + Htrans + 9 [Vint™, Sees. (3.25) 


where the suffix n.d.p. means that the non-diagonal part is to be taken. 
This equation is correct to second order in S* in the following sense. 
If we count H° and Hpn” as being of zero order, and V jn¢"? as of first order, 
then (3.23) is consistent. Then Hp,” —Hy»® must be taken to be of 
second order if eqn. (3.24) is to be consistent. 
If Sy is defined by the equation 
SP =% 2, O* (Pk, f)qg, —iG'(k, f) p¢), ns euee sag) 
k, 

and F and G satisfy the integral equations 


Y (2V 5° — Ve_,)(m(k’ —f) — nk’) ) F(R, f) + (Wie — Wi) F(k, ff) 


+ (a ee) F(k, f) +o 2O(k, f) + V=0, ele ea Bee) 
and 


YT OVE=VE gill’ A) mk, + (Wie Wad 
Sb eg dk, f) + FU f)=0, (8:98) 
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“a 


where W,, is the usual Hartree-Fock exchange energy, 


Wy='d Vouk =f) oe 2 oes 
- 


then eqn. (3.23) is satisfied as a matrix equation between any pair of 
eigenstates of H. In eqns. (3.27), (3.28) and (3.29) n(k) is the number 
operator aa, and should be treated as a diagonal operator in evaluating 
the commutator in eqns. (3.23) or (3.24). We now examine the transformed 
Hamiltonian H,’. The frequencies w, are given implicitly by the 
eqn. (3.24) which can be written as 


wP—-OQP=> E(k, f)V[n(k) — n(k—f) J. -/ son (S280) 
k 
No exact solutions of the integral eqns. (3.27) and (3.28) are known. The 


approximation that Nakajima made is to ignore the exchange energies W,, 
and the other exchange terms ; 


> Vie—u(m(k’ — f) —n(k’))F(k' f) | 
and Peeters ee es) al 
2 Ver—ulk—f)— ake NGA). 


When this is done the equations can be solved to give 


Gk, f) = — a<S - (3.32) 
and sea 
(<x — €_f)U¢ 
F(k,f)= 
USE (cae yaaa a (3.33) 
Here 
ve=V [1+ VES(f)]A w scan Mays “Gi Uae arene 
with 
S(f) = Shree) (3.35) 


ie Ske ae Cee 


The function 2, is usually referred to as the screened matrix element and 
the factor multiplying V, on the right-hand side of eqn. (3.34) as the 
screening factor. S(f) is calculated explicitly in Appendix (3) where it 
is shown to be a rather slowly varying function of |f]. When |f| tends to 
zero S(f) tends to a constant, and v, is for small |f] given by the 
equation : 
Alf * 
= Be (3.36) 
_where A and B are constants. This is exactly the form we would expect 
for the matrix element if the interaction between the electrons and ions 
were screened Coulomb interactions rather than bare Coulomb interactions. 
In § 3.3 we shall see that this approximation to eqn. (3.27) is closely related 
to the Hartree self-consistent field treatment of the problem. 
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With the same approximation we find 


we = OP —0,V_OS(f) e e ° e e (3.37) 
and for small |f| 
QQ? 
Beige My Uae 38 
Wr TG) (3 3 ) 
This equation implies that w, tends to c,|f| for small |f|. In a free electron 
model this leads to the simple formula, 


Le fam \82 
es= 5 (Fz) ae Ms tor weet et (S30) 


where v, is the velocity of the electrons at the Fermi surface and m is the 
mass of the electrons: This expression for c, is certainly of the correct 
order of magnitude. We conclude that this simple approximation leads 
to results which are certainly qualitatively correct; whether one could 
trust this approximation quantitatively is a much more difficult question. 
In §4.1 and §4.2 we shall discuss the quantitative comparison of the 
above results with experimental data. 

The main difficulty in proceeding further with eqns. (3.27) and (3.28) 
is first that W,, is a fairly complicated function of |k| and secondly that 
the terms (3.31) make them non-trivial integral equations. It is necessary 
to take both these terms into account, to obtain a realistic improvement 
on the ‘ Hartree-like’ approximation. There have been two attempts 
to improve on this approximation (Hone 1960, Bailyn 1960). Hone 
starts from the full set of eqns. (3.27) for F and G. He then makes the 
approximation of neglecting w, compared with («,—«_,). This is for 
many purposes an accurate approximation because. «,—«_, will be of 
the order of the Fermi energy for most transitions whereas iw, can never 
be greater than a few per cent of this. It is however unlikely that this 
approximation can be valid always; for example, it is almost certainly 
a poor approximation if one wishes to calculate an effective electron— 
electron interaction to introduce into a theory of superconductivity. 
We shall discuss this point in more detail in§ 6.1. When this approximation 
~ is carried out eqns. (3.27) reduce to a single integral equation. If we define 
an effective matrix element v(k, f) by the equation . 


v(k, f) 
= 3.40 
Fk, =p (3.40) 
where 
i= Clea Wi: 5 . . e 5 (3.41) 
Then ; 
v(k, f) = VP + > A(k, k’, f)v(k’, f), = ie: «+, (3.42) 
“ip 
where 


Cea Tc PS ; k’ 
Aone ee teen (BRAG) 
kf — SK 
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These equations are identical with those derived by Bailyn. His method 
is to extend the Hartree approximation (Bardeen 1937) for the problem 
to include the exchange terms. In eqn. (3.42) the sum over the wave 
vector k’ is weighted with the exact occupation numbers n(k’) for whatever 
exact state of H® we are interested in. In Bailyn’s treatment the 
corresponding sums are to be taken over the unperturbed Fermi 
distribution %(k’). It is of course not known how different the exact 
occupation numbers 7(k) are from the unperturbed ones. The reader is 
referred to a recent paper by Daniel and Vosko (1960) for a comparison 
between the unperturbed n(k) and their values for a dense electron gas. 
It is possible to express the renormalized frequencies directly in terms of 
v(k, f): 


We = OQ? - lV MPRA), Fy 4 ° . < (3.44) 
where 
= k — f) —n(k) 
R(f)=|V|-2 idl Saat rete! oes) 
(f) | f | 2, Ee ( ) 


This leads to the following expression for the velocity of sound, 


[ 1— VPR(f) A. 


2071 
CS eG 


I fl>0 


Both Hone and Bailyn have found approximate solutions of eqn. (3.42). 
It is not clear how accurate these solutions are in fact; the reader is 
referred to these papers for the details of the methods used. As we have 
remarked above, the equation does in principle contain the unknown 
n(k)’s. Hone replaces these by their free electron values and so his 
equations are finally identical with Bailyn’s. One could question whether 
any major effort should be expended on eqn. (3.42) as long as the n(k)’s 
have to be replaced by what appear to be rather crude approximations. 

We next turn to the diagonal part of [Vint, Sy] which is the remaining 
term in the transformed Hamiltonian. We have 


a 
m3 [ Vents aa py = 2 Vg Pee ge yy" Cy ¢G(k, f) sh > Meade + > NeedsPe- 
iF ff’ 


The second and third terms on the right-hand side of this equation are 
two-phonon terms and are likely to be small; no discussion has been given 
of them in the literature. The first term can be written in the form 


vs|?n(k)(1—n(k —f : 
et Mk eS Vee eax sk, f). (3.45) 


ik (Exe 5)? — we? k+k,f 


The first diagonal term in this expression is exactly what we would have 
found if we had done second-order perturbation theory with a screened 
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matrix element %,. It is fact identical with the expression used by Frohlich 
(1950, 1951) in an attempt to construct a theory of superconductivity. 
We shall discuss the effect that this term has on the thermodynamic 
properties of non-superconducting metals in §4.1. The off-diagonal term 
has the form of an interaction between the electrons. The interaction 
arises from the absorption and emission of virtual phonons; is thought 
(Bardeen 1956) to be responsible for the transition of a metal from the 
normal to the superconducting state. We remark at this point that it 
depends, through the frequencies w,, on the mass of the ions in the lattice 
and can therefore in principle at least lead to an isotope effect (Swihart 
1957). ‘This term will be discussed further in § 5. 

It is clear from the structure of the functions / and G in the Hartree 
approximation that all the higher order terms in the transformed 
Hamiltonian will contain energy denominators of the form 


Ex — €x_# + hows. 


It is possible to transform the integral equation for v(k, f) to show that 
it is a well-behaved function when «,—q&_,= +hw, and consequently 
even if the more complicated equations for F and G were solved exactly 
we should still find that all the higher order terms would have the same 
energy denominators as the Hartree approximation. The appearance of 
these denominators raises two questions: (a) how can we hope to find terms 
in the transformed Hamiltonian which correspond to real one-phonon 
processes in which «,—«_,= +a, and (b) how do these denominators 
effect the convergence of the series which defines the transformed 
Hamiltonian The real one-phonon unscreened interaction Vin¢® has 
been entirely removed from the Hamiltonian, and since both functions 7 
and G become infinite for these transitions it is impossible that any 
analogous term should appear in the higher order terms in H’. A possible 
way round this difficulty was first suggested by Frohlich (1952) and later 
by Bardeen and Pines (1955). As it has not yet been fully developed in 
the literature we shall outline the procedure here. The basic idea is that 
the canonical transformation we have discussed should be carried out for 
only that part of Vint that lies off the energy shells given by the 
equations : 


Ex — eat = we. 


It is convenient to split Vint? into two parts. We do this by defining a 
function 0(A, e&.— &_¢+hw,) such that 


O(A, &e— &e_gtho,)=1 oN ere e yg b Th <= Al 


3.46 
=( otherwise. (9-26) 


The parameter A defines the width of the shell; we leave it arbitrary for 
the moment. 


P.M.S, 2E 
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Then we can write 
Vint? = Us + Ugg, . . . ° . (3.47) 


where 
Ug = > Vay eyed HA, ie — e—# + rs) . -. (3.48) 
k,f 
and 
Oyg= »S V Pay *ay_194(1 — (A, €y — &e_¢ + hews)). . (3.49) 
k,f 


We now introduce a modified Nakajima transformation S,°*) which 


satisfies the 
i [Syo*, H® + Hpn™)] + Uo, =0 os ee OD 


and introduce a new transformation S,* on the energy shell which satisfies 
i[Sy5, fO1+ 0.82055” 9) ees 


If we set 
Sy =i Dd aykay_-(F(K, f)q¢— 1G5(k, f)p_¢) 
kf 


then eqn. (3.50) is satisfied if F°% and G°* satisfy the equations 


Y (2V 0 — VE_,)(n(k’ — f) — n(k’)) FOS(k’, f) + (Wir — Wi) F%(k, f) 
= 
+ (ee — €_¢)FOS(k, f) + w2@5(k, f) + V1 — 0) =0 aaa 
and 
¥ (2V C= Vy Sp )(ra(e’ — f) = 0(k'))G95(k', f) + (Wye — We )G°S(k; 
= 
+ (€x — €x_¢)G9S(k, f ) + FS (k, f) =0. cm ye ee 


On the other hand we can satisfy eqn. (3.51) with 


Sy => >, FO(k, f)q Oy *Ay_¢ ' - + 5 (3.54) 
k,f 


if F) satisfies the equation 
{len— Wi) — (ee_1— Wig) JPA (ke, f) + 2 2 V¢°(n(k’ — f) — n(k’)) PS(k’, f) 


+ > VE_,(n(k’) —n(k’—f))FS(k’, f)=V~MPO. 2. (3.55) 
- 
We can now define the renormalized frequencies &, by the equation 
v (Jos 8 0S a Us y § (Db) (r) 2 BR 
5| mo) By’ dap, + 5 | 1") Swe lap. + (Apn® — A pn) = 0 (3.56) 
or more explicitly, 


HP-OP= ¥ F(k, AV _ 1-0) + D Pk, f)V_ «8. (8.57) 


k k 
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From this equation we can see that &, differs from w, by terms of the 
order 
n(k)@ A 
ee Pree yk S580 5. (3.88) 
iof k— &e_f— hw, 
Consequently as long as we choose A<» the new frequencies will differ 
negligibly from those defined by eqn. (3.37). 
The remaining terms in the transformed Hamiltonian are the 
off-diagonal parts of 
[Sy, Ugg, [Sy®, U gl, 


and a one-phonon interaction term 


[Sy pO l=t>- a", F5(k, fp et. 5... (9:59) 

kf 
The off-diagonal parts of 
: [Sy, Us| 


are of exactly the same form as those of 
[Sy®, Ugg] 


_ except that they are only non-vanishing on the energy shell defined by 
eqn. (3.46) and have energy denominators of the form «,— «,_, instead of 
€e— &_+ thw,. This implies that the residual interaction term between 
the electrons is quite large. It is however no larger than that given by 
eqn. (3.45) when «,—e&_-<hw,. The most significant point about the 
residual electron—electron interaction terms is that the energy denominators 
now never vanish. In fact we can let the width A of the shell tend to zero 
in both terms and the final interaction will then be identical with that 
given by eqn. (3.45) except that the singularity due to the vanishing 
energy denominator will now be defined to be a principal parts integral. 
We now turn to the one-phonon scattering term, eqn. (3.59). It is 
usually assumed (Wilson 1953) that it is justified to introduce a one-phonon 
interaction term, with a screened matrix element, into the Hamiltonian 
of a metal. If we call the screened matrix element é, then the form 
assumed for the interaction is 
Dd dye. Pere aetesry. oP a( 3.50) 
fk ; 
The interaction (3.59) differs only slightly from this. If we solve eqn. (3.55) 
in the Hartree approximation and notice (a) that on the exact energy shell 
€x— &_¢= hw, and (b) that the matrix elements of ip*,/@, are the same 
as those of g; then the interaction becomes 


Y 1 Ay, Aaah sco ea(S.61) 
k,f 


which is identical with (3.61) on the energy shell if we identify v, with %,. 

Since V, is (see eqn. (3.36)) a screened interaction our one -phonon 

interaction term has just the correct structure, There is however no term 
2E2 
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in the transformed Hamiltonian corresponding to those parts of (3.61) . 
which lie off the energy shell and consequently there is no justification 
for the introduction of this interaction term except for energy conserving 
processes. 

We givé now a brief discussion of the convergence of the series that defines 
the tr ansformed Hamiltonian. The convergence of this series is of course 
the crux of the whole method—if the low-order terms in the series give 
an adequate description of all the physical phenomena we are interested 
in and high-order terms merely lead to small corrections then the method 
can certainly be said to be satisfactory. Now we have seen that the terms 
up to second order (a) define a reasonable set of renormalized frequencies, 
(b) give rise to a residual electron-electron interaction that might be used 
to provide a basis for a theory of superconductivity and (c) can lead to a 
one-phonon scattering term with a screened matrix element. There is 
thus a distinct ogabiliby that these terms could provide an adequate 
basis for the phenomena we are interested in. The convergence of the 
series has never been investigated rigorously and this would be an 
extremely difficult problem. We can however give a simple qualitative 
discussion; a similar discussion has been given by Bardeen and Pines 
(1955). The essential question we have to decide is the order of magnitude 
of the transformation function S. We can write (compare eqn. (3.26)), 


S= ¥ S(k), 
k 


where S(k) is the transformation function for an electron with wave 
vector k. Now the order of magnitude of S is given by 


elas 


————— | n,(1—m_+). 


This estimate follows directly from eqns. (3.26) and (3.33), and we have 
neglected the small term @, in the denominator. We have already seen 
that we can define the singularity in sums like this by a principal parts 
procedure. If this description is used then we find |S|~ (h@max/n)!. 
This crude estimate is therefore quite satisfactory from the numerical 
point of view. The whole point however is that in proceeding in this 
very naive way we have may to treated quite incorrectly those terms S(k) 
in the sum for S which have small denominators. These terms come of 
course from electrons in the danger zone. If we have treated those terms 
incorrectly how reliable is our estimate of the magnitude of S? This 
question cannot of course be answered rigorously until we are able to 
treat the electrons in the danger zone correctly. This we cannot do at 
the moment; what we can do is to show that if these special terms, when 
treated correctly, do not contribute abnormally to S, then our estimate 
is at least numerically correct. We show this by simply calculating S 
with all the terms, whose energy denominators le, —€,_,| are smaller than 
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a certain amount AH, removed from it. If we call this new function S 2 
then it is easily shown (see eqn. (3.58)) that, 


[S4/S]~1+0(AB/y). 


If A<» the magnitude of S4 is identical with that for S. This conclusion 
is in agreement with what we found earlier for the frequencies; eqns. (3.57), 
(3.58). It is also in agreement with our qualitative discussion in $2. On 
the other hand it should be clear that any prediction we make about the 
detailed behaviour of this particular group of electrons must be treated 
with extreme caution. It is in fact believed (Bardeen ef al. 1957) that it 
is just these electrons that are involved in the transition from the normal 
to superconducting states of a metal. If this conjecture is correct then it 
would be doubtful whether one could base a theory of superconductivity 
on the lowest order electron-electron interaction term (3.45). We shall 
discuss this point in more detail in § 5. 


3.3. The Born—Oppenheimer Method 

In this section we shall discuss the application of the Born—Oppenheimer 
method (Born and Huang 1954) to the electron-phonon coupling problem 
in metals. We have already discussed in §2.1 why we believe that this 
method is likely to be applicable in spite of the breakdown of the standard 
criterion for its validity. The Born—Oppenheimer method provides a 
scheme for calculating the wave functions and eigenvalues of a system of 
nuclei or ions and electrons. It is essentially a type of perturbation 
expansion in which the expansion parameter is (m/M)"*; the eigenvalues 
generally contain only even powers of this parameter while the wave 
functions contain both even and odd powers (Born and Huang 1954). 
Since this parameter is never greater than 10-' we may expect that these 
expansions will converge at least numerically. As far as we are able to 
see this is the only small parameter that is available in the electron-phonon 
problem and it is for this reason that we believe it is worth while to explore 
the application of the method to metals. 

Again however all our usual reservations must be made about its universal 
applicability. In particular it is most unlikely to be applicable—at least 
in its naive form—to the problem of superconductivity. If we take « 
as our basic parameter then we shall certainly avoid any reference to the 
electron-phonon coupling constant which can of course be quite large. 
We shall see in fact that each term in the Hamiltonian can be easily 
characterized by a definite power of «, consequently there is no difficulty 
at all in generating the formal expansions we require. We should perhaps 
emphasize that the method is entirely adiabatic until rather high powers 
of « are reached in the expansion. To be precise the wave functions are 
adiabatic up to and including the terms of order x*; the eigenvalues up 
to order x4, This implies that the method is quite unsuited to dealing with 
problems in which non-adiabatic motion is especially important. 
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The standard method of generalizing the Born—-Oppenheimer expansion 
is rather inconvenient for our purposes because we want to compare it 
with Nakajima’s method. We therefore prefer to generate the low-order 
terms in the expansion by means of a canonical transformation of the 
Hamiltonian. The complete expansion can be generated by applying a 
sequence of canonical transformations to the Hamiltonian. We shall 
begin by working to order «?; we can therefore start from the Hamiltonian 
H,. If we transform this Hamiltonian with a canonical transformation Sz 


defined by the equation 
i[Hei, S3]+ Vint =9, vce baw ehaoen 


then to second order in « we find that the transformed Hamiltonian H3,, 
is given by 


1 - a Y . . 
Hy,y= Het Apn + a [V int”, Sp]+ Atrans- . a (3.63) 


Here we have again dropped the term Vint. Equation (3.63) follows 
directly from (3.62) when we realize that S, must be of order « if eqn. (3.62) 
is to be consistent. 

It is a straightforward matter to show (Chester and Houghton 1959) 


that to order x?, H,, gives the same eigenfunctions and eigenvalues as Hy. 
If 


Sx =4 3 Oy Ay _¢ V ard (k, f). . . . (3.64) 
kf 
Then (3.62) is satisfied if 
(ee = Wa) = (ee — Wir) }P(k, f) +2 SV 9O(n(k’ —f) — n(k’))6(k’, f) 
= 


+ ¥ Vie P(n(k’) — n(k’ —f))d(k’, f) = 1, 
- 


(3.65) — 
and then 


5 [Vin®, Sp]=2 2 > VeOV_ age (k, f)[n(k) — n(k — f)] 


+4[Vint®, SeJnon, diag.” Seno eI EN 
The first term on the right-hand side of this last equation can be combined 
with Hpn to give a set of renormalized frequencies w ,;: 
wp =F + Y VOV_Md(k, £)[7(k) — 2(k —f)]. (67) 
k 


We can rewrite H,,, as 


Ay,=Heat+ Ap” + Att bY Vie V ededeb(k, f)(y dy ¢_¢ — Uy +7" My_¢) 
eet 
(3.68) 
If we compare eqns. (3.64), (8.65) and (3.67) with eqns. (3.26), (3.27) and 
(3.30), we see that they are very similar indeed. In fact they become 
identical if we neglect the terms proportional to w,2 in eqn. (3.27): this 
is the approach used by Hone in his discussion of the effective electron 
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phonon matrix element. As we have pointed out, since w,<lex— €x_| 
for most transitions, this approximation is likely to be quite accurate. 
We can easily verify this by solving the eqn. (3.65) in the simple Hartree 
approximation. We find that 


il 
CALS INCE A) Ae aoe an ee ee (3.69) 
eave 3 [a —#)—ni(k I 
k’ EK — Eye _f 
and hence 
wp? =OQ2+ Voy ay RO: Whee ERTL) 
k Ek — kf 
The sum 
i(k) — 2(k —f) 


PHAR DS 
k Ek — €k—f 

is calculated in Appendix 3, where it is shown that it differs by terms of 
order fi@,/n from the sum in which the energy denominators are 
Ex — €&k_-4—hwp. Consequently when we compare the equation for w, 
with eqn. (3.37) for a we see that they differ by terms and of order 
hwp,/yn; that is to say by at most a few per cent. The screened matrix 
element v, is given now by the equation 


n= VO Ve —— | ‘ eietiner rt 
k €k — €k_f 

and differs by a few per cent from that defined by the Nakajima 

_ transformation. 

We see from eqn. (3.68) that to order x?, H,, contains no term that 
could give rise to real one-phonon processes and no terms that represent 
electron—electron interactions via the exchange of virtual phonons. 
However, both these terms arise in higher-order terms in H,. If the 
Hamiltonian is calculated to third order we find a term 


1 > Oy dy V (MV A}(K, f)p¢*. - F , (3.72) 
k, f 


This term gives rise to energy-conserving one-phonon processes and on the 
energy shell is identical with the usual phenomenological interaction 
(see eqns. (3.60) and (3.69)). When «,—«,_,= thw, the matrix elements 
of this term are of order x; it is easily seen that all other terms in the 
Hamiltonian that give rise to real one-phonon processes have matrix 
elements of higher order in «. It is therefore consistent to use the 
interaction (3.72) to calculate the one-phonon scattering. The Born— 
Oppenheimer method therefore has the advantage that it automatically 
gives a one-phonon scattering term of the correct structure. This result 
was first noticed by Ziman (1955) and Goodman (1958). 

To order «? no term of the type given by eqn. (3.72) arises but at order x4 
we do find a term of exactly the type; namely 


2 
n(k)(L—n(k — f))ve + > (off-diag. terms), . . (3.78) 
ict (ee x4)” kif 
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where we have separated the diagonal part of this term from the scattering 

part. The only difference between these terms and those derived by 

Nakajima is that the energy denominators contain no phonon energies. 

We see at once that when e,—«,_, is small these terms become very 

large—larger than those in eqn. (3.45).. In particular if we believe that 

interactions, in which both «, and «,_, are very close to the Fermi energy, 

are important, then these terms will become exceedingly large. This 
difficulty is of course the counterpart of the difficulty we found in 

Nakajima’s method where we saw that the method broke down when 

€x— €k_-¢ —w, became very small. 

The convergence of the method is again very difficult to assess; the 
main difficulty lies in the fact that the amplitude of the transformation S; 
becomes very large when «&,—«_, is small. By the same argument that 
we used in § 3.2 if we cut out an energy shell of width A then we can expect 
reasonable convergence if A <7. 


3.4. The Hartree-Fock Method 

In this section we shall show how the Hartree-Fock method can be 
applied to the problem we are interested in. This method was first used 
by Bardeen (1937) to calculate the effective electron-phonon matrix 
element. Toya (1952) extended it to obtain the renormalized frequencies. 
The method—in so far as it has been used to date—really assumes that 
the adiabatic approximation is valid. We can therefore easily formulate 
it in terms of the Born—Oppenheimer expansion. The Hartree-Fock 
method is very easily set up in configuration space and it is therefore 
convenient to derive the Born—Oppenheimer expansion by the conventional 
approach (Born and Huang 1954) rather than use a second quantization 
scheme of § 3.3. 

The expansion can be generated very simply using eqns. (2.1)-(2.11). 
We proceed as follows. First we solve eqn. (2.1) by a perturbation 
expansion in powers of x. Thus 


hy(F[R) = > 1p, ey Meee ay et 
E,= > «*H,\. it gad, Meat 


s 


Next we use this expansion to obtain the expansion of the matrix operators 
wa ; ’ : 

Cnn’ mn eqn, (2.4) and we also expand Vion in powers of «x. The final step 

is then to solve eqn. (2.4) by a perturbation expansion in «; thereby 

obtaining an expansion for the exact eigenvalue &: 


En = > KE (ne EE age aie aT 
& 


and wave function 


@, = by K 5, 8) F . , : . (3.77) 
We find that : 


C 


P’ (0) — FF (0). C iy est . Oo 
Om = KS uP 7) e — AO Pe ee Ee aS (3.78) 
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mm Z) Vr cs . * . . 
The second-order term &,,° is the eigenvalue of the harmonic approxi- 
mation equation : 


(Lion + Vion® + £,?) ¥nn™ = Gey ee on) 


As we have already mentioned, the whole difficulty in this procedure stems 
from the fact that we have in fact to solve eqn. (2.1) for the motion of the 
electrons in the presence of the fixed ions and interacting with each other 
via their mutual Coulomb interactions. The simplest approximation to 
this problem which retains some (at least) of the essential features is 
provided by the Hartree-Fock approximation. .The only difference 
between eqn. (2.1) and the familiar many-electron equation is that it 
contains an extra interaction energy between each electron and the fixed 
ions. This extra interaction energy is a sum of single electron terms and 
the Hartree-Fock equations for this problem merely have an extra potential 
energy term in them arising from this interaction. If we set up the 
variational equation using a Slater determinant of single particle wave 
functions ¢,(r|R) which are allowed to depend parametrically on the 
coordinates & of the ions we find that they must satisfy the equation 


KH by t+ Ving(t[R) P= Pre Wie ncke er) 
where 
KH y= — v4, > eso ped me ial Pao od) 
noes 


and is the usual Hartree-Fock Hamiltonian; Vint(r|R) is the interaction 
energy between one electron and the entire array of ions in an arbitrary 
configuration R,, R,...Ry. We have therefore reduced the many-body 
eigenvalue (eqn. (2.1)) to a simple Hartree-Fock equation for the single 
particle function ¢,(r|R). We now solve eqn. (3.80) by a series expansion 
in x; we write 


P= Px + 5 Ss Di Pye ie cg Sibel 
n>1 k’ 
and 
€, = 42+ > KE, 
; n>1 
Here ¢,° satisfies the equation 
(A+ Vint) bo = Pb? Pe ay he tee (0d) 


and is of course a plane wave solution with eigenvalue 
=H. = &&— W,. . ; : ‘ . (3.84) 


Here «, is the usual Bloch energy of an electron moving in the static 
potential V(r|R,), and W, is the exchange energy. 
To calculate the frequencies of the lattice in this scheme we require 


LP= > nlkje ne, ay Stes aia Ais SEE (ScD) 
k 
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The n(k)’s are the usual unperturbed occupation numbers 


n(k)=1 |k| <ks, (3.86 
rise, Mee ;) eT ase 


To calculate «,2 we need ¢, to first order in «; we find that 
e, = , >, V inte Ot de ae . . . . (3.87) 
If we set k’=k+f and write 


bad =b,, (= BE Me ee 
iat k— 


then the first-order approximation to eqn. (3.81) can be written as: 


n(k’)(Vs0 = ae Kn View —)(k’, f) ‘ 


o(k, f) = V+ +— (Ho «= H.) 


(3.89) 
This is identical, apart from the definition of the n(k)’s, with eqn. (3. 34). 
If the exchange energy W, and the terms containing the factors Vy_,. 


and V{_,,, are neglected this equation becomes identical with that used 
by Bardeen (1937). We then find that 


k, f) V_s®n(k) : 
FOP Oy yeni PE 
: 2 (EL, — Ex) ck 


and lattice frequencies w, are given by the equation 


v(k » fyn(k) 
By — Eyag 
The only difference between these equations and the results we found in 
§3.3 is that the sums over occupation numbers are now over the usual 
unperturbed Fermi distribution (3.86). 
The lowest order one-phonon interaction term that can lead to real 
transitions is now given by the lowest order term in A this is 


of@=02— Vo oop oe (3.91) 


Alen - 5 [tn oh RiVin OdtVa; ‘ . = . (3.92) 


[f we introduce normal coordinates and momenta for the lattice degrees 
of freedom and express y¥,"" in terms of o(k,f), eqn. (3.88), we find that 
A, has the form 


v(k, f) ° 
ate eee toe nee 8) 
k k— 


This differs from the expression given in eqn. (3.59) in that the energy 
denominator has the Hartree-Fock energies £, rather than the Bloch 
energies «,. On the energy shell it is of order x. 

It is clear that by proceeding in this way we can calculate the 
Hartree-Fock approximation to every term in the Born—-Oppenheimer 
expansion. ‘The method therefore has the same limitations as the 
expansion itself, 
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$4. THERMODYNAMIC FUNCTIONS AND TRANSPORT COEFFICIENTS 


4.1. Thermodynamic Functions 

In this section we shall discuss the influence of electron—phonon 
interactions on the equilibrium properties of a metal. We first consider 
two properties, the cohesive energy and the normal mode frequencies 
that we believe are not very sensitive to electrons in the ‘danger zone’. 
The cohesive energy or, what is essentially the same thing, the eround 
state energy of the system can be written as 


6y=Hyt+ > Maw, + AZ, ke BAe tae 4a} 


where H, is the lowest electronic eigenvalue for the static lattice problem, 
Ddthw, is the zero point energy of the lattice vibrations and AZ is the 
correction term that we are interested in. This expression is correct 
no matter which of the methods we use; there is only a small difference 
in the definition of the frequencies in the various methods. In Nakajima’s 
method the lowest approximation to AH is 


,( _ Ny | V ei? 


AH=AENy=> (4.2) 


Ex — €k—#)2 — (he)? 
while in the Born—Oppenheimer and Hartree-Fock methods the lowest 
approximation gives 

(1 — _+)|V «|? 

(ep eo)" 

Both the expressions contain formally divergent sums due to vanishing 
energy denominators. We have seen that a possible description for such 
sums in the Nakajima method is to define them by taking a principal part. 
If this is done then we in fact find that AZ, ~ (m/M)"?> hw, and is very 
small indeed. The divergent sum in eqn. (4.3) can be defined by taking 
principal parts—though we do not now have a very good reason for doing 
so. If the sum is defined in this way then we again find that 


AE, ~ (m|M)#2 d Mica. 
i. 


A= A= Ss > . . - . (4.3) 


Thus both A#y and AH, are about 10-*ev in magnitude. While this is 
of course far beyond the accuracy with which we can hope to calculate 4, 
it is very large indeed compared with the energy difference between the 
normal and superconducting states of a metal which is of the order of 
10-7"ev. We mention this at this point to emphasize the enormous 
difficulty that would arise if we were to attempt to calculate these small 
energy differences from first principles. The divergences that arise in 
eqns. (4.2) and (4.3) of course come entirely from electrons in the danger 
zone. We could in fact argue that even though the methods cannot treat 
these electrons properly it is very unlikely that they contribute a larger 
amount to &, than all the rest of the electrons do. Consequently our 
numerical estimate of AZ is likely to be reliable. If this is so then we can 
conclude that it is in fact quite negligible. 
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We have already given formulae for the normal mode frequencies 
eqn. (3.44); we shall work in the framework of the simple Hartree 
approximation, eqn. (3.37). This leads to the following expression for 
the velocity of sound 


N\i2 1 
(Oh (=r) Sih? es 5 . 5 3 (4.4) 
M 4/ So 
where S, is defined by the equation 
k n({k 
Siceslin ee ee As ees, (4.5) 
Iflo0 k &k— €&k-#- MF kk Ek FE TOF 


In the Hartree-Fock method the phonon energy w, would be absent in 
the denominator. This leads to a very small difference in the velocity 
of sound in the two methods; see Appendix 3. We can therefore write 


§= lim Bye 

1 f1s0 €k — &k—F 
= lim pguabee Sl! 

; 1f1>0 Ex — €k—f 
This last expression for S, shows that it is only dependence on these terms 
in the sum for which e, and «, _, lie very close to the Fermi surface, and 
are in fact in the danger zone. On the other hand the general expression 
given by eqn. (3.44) for the normal mode frequencies shows that for most 
values of |f| the important terms in the sum for S(f) come from electrons 
lying. well away from the Fermi surface. We conclude that the normal 
mode frequencies are sensitive to electrons in the danger zone for small 
wave numbers. For large wave numbers they are not. For a spherical 
Fermi surface we can write 


(4.6) 


See n(k) — n(k —f) 
Ifls0 k Ex — €x_f 
where V(Z,,) is the density of states at the Fermi surface. This expression 
for Sy leads to the simple formula for the velocity of sound given by 
eqn. (3.39). The fact that this simple calculation leads. to results of the 
correct order of magnitude encourages one to believe that in this instance 
we may have given a reasonably accurate description of the motions of 
the electrons in the danger zone. A more sophisticated comparison of the 
normal mode frequencies with experiment has been carried out by Bardeen 
and Pines (1955). They conclude that the elastic constants calculated 
from the frequencies agree with the experimental values in the alkali 
metals to within about 20%. 
It has been pointed out a Kohn (1960) that the equation for omission 
implies that dw(f)/d|f] has singularities at certain values of |f|, for a free 
electron model we find that 


=N(EL 2. ee 


dev/OlF[ = — ba? /Q vie asi) 
3A 4ky?+|fP 1 | ake + IE 
= ~$w,2/O;V °— 2 5 8 
on ELK 4kylfP “la | ve 


and hence has a logarithmic singularity at ene - 
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If a more general expression for the single particle energies is used then 
Kohn has shown that the singularities in dw(f)/dl|f| can give information 
about the shape of the Fermi surface. These singularities have recently 
been detected in lead by Brockhouse (1961). The most prominent of them 
does occur at |f[=2k,. Lead is probably an especially favourable case 
because the ‘magnitude’ of the singularity clearly depends on the size 
of the matrix element and this is known to be large in lead. 

We now turn to two quantities that may be very sensitive to electrons 
near the Fermi surface; namely the specific heat and the spin para- 
magnetism. The lowest order corrections to the specific heat can be 
calculated quite easily by treating the correction terms, eqn. (3.45), by 
-perturbation theory. If we calculate the lowest order correction F° to 
the free energy in Nakajima’s method, then we arrive at the same 
expression as Buckingham and Schafroth (1954), namely 


Tes. M(1 — Hy_1)0" 


$< 4.9 
lex — €, «|? — (hws)? oe) 


Now F° is temperature dependent through the Fermi functions 7, in the 
numerator of the sum in eqn. (4.9). These imply that the temperature 
dependent part of the product 7,(1—%,_,) will be non-zero for transitions 
for which e, and «,_, lie within k7' of the Fermi surface. A rough estimate 
of the magnitude of the energy denominator for these transitions is 
(kT)? — (hw)?, where @ is an average value of #. At temperatures such 
that kT’<ho& we can therefore approximate F° as far as its temperature 
dependence 1s concerned by the expression 


cw yy Ml L= M4) |, I2 
E> aa ayes Mol Mike re oie octegtl aL O) 

It is clear that this term can be quite large because of the very small 
energy denominator. In fact a simple order of magnitude estimate of 
(4.10) indicates that it is of order C(&k7/H,,)?, where C is the conventional 
(Wilson 1953) electron—phonon coupling constant. Since C ranges in 
value from (4—1)H, it is clear that the correction to the temperature 
dependent part of F'© is quite comparable with the ordinary temperature 
dependent part of the unperturbed free energy which is of order of 
magnitude #,(kT/E,)?. A detailed calculation of the specific heat 
contribution from F° has been made by Buckingham and Schafroth 
(1954). They find that at low temperatures the specific heat has an 


extra term O\2 /1\ U8 
2) [ 
— — . . . Le 3 b a 11 
Sc=yeT" Ga ) (5) ky 


where yg7" is the usual Sommerfeld specific heat. At temperatures above 
the Debye temperature we can no longer argue that ¢,—«€,_, 18 very small 
compared with the average value of hw, The energy denominator will 
now be quite large on the average and we must expect that the extra 
contribution to the specific heat will be greatly reduced. This expectation 
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is confirmed by the calculations of Buckingham and Schafroth who find 
that for 7’> 6 the contribution (4.11) is reduced by a factor k6/H,~ 10~* 
and is consequently quite negligible. It is clear from this discussion 
that the very large contribution to the specific heat came from transitions 
in the region of the Fermi distribution that is thermally excited. That 
is a region of width k7’ about #,,; consequently if &7' is small we are led 
to very small energy denominators. The importance of this result lies 
in the fact that the effect comes entirely from electrons near the Fermi 
surface. If therefore it were possible to check it experimentally, then 
we should have a direct check on the validity of the method for this group 
of electrons. So far no experimental check has proved possible—one 
should note that C is independent of the mass / of the ions and so there 
is no dependence of ¢ on M. There are of course small corrections to the 
specific heat from the dependence of the normal mode frequencies on the 
distribution of electrons. These corrections have been calculated by — 
Chester and Houghton (1959) who find that they amount to only a few 
per cent. 

The effect of the correction terms to the free energy in the spin 
susceptibility can also be computed. This is perhaps of more interest 
as it is possible to measure the susceptibility over a wide range of 
temperatures. In contrast the ‘electronic’ specific heat of a metal is 
only directly observable at low temperatures when the lattice specific heat 
is negligible. 

In Appendix 2 we outline the calculation. We find a correction dy, 
to the temperature independent susceptibility yp», given by, 


ke 
byo~ O (F.) Xo: . - . . : . (4.12) 
whereas the temperature dependent part y,, has a correction 
8xp~O(C)xp (T’ <6), 
a ate tds) (8 S. 
~0 (lH) x ("> 6), es 
Ey 


where C' is again the electron-phonon coupling constant. These results 
are really very similar to those we found for the free energy. There the 
temperature independent part—the ground-stage energy—had a negligible 
correction while the temperature dependent part had a correction of 
exactly the same form as y,».. The most important conclusion to be drawn 
from this result for y is that the effects of the electron—phonon correction 
terms cannot be simply described in terms of a modified density of states 
for the electrons. For if this were so we must expect that y and C would 
be modified in the same way. The physical reason for the small modification 
of yy appears to be that the distortion of the Fermi distribution caused 
by the application of a magnetic field is really quite different from that 
caused by raising the temperature of the metal. For in the latter case 
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electrons are excited above the Fermi level and holes are left below; 
consequently transitions can take place quite freely among states within 
kT of Hy. When a magnetic field is applied electrons of one spin orientation 
are ‘excited’ to states above H, and electrons of the opposite spin 
orientation are ‘removed’ to states below #,. One if therefore tempted 
to say that a certain number of ‘electron-hole’ pairs have been created. 
This is not really correct however because all the excited electrons have 
one spin orientation and all the holes, the opposite. Since the electron— 
phonon interaction cannot change the spin of an electron no transitions 
are possible between the ‘electrons’ and ‘holes’ and y, is not appreciably 
modified by the electron-phonon interaction. Again there are other 
small corrections to x) from the normal mode frequencies. They amount 
to only a few per cent (Chester and Houghton 1959). It should perhaps 
be mentioned that they are temperature dependent; so a very accurate 
measurement might detect their presence. No direct calculation of the 
effects of electron—phonon interactions on the Knight shift have been 
made. We have little doubt that qualitatively the results will be very 
similar to those we have described for the spin susceptibility. 

It should be clear from this discussion that while many properties of a 
metal are really quite insensitive to the behaviour of the electrons near 
the Fermi sea, there are others which are really very sensitive indeed to 
the way we deal with those electrons. Any experiments that would give 
information about this very special group of electrons would be extremely 
valuable. 

In conclusion we should like to mention that it is likely that for 7’ > 6 
an entirely different approach to the problems we have just discussed 
can be found. For at such temperatures it is plausible that we can expand 
the partition function in powers of 0/7’. Since docM~'? this would be 
an expansion in powers of «? and would be closely related to the 
Born—Oppenheimer expansion. However, we suspect that it might have 
a wider range of validity than that expansion. The reason for this is that 
in the expansion of the partition function we would only require that 

¢/kT' <1 rather than that ¢/AH <1, where ¢ is a measure of the strength 
of the ‘perturbation’ we are poasiderins: In our case C would be the 
average value of the non-adiabatic terms C’,,,, in the equations of motion. 
The energy difference AH is the energy separation of the states we are 
interested in. Now we saw that we could satisfy the criterion C/AH <1 
except possibly for a group of states very close to H, for which AE is very 
small. However, if this criterion is replaced by C/k7' then this particular 
group of states no longer plays any special role. This is an example of a 
familiar situation in statistical mechanics. The difficulties that arise in 
Rayleigh-Schrédinger perturbation theory due to small energy 
denominators do not occur in the perturbation expansion of the free 
energy at a finite temperature (Schrafroth 1951). If we use our estimate 
of C,,,, then we require that 

ho 


ca 
ik ~ 
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Now v is close in magnitude to the velocity of sound, and /-'~a™ where 
a is the interatomic spacing. Consequently hil ~hea~1~ hw max ~ k8, 
and we then require that 

O/T <1 
as we guessed earlier. The low-order terms in this expansion have recently 
been obtained (Major 1960, unpublished) and the expansion appears to 
converge well for 7'> @. 


4.2. Transport Coffiecients 


In §3 we have discussed in detail how the various theoretical approaches 
handle the real one-phonon processes that we believe are responsible for 
the electrical and thermal resistance of a metal. Our conclusions in that 
section were (a) that the Nakajima method can be modified to handle these 
processes in a satisfactory manner, (b) the Born—Oppenheimer and 
Hartree-Fock methods lead directly to interaction terms in the 
~Hamiltonian which give rise to this type of transition. We showed in 
fact that in all three methods the interaction term on the energy shell 
could be written in the form 


> ta a,.*a,_ ¢v(k, f). i) 5 yeh key pee 
k,f wh f 
Since the matrix elements of p,/hw, are the same as those of q, this 
interaction term is very similar to the one that is introduced phenomeno- 
logically in the theory of transport coefficients (Wilson 1953). The only 
difference being that the effective matrix element v(k,f) depends on the 
wave vector k as well as on the difference of wave vectors between the 
initial and final states. If we make the simple Hartree approximation the 
v(k, f) is replaced by Bardeen’s (1937) matrix element v(f). A detailed 
discussion of this matrix element has been given by Ziman (1960). In $3 
we referred to two attempts by Hone (1960) and Bailyn (1960) to compute 
v(k, f) from the more general eqn. (3.42). They both compute averages of 
v(k, f) over all values of k. These average values turn out to be comparable 
in magnitude in v(f). The eqn. (3.42) is a three-dimensional integral 
equation in which one of the variables (the azimuthal angle) can be disposed 
of. ‘There would probably be little difficulty in solving it numerically. 
It would of course be very difficult to use the computed matrix element 
to calculate the transport coefficients. Nevertheless it would be of 
considerable interest to compare v(k,f) with the much simpler Bardeen 
matrix element o(f). At this stage we should point out that the calculation 
of both v(f) and v(k,f) becomes very much more difficult if we start to 
consider the detailed band structure of the metal; results can only really 
be obtaiend for very simple models (Ziman 1961). The qualitative 
agreement with experiment when the simple Bardeen matrix element is 
used is quite good; this gives us some confidence in the basic methods 
employed. The thermoelectric power of a metal does seem (Ziman 1960) 
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to be very’sensitive to the ‘shape’ of the matrix element and this property 
may therefore provide some additional information about the assumptions | 
on which the theory. is based. 

We conclude this section by pointing out that the effect of the diagonal 
self-energy terms in eqn. (3.45) on the transport coefficients has never 
been investigated. Since the electrons involved in the transport of energy 
and charge are again those very near to the Fermi surface these self-energy 
terms might give a significant modification to some of the transport 
coefficients. Perhaps the simplest quantities to consider in this connection 
are the anomalous skin effect (in the extreme anomalous limit) and the 
Wiedemann-Franz ratio. 


§ 5. SUPERCONDUCTIVITY 


The superconducting transition temperatures of mercury (Maxwell 
1950 a), tin (Maxwell 1950 b), lead (Olsen-Boér 1951) and thallium (Maxwell 
1951), depend on the mass of the ions in the lattice. The transition from 
the normal to superconducting state of these metals must therefore be 
influenced by the interactions of the electrons with the vibrating lattice. 
Quite independently of these experiments Frohlich (1950) suggested that 
the exchange of virtual phonons by the electrons could lower the energy 
of a metal and cause a transition of the superconducting state. It is 
tempting to generalize from the rather scanty data we have mentioned 
and to say that the transition temperatures of all superconductors will 
depend upon the mass J of the ions in the lattice. However it has recently 
been shown (Geballe et al. 1961) that the transition temperature 7, of 
ruthenium has at most only a very weak dependence on M. For the four 
elements we have mentioned earlier 7'\,ocM-!?; 7, for ruthenium shows 
less than 10% of the variation that would follow from this dependence. 
The most that can therefore be said with certainty at the present time 
is that the electron—phonon interaction plays some part in the transition 
to the superconducting state. It is by no means established that it plays 
any very general role. What is obviously required are more extensive 
measurements on a much wider range of metals to see to what extent the 
dependence of 7’, on the isotopic constitution is a general phenomenon. 

We do not intend to review the current theories of superconductivity 
(Bardeen et al. 1957, Bogolyubov ef al. 1959) in detail. The reason being 
that the current theories only require that there be some effective 
interaction between the electrons which is attractive close to the Fermi 
energy. The theories are not sufficiently well developed to determine 
whether this interaction must come via the phonons or could arise from 
another source. One can in fact legitimately regard the theories as models 
of the superconducting state in which it is postulated that an attractive 
interaction exists between electrons with energies close to the Fermi 
energy. The importance of this type of attractive interaction was first 
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realized by Cooper (1956). He showed that if two electrons moving above 
a Fermi sea of non-interacting electrons had an attractive interaction 
then they could form a ‘correlated state’ with an energy lower than the 
Fermi energy of the system. The motion of the two electrons in this state 
is correlated over very large distances (of the order 10~*cm) and there is 
an energy gap in the spectrum of the excited states of the system. It is 
now believed (Van Hove 1959) that it is a quite general feature of Fermi 
systems that the normal state, in which the electrons are pictured as 
filling a Fermi distribution up to a definite level, becomes unstable in the 
presence of an attractive interaction. We say ‘believed’ because there is 
no general theoretical proof of this assertion. If this is the case then we 
are immediately faced with the problem of explaining why only some 
(40%) of the metallic elements are superconductors. This is of course an 
extremely difficult question to answer. The usual reason put forward is 
that the Coulomb (screened or otherwise) interaction is certainly repulsive 
and must to some extent balance out any other attractive interaction. 

There is one other point that is worth discussing qualitatively. This is 
whether it is possible to decide which, if any, particular group of electrons 
is involved in the transition from normal to superconducting state. There 
are quite strong arguments to support the view that only a small fraction 
of all the electrons in the metal are appreciably affected by the transition. 
The strongest argument being the very small energy difference between 
the two states. This energy difference is typically about 1 part in 10-%. 
If all the electrons were involved in the transition this would imply 
a change of energy of about 10-7 ev per electron. While such a very 
small energy change for all the electrons cannot be ruled out it seems 
much more plausible to assume that a much smaller number is involved 
and that each suffers a comparatively large energy change. Now we have 
already given a qualitative argument (§2) that the bulk of the electrons 
probably follow the motion of the lattice in a simple adiabatic fashion. 
We estimated that only those electrons within hwmax of the Fermi energy 
are likely to be affected in a non-adiabatice way by the motion of the lattice. 
If only these electrons were involved in the transition then the energy 
change per electron would be about 2 y/h@ max X 10-7 ev~ 10- ev. 
Since this is of the order of k7'g where 7, is the transition temperature 
it seems to be a much more reasonable quantity than 10-7ev. We believe 
that the qualitative argument we gave in §2 does lend some support to 
the view that it is in fact electrons in this special region that are involved 
in the transition. 

_ If this conclusion is correct, then we can understand at once why 
quantities used as the cohesive energy, normal mode frequencies and the 
effective matrix element 1, are changed very little by the transition. 
For all of these quantities involve all the electrons in the Fermi distribution 
and only a small fraction (hw,,,./H,,) of these are involved in the transition. 
On the other hand the specific heat, and spin susceptibility should be 
considerably modified as these involve electrons in the region within 
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homax of Hy. It is of course quantities just like these that are found 
experimentally to be quite different in the two states. 

We must now consider whether any of the methods that have been 
described in §3 could provide a satisfactory basis for a theory of super- 
conductivity. As we have seen, what is probably required is an interaction 
between the electrons which is attractive in the neighbourhood of the 
Fermi energy. Unfortunately it is just in this region of energies that the 
theories are least reliable. Examination of eqns. (3.27), (3.65), which 
define the transformations S* and S?, shows at once that they become 
large when |e,—«_,|</iwmax. One cannot have much faith in the 
convergence of the perturbation expansions for these electrons. This 
point has been made more precisely by Bardeen (1956). 

In spite of this fact Bardeen et al. (1957) did attempt to base their 
original theory on the effective interaction derived by Nakajima (1953). 
This interaction, eqn. (3.45), becomes large and negative for transitions 
in which|e,— _,|</iwmax. In their original paper Bardeen et al. assumed 
that this was the most important part of the interaction and used it to 
develop the theory. However, since none of the detailed properties of the 
interaction was used in the subsequent development it is really left quite 
open whether this particular form of the interaction is vital to the theory. 
In other words only the more detailed quantitative features of the 
behaviour of superconductors can help us to determine the structure of 
the interaction terms in the Hamiltonian. The only attempt that has 
been made to correlate the interaction with experimental properties of 
superconductors is due to Pines (1958) and Morel (1959). 

Bogolyubov (1957) started from the usual phenomenological screened 
interaction (Frohlich 1950, Wilson 1953). As we have seen, there is very 
little justification for introducing such an interaction except for real 
- one-phonon processes. Indeed if this interaction is introduced then it 
is very difficult to introduce Coulomb forces between the electrons and 
treat them consistently ; for a ‘large part’ of these forces has already been 
taken into account in the screened interaction. It is in fact quite easy to 
make the mistake of screening the electron—phonon interaction ‘twice 
over’. As far as we are able to see the later work by Bogolyubov et al. 
(1959) on the Coulomb effects does contain this error. However, if his 
basic Hamiltonian is ‘reinterpreted’ in terms of bare rather than screened 
interactions, then his final results (eqn. (6.28), Chapter 6) for the 
renormalized frequencies and matrix element are identical with those we 
derived using the Hartree approximation. For this reason we do not think 
that his method of handling the electron—phonon interaction is any more 
convincing that that of Bardeen e¢ al. 

We conclude this section with some suggestions about how the current 
methods could be modified to provide a satisfactory basis for a theory of 
superconductivity. The major defect of the methods we have described 
is that they are all in some sense ‘perturbative’. This is clear from the 
fact that the final physical properties are expressed in terms of the Bloch 
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(or exchange) energies of the electrons. In particular the effective 
electron-electron interaction term contains these energies in the denomi- 
nator. Now if we believe that the excitation spectrum of a superconductor 
is quite different from the Bloch energy spectrum, then we might expect 
that the actual energies of the excitations should appear in such quantities 
as the velocity of sound or the electron—phonon matrix element. However, 
the excitation spectrum of a superconductor is not the solution of any naive — 
zero-order approximation. This suggests that a satisfactory theory can 
only be developed if the zero-order approximation is correctly chosen. 
A possible way of choosing the zero-order approximation would be to 
suppose that there is an attractive interaction between the electrons and 
to add and subtract this from the Hamiltonian. The zero-order 
approximation could be chosen to include this attractive interaction 
and we would then obtain a superconducting state as the solution to the 
lowest order equations. If now a canonical transformation similar to 
Nakajima’s is applied to the Hamiltonian, then the frequencies would be 
renormalized in terms of the excitation spectrum of the superconducting 
state. The attractive interaction that we have also subtracted from the 
Hamiltonain could be eliminted in a self-consistent manner by demanding 
that it be cancelled by the interaction term (eqn. (3.45)) that is generated 
by the canonical transformation. Ifa programme like this could be carried 
through, the physical quantities such as the frequencies and effective 
matrix elements would be expressed in terms of the actual excitation 
spectrum of the system. This would in turn imply that the energy 
denominators would never vanish if there is a gap in the excitation 
spectrum. An alternative scheme would be to transform the Hamiltonian 
in terms of quasi-particle operators (Valatin 1958, Bogolyubov et al. 
1959) before attempting to renormalize the bare frequencies and matrix 
element for the electron—phonon interaction. This is essentially the method 
used by Bogolyubov in his original treatment of the problem—except that 
he started with a screened electron—phonon interaction and omitted the 
Coulomb terms. As we have already mentioned one can ‘reinterpret’ 
his basic Hamiltonian in terms of bare interactions rather than screened 
ones. If this is done then the renormalized frequencies and matrix elements 
are indeed expressed in terms of the actual excitation energies of the 
system. This method does of course lead to a complicated set of equations 
that have to be solved self-consistently. So far only the simplest 
approximations have been made to this problem. These lead to the usual 
Hartree approximation (§3.4) to the matrix elements and frequencies. 

We conclude this section with some comments on the low-temperature 
specific heat of superconductors. It has recently been found (Bryant 
and Keesom 1960, Borse et al. 1960) that at very low temperatures 
T/T',<1/7, the specific heat of superconducting niobium and indium » 
does not tend towards the lattice specific heat of the normal state but 
to a value lower than this. The data can be fitted to a 7* curve which 
we can write in the form 
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where the term A,,7" is the normal lattice specific heat. The term — BT? 
is about 10-15% of this at the temperatures of interest. The ‘constant’ B 
may be temperature dependent, Bryant and Keesom suggest that 
B~we?, Tt would seem that both terms in eqn. (5.1) could be combined 
into a single term and one could then say that the lattice specific heat has 
altered by some 10%. This would imply a change in the velocity of sound 
of about 3%. Direct measurements by Chandrasekhar and Rayne (1961) 
show that the velocity of sound in indium changes by only 1 part in 10-5, 
so that this simple suggestion cannot be correct. An alternative suggestion 
has been put forward by Daunt and Olsen (1961). They point out that the 
free energy of the system always contains the ‘zero-point energy’ of the 
lattice modes >,4iw,. Now we have seen that the frequencies o/(f,s) 
depend on the electronic configuration through the presence of the factor 
mk) in eqn. (3.35). In §4.2 we pointed out that because of this there is a 
contribution to the specific heat of a few per cent in the normal state of a 
metal. Daunt and Olsen suggest that a similar effect should arise in the 
superconducting state and would lead to an additional term in the specific 
heat. This term really arises because of the coupling between the motion 
of the ions and electrons. It is not therefore a contribution arising from an 
additional mode of motion that has been neglected—for this reason it can 
have a negative sign. No actual calculation of the size or sign of this 
contribution has been published. The effect certainly exists and the 
calculation is perfectly straightforward. It is perhaps rather strange to 
talk in terms of a temperature dependent zero-point energy. One normally 
talks in terms of a zero-point energy that is associated with the ground 
state of the entire metal. However, since the frequencies w(f,s) depend 
on the electronic state of the system the ‘zero-point energy’ $hw, will 
also depend on the electronic state. Consequently the true zero-point 
energy must be defined in terms of the set of frequencies which are defined 
in terms of the ground state of the electronic system. It is not however this 
zero-point energy that really occurs in the free energy but rather the 
statistical average of >,t/w, over the electronic states that are occupied 
at any given temperature. This implies that >,}hw, will be slightly 
temperature dependent; it is for this reason that we speak, somewhat 
loosely, of a temperature dependent zero-point energy. 
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APPENDIX 1 


THe MAGNITUDE OF THE NON-ADIABATIX Matric ELEMENTS 


The matrix element C,,,,, is defined (eqn. (2.5)) as the sum of two terms 
A,, and B,,,. We estimate the magnitude of A,,,. Exactly the same 
methods can be applied to B,,,,.. Now 


He 
Any =~ 57% [Ua* Vajlae OV oie AES 


The ‘matrix element’ is of course really an operator because only the 

electron coordinates are integrated over. What we want to do is to 

estimate the magnitude of the product of A,,,, and a typical zero-order 

wave function yn» (see eqn. (2.1)). The velocity operator of the jth ion 
= — (th/M) Vp; So we can write 


Age Xn =th Sn" VaiVn’ i O;Xmn- . i (A 1.2) 
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When the operator v; acts on x,,,° it can be replaced by 6, where @ is the 
average velocity of the ions in the lattice. A more unfavourable estimate 
would be to replace v; by the maximum velocity of the ions, this would 
lead to only a small numerical change in our estimate. We can therefore 
write for the magnitude of A,,,,: 


Ann S401 S[px* Vande]... (AL) 
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The Schrédinger equation for %,,, is 
Heh: =H nn: . . . . (A i 4) 
If we apply the operator > Vie ; to both sides of this equation, multiply 
by #,,* on the left and iapeeene over the electron coordinates, we find that 


(V R; Vint)an’ 
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where Vint is the interaction potential of the electrons with the ions. If — 
the ions suffer displacements whose average magnitude is 1, then the 


distortion dy, in x, is given by 
(= Va,in) | 
bb, =1 avatar he a n sah ee 
Pn > iE, bn ce AST.) 
For an appreciable distortion we require that / satisfies 
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Combining this relation with eqn. (A 1.5) we find that 


1 
X [Yor Vaghar@r | ~ 7, ee as) 


where / is the average distance the ions have to be moved to produce an 
appreciable distortion in %,,. 


ae N DEX 2 


THE SPIN SUSCEPTIBILITY 


In a uniform magnetic field H we can write the lowest order correction 
to the free energy in the form: 


+(1—n,-¢7) mM (l—™%_+7) 
iff ce SSL EL A a ee 
>, a ae > GHG ( ) 
Here 
me = [exp [ee —n*/kT] +1], tee (499) 
where 
n*=n(H) + pl eee re GA 3) 


and 7(H) is the chemical potential of the system in the magnetic field. 
It is determined by the equation: 


N => (mt +00); Reape ai (A 2.4) 
‘ k 
a straightforward calculation gives 
1 p22 
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This expression is accurate to second order in H. The correction dy, 
to the susceptibility is given by 


o2Fe 
Sp He eee So eg tele) 
ae GDh, 
sored otis + (SR) ii pees (A977) 
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after some simple manipulations we find that 


GHW) + ae | (A 2.8) 


On? On? 
In this last equation F°(0) is the correction to the free energy at zero 
temperature in zero field, F,°(0) is the correction to the free energy to 
order (k7'/n)? in zero field and 7 is the chemical potential in zero field. 
Now the dependence of F° on 7 is given by, 


Fo=NnO| eyns(2m)n)-+ (=). (A 2.9) 
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Here n, and n, are two pure numbers of order unity, c, is the velocity of 
sound and C is the electron—phonon coupling constant. The first term in 
the expression is of the order of k@ per particle—as we mentioned in § 4.1. 
The second temperature dependent term leads to the large correction to 
the specific heat that was also discussed in §4.1. If we substitute this 
expression in (A 2.9) then we find that 


_ (ke 
5xo=Pite(—) + Pox’, +? yeaa 


where p, and p, are two numbers of order unity, xy) is the temperature 
independent Pauli susceptibility and y, is the temperature dependent 
term in the susceptibility. This is the expression on which we based our 
discussion in § 4.1. 


APPENDIX 3 
THE Function S(f). 


In this appendix we shall calculate the function 


ie raeres 


for a free electron model. We shall show that it is negligibly different 
from the function 


(A3.1) 
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The z interaction is trivial and we find: 
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The remaining integration is elementary and we finally find that 
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Now 
a=2melh~ a Kp~ LO str : 


consequently, for all values of f, « is negligible. 
We find that in the limit | f |+0, 


3N 
S(f) = —- (1+ 0(k6/Z)?) Pipe fistaL pt (ao) 
Di 
whereas the function S°(f) tends to the limit: 
3N 
oie (A 3.10) 


The correction terms are therefore of order 10-4. 


Since « is negligible for all |f| we can set it equal to zero in our expression 
for S(f). We then find that 


BN [1 4ep®@—f2) [Qe tf 
De \hoeberacs fy, og | etd |. eA 


The largest value of f is 2x, and for this value 
bless aire es rs ok (A312) 


This is exactly one-half of its value in the limit |f|+0. It is a simple 
matter to show that 


Sif) = 


S(2kp) < S(f) < S(0) Rey (A 3.13) 


and is a very smoothly varying function of f. 
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Crystal Stability and the Theory of Ferroelectricity 
Part II. Piezoelectric Crystals 


By W. CocHran 
Crystallographic Laboratory, Cavendish Laboratory, Cambridge 


ABSTRACT 


In Part I the dielectric properties of cubic crystals with atoms in special 
positions were discussed, with particular reference to crystals which become 
unstable against a transverse optic mode of vibration of long wavelength. 
These considerations are now extended to crystals of general symmetry, and it 
is shown that the anomalous temperature dependence of certain dielectric, 
piezoelectric, and elastic constants in the paraelectric phase of ferroelectric 
crystals can be described and correlated in a very natural way in terms of 
crystal stability. A formula relating the dispersion frequency of the anoma- 
lous dielectric constant to the spontaneous polarization and the displacements 
of the atoms in a ferroelectric transition, previously derived for barium 
titanate, is generalized. The transitions in KH,PO, and in NH,H,PO,, and 
the ordering of the hydrogen atoms in these crystals, are discussed qualita- 

‘ tively in terms of the stability of normal modes of vibration. 


§ 1. INTRODUCTION 


In Part I (Cochran 1960) a ferroelectric transition in a crystal was considered 
in terms of the normal modes of vibration. Attention was confined to 
crystals which are cubic in the paraelectric phase and have atoms in special 
positions in the unit cell. This approach provided an interpretation of 
certain of the properties of crystals of the barium titanate type; in 
particular it proved possible to correlate the dielectric properties with 
the displacements of the atoms at a transition and with the frequencies 
of certain normal modes of vibration of long wavelength. In Part I it 
was tacitly assumed that in the cubic phase each atom occupied a centre 
of symmetry in the crystal, which was therefore not piezoelectric. The 
elastic constants were not derived in terms of the parameters used to 
specify the short-range and the long-range (Coulomb) interactions; this 
could readily have been done and it would have been found that no 
anomalous variation with temperature was predicted. The piezoelectric 
_and elastic constants of the ferroelectric phases were not considered. In 
this paper these considerations are extended to include crystals which are 
piezoelectric in the paraelectric phase, and which are not restricted as to 
structure or space group symmetry. It is however assumed in the 
meantime that we have to deal with a perfect crystal, that is, one which 
does not exhibit the kind of positional disorder found in KH,PO,. It is 
found that a lattice dynamics treatment correlates satisfactorily the 
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properties considered in Part I, and also the anomalous temperature 
dependence of certain elastic constants and piezoelectric constants in the 
paraelectric phase. While many of the conclusions arrived at are no 
different from those given by the phenomenological theory of Mueller 
(1940) and of Cady (1946), it is of some interest to show that they can be 
reached by a more atomistic approach. The relations between macroscopic 
properties established by the thermodynamic theory and the fact that 
these relations are confirmed by experiment are sometimes summarized 
in the statement that ‘the piezoelectric and elastic anomalies can be 
considered to be consequences of the dielectric anomaly”. Our object 
is to show that all three may be given a more fundamental interpretation 
in terms of crystal stability. 

In §2 we derive general equations which govern the lattice dynamics 
of a crystal in which the atoms are polarizable. By way of illustration, 
§3 deals with the dielectric and elastic properties of the simplest type of 
crystal which exhibits piezoelectricity, namely the zinc-blende type. It 
is shown that a transition which exhibits all the characteristic features 
of a ferroelectric transition is in principle possible even in this type of 
structure. The discussion is extended to include crystals of general 
symmetry, in§4. The way in which a lattice dynamics treatment must be 
extended to take into account special features of KH,PO, and isomorphous 
crystals is discussed qualitatively in the final section. 


§ 2. Larrick DyNAMIcs OF A CRYSTAL WITH POLARIZABLE IONS 


The treatment of this problem is based on an extension of the 
Born—von Karman theory made by Woods e¢ al. (1960) which in turn 
can be related to earlier work by Mashkevich and Tolpygo (1957). It is 
assumed that the energy perturbation in a crystal in which the atoms are 
given small displacements is a quadratic function of the nuclear displace- 
ments and the electronic dipole moments. The most general expression 
of this form is the following: 

* @=-}> { Y [O,y (UK, UK')u, (Kyu, (UB) 
Ke (UV Ky 
+ YO, OK, UR’ Ju, (UK) p,(UK’) 
+ Y p10, UK", LK) p(1K)u,(UVK’) 
+ ¥gA¥ p71, (9K, UK")p,(IK)p,(UK’)] 


~ cx-1p, AIK) + (p_(LK) + Zytta(IK))B AUK) } nee tio 


Here u(/K) denotes the displacement of the Kth atom in the Ith cell, 
p(/K) its electronic dipole moment, a, its electronic polarizability and Z, 
its ionic charge (taking +4-8 x 10-e.s.u. as the unit of charge), while 
E(/K) is the field due to the remainder of the crystal. In the Born theory 
of ionic crystals (Born and Huang 1954) the tensors ®” and ®©) do not 
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appear, and in much of the theory hitherto it has also been assumed 
_ that «,=0. The additional terms involving ®® and ®) in eqn. (2.1) 
take account of short-range forces between polarized atoms. They can 
be interpreted witltin the framework of the Born—von Karman theory 
in terms of a shell model in which the outer electrons are represented by 
a shell carrying a charge Y,. This shell exerts short-range forces on both 
cores and shells of peer pouring atoms, and p(lK)= Y ,w(lL) where w(lK) 
is the relative displacement of shell and core. It should be emphasized 
that this model does no more than provide a convenient nomenclature 
and a link with the Born—von Karman theory; the eqns. (2.11), at which 
we eventually arrive, can be written so as to involve only atomic displace- 
ments and dipole moments, without mention of cores and shells. The 
equations of motion to be used in conjunction with eqn. (2.1) are 

—00, 0® 


mit (IK) = sa» 0= a) Nees (222) 


(Mashkevich and Tolpygo 1957). These lead to 
Mrt(LK)= > (O,,°(K,UK'\u,(UK') 
4 @ PUK UK ww (UK) + pH (UK), 
0O= > (0,,P UK’, K)u,VE)+0,,O(K,UK'\w,(VK’)) 
yen en OK) ee ea. oe (2.8) 
These may now be written in terms of travelling waves by introducing 
u(lK) =U(K) exp i(q. r(IK) —t), | 
w(K) =W(K exp 2(q . r(LK) — at), 
E(l/K)=&@(K) expi(q.r(lK) —at), | 


(2.4) 


where ¢=27—~wavelength and w = 2n x frequency of a normal mode of 
vibration. The equilibrium position of an atom is 
r(UK)=r(l)+r(K), SA tek es ose 2D) 
where r(/) is a lattice translation. 
We also define the ‘coefficients’ of the short-range interaction by 
R,,(KK’)= — > Op, (1K, VK") exp iq. (r(UK") — r(LK)) 
_ 


and similarly for 
TEE AIRS ceri | KC); aetovesat (2.6) 


while the corresponding coefficients of the Coulomb interaction are defined 
by 


C,,( KK’) = —expiq. (r(K’)— r(K)) lim 
ror (K)—r(K’) 


«1 Daeg Poe “iO au) Pia 1} ies eo) 


(Kellermann 1940). 
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The amplitude of the electric field is then found to be given by 
8 (K)=— 3 Cyy(KE'\ZpU (K+ ¥ pW, (K')). (2.8) 
Ky 


and the equations of motion reduce to 
wm ,U(K)= > (Rz,(KE)+2Z SOAS OR AUG IE AG) 
a pa WKK')+ ZV pC y(KK’))W (kK), 
= xis (K'K)+¥yZ,0,,(KK'))U,(K’) 


ye 


a Le A(KK')+ V_¥ p0q,(KK'))W(K’), (2.9) 


where 

oS py( KK’) = Sqy(KK') + 82,8 e7-%R TY x. Ree ee) 
The eqns. (2.9) apply as they stand for crystals which are diagonally cubic. 
For crystals of lower symmetry the coefficients C,,(KK’), (eqn. (2.7)), 
and o%,,(KK’) (eqns. (2.6) and (2.10)) must be re- -defined for K=K’. 
The details are however not important for our present purpose and will 
be given elsewhere. 

Equations (2.9) can be written in matrix notation as 


w*mgU.=(R+ ZaCZqa)Uet+ (T+ ZaCYa)We, 2 11) 
O=(T*+ YaCZaUc+(S+ YaCYa) We, 


where T is the transposed matrix obtained by interchanging rows and 
columns of T. As in Part I, subscript d denotes a diagonal and c a column 
matrix. Other matrices such as R are square matrices of order 3n x 3n, 
where n is the number of atoms per unit cell. These equations are nearly 
identical in form with those given in Part I (p. 406), they are, however, 
much more general in that they now apply to a crystal of any symmetry 
and for any value of the wave vector q. The corresponding equations 
of Part I involved matrices of order n x n, they applied only for q>0 and 
to purely transverse modes of vibration. A simplifying assumption made 


in Part I also resulted in mean in the present notation. 

It has now been shown by Cowley (196la) that eqns. (2.11), which 
have been derived here from the potential function (2.1), can also be 
derived by a purely wave-mechanical treatment of the motions of the 
electrons and nuclei, provided that quadrupole and higher moments of 
the electron distribution are neglected. 


§ 3. DIELECTRIC AND Exastic PROPERTIES OF THE ZINC BLENDE 
STRUCTURE 
In Part Lit was found useful to discuss the lattice dynamics of a diatomic 
crystal of the sodium chloride or of the caesium chloride type. It was 
shown that with a suitable choice of the parameters governing the lattice 
dynamics of such a crystal, it would exhibit dielectric and other properties 
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closely paralleling those of barium titanate. In this section we show that 
the same choice of parameters for a crystal of the zinc-blende type will 
cause it to exhibit properties typical of the paraelectric phase of a piezo- 
electric crystal approaching a ferroelectric transition. A crystal of this 
class has three elastic constants C,,, Cy, and C,,. It can be shown from the 
theory of elasticity, and also by the methods used here, that O,,=0,,°, 
Cyg® =C yy? but Cyy"4Cy?. There is only one piezoelectric constant Cra 
The clamped and free dielectric constants ¢,” and ¢,’ differ, there is 
however only one dielectric constant <e at optical frequencies. The two 
atoms per a) unit cell will be identified by K=1 at the origin, 
K=2 at (—4, —4, —4)ry where 2r, is the (cubic) unit cell dimension. 
When the wave-vector q is in a direction such as [100] the 6 x 6 matrices 
of eqn. (2.11) can each be separated into three 2x2 matrices which 
separately describe the longitudinal and two transverse modes, for each 
of which the form of the equations of motion is still as given by eqn. (2.11). 

The macroscopic properties of the crystal can be derived quite simply, 
from these ‘one-dimensional’ equations. The dielectric properties follow 
from equations already given in Part I. The static polarizability is given 
by ags=aj+oe, where ae is the electronic polarizability (of one unit cell) 
and the ionic polarizability is given by 


Ret Zhe Ain Cees (351) 


In this equation the effective charge Z,’, and the effective short-range 
force constant Ry’ between Bravais arrays of atoms, are as defined in 
Part I. The folam ped) dielectric constant is then given by the Clausius— 
Mossotti relation 
Bs isl 22) (3.2) 
1 — 47rag/(3v) 
where v= 2r,* is the volume of the unit cell. A similar relation gives the 
high-frequency dielectric constant €e in terms of we, which is in turn related 
to atomic parameters as in Part I. The (circular) frequencies w;, and wy 
of longitudinal and transverse optic modes for which the wave-vector q 
tends to zero are given by 


[Low 2 = bg + 82(€e+2)(Z, ’)P/(9vee), 3 
eR — 41(€e+2)(Z, ‘2/(9v v), } ere, iis 


again as in Part I. 

The elastic constants can be derived by considering the limiting value 
of the ratio w/q for acoustic modes propagating in different symmetry 
directions in the crystal. To obtain explicit expressions for the elastic 
constants some assumption must be made about the short-range force 
constants between the atoms. We shall make the simplest one, namely 
that there are non-Coulomb forces only between nearest-neighbour atoms 
in the crystal. There are then only two independent force constants, 
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®,,.°(12), O,,%(12) and of course corresponding quantities with super- 
scripts (Z’) and (8) (eqn. (2.1)). Comparatively simple expressions for Cy, 
and for e,, are not obtained unless it is assumed that 

D,(12)/Ogq(12) = Dp yP(12)/ zg (12). cater, (aca) 
It was found necessary to make this assumption in order to account 
for the measured elastic constants of germanium (Cochran 1959), it is 
therefore a reasonable one in the present context and is included in 
what follows in this section. The elastic constants may be more con- 
veniently expressed in terms of the ratio (3.4), for which we shall write y, 
and a parameter 
Ry =49,,(12) Oa Tae Mee att 
(R, is also an element of the matrix R for q=0). The elastic constants 
involved, for acoustic waves travelling in different directions, are shown 
in the following table. 


Direction of 


Meet : Direction of Elastic 
BASLE OI Mie Sos electric polarization constant 
100 100 — Cy 


110 110 = : 
100 001 010 Cu 
110 001 D 


The explanation in lattice dynamics terms of the different electrical 
conditions which result in the appearance of C,," in one instance and of 
Cy,” in the other is as follows. Although the corresponding modes are 
described in elasticity theory as transverse, the polarization, in the sense 
of the motion of the atoms, is in fact elliptic with the ratio of minor to 
major axis proportional to wave-number qg. The axes of the ellipse are 
oriented as shown in fig. 1. Atoms K=1 and K=2 describe ellipses in 


Fig. 1 


‘| 


ip i 


Displacement and direction of motion of successive planes of atoms near a node 
of an acoustic wave in the zine blende structure. When q is parallel to 
[100] the minor axis of each ellipse is parallel to [010] but when q is 
parallel to [110] the minor axis is also parallel to [110]. The two types 
of dot denote the two types of atom, ‘ 


Crystal Stability and the Theory of Ferroelectricity 407 


opposite senses and since they are oppositely charged a wave of electrical 
polarization results which in one instance (third line of the table) is 
transverse but in the other (fourth line of the table) is longitudinal. The 
former gives no macroscopic electric field, the latter produces a macroscopic 
field of amplitude —47Y, where F is the amplitude of the polarization 
wave. A detailed derivation of the elastic constants will not be given here, 
the method is similar in principle to that used to derive the elastic constants 
of germanium (Cochran 1959). The correctness of the results has subse- 
quently been checked against the more general results obtained by Cowley 
(1961b). The elastic constants C,, and C,, are found to be independent 
of the electronic polarizability and are given by 


1 
C= 7, (Bo/8 +. 0-18 Zi 10), . ° . 5 e (3.6) 
0 


Cra = — (Ro(2y—1)/8—1-32 22/0). a3 7) 
0 


The expressions for C,," and C,,? are considerably simplified by writing 
them so as to involve ¢,”, thus: 


1 
Osa? = — (Ro(1 = y*)/8 + 2-962,2/0—2a( es" +2)(1-20— }y)°Z2/(90)), (3.8) 
0 
; 1 
Oe "I (Ry(1 — y?)/8 — 1-57Z,?/v + 4r(es” + 2)(0-60 + dy)?Z,?/(9veEs”)). (3.9) 
0 
The piezoelectric constant ¢,, can be derived in a similar way by 
considering an acoustic wave propagating in the [100] direction, with the 
atomic displacements along [011]. In this case the electric polarization 
is also along [011] and arises from the fact that the atoms K =1 and K =2 
do not move in phase. The ratio of polarization to strain gives the piezo- 
electric constant ¢,, since the macroscopic field is zero. The result obtained 
is 
C14 = (Qyro(Es” + 2)/(3v))(5-03as/v — Fy), Se eee) 
where as and «,” are related by eqn. (3.2). 
We now have rather more than the necessary number of relations 
between atomic parameters and macroscopic constants. The latter 
should be related by the following results 


wy? /wp? = €s"/€e aera epa et e311) 
(Lyddane et al. 1941) and 


Cue == Opies = Are, ,"/ €a” . eit ins. . . (3.12) : 


(Cady 1946). The first of these is known to be satisfied for this model of 
a crystal (Woods et al. 1960), and the second is also found to hold, which 
provides a useful check on the derivation of eqns. (3.8), (3.9) and (3.10). 
Other macroscopic constants could be calculated directly in terms of 
atomic parameters, but as in the above two instances, they cannot differ 


P.M.S. ae 
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from those given by the relations of the phenomenological theory. It is 
therefore simpler to obtain ¢s’, the free dielectric constant, from the 


relation 
“PY = €;°C,°/Gy" . . . . . (3.13) 


(Cady 1946), since all quantities on the right-hand side of (3.13) have 
already been calculated. 


Fig. 2 


Te To 


The variation with temperature of certain quantities when the frequency wy 
varies as shown in the lowest graph, for the zine blende structure. 


Now let us suppose, exactly as in §2 of Part I, that the force constants 
happen to be such that wp (eqn. (3.3)) is close to zero, but other frequencies 
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are comparatively high, and furthermore that the force constants and/or 
unit cell volume, etc., vary linearly with temperature so that 


L—4a(€o + 2)(Z,')*/(9vR, ) = pon?/R, =B(T =T"), (3.14) 


where f is a temperature coefficient of the same order of magnitude as 
the volume coefficient of expansion. The temperature 7',” is then that 
at which the crystal would become unstable against a transverse optic 
mode for which the wave-number is near zero. As in Part I, ¢,” will then 
follow a Curie-Weiss law with the Curie constant 


6 =(€e+2)/(478). octamng lactis fag NE) 


It also follows from the expressions derived above for the various 
other macroscopic quantities that their temperature variations will be 
as shown qualitatively by the full lines in fig. 2. In terms of crystal 
stability, what happens is as follows. Before the crystal can become 
unstable against a transverse optic mode, it becomes unstable against a 
transverse acoustic mode, since the expression for C,," contains a term 
which depends on ¢«s”. This term causes C,,® to vanish at a temperature 
T’.’, greater than 7'.”, and given approximately by 


8n2GZ,2(1-20 — dy)? 


TT , Ens Tt u a 
ORGS 9ur Or? 


(3.16) 


The free dielectric constant therefore becomes infinite at this temperature 
(eqn. (3.13)) and the crystal shears spontaneously about a crystallographic 
axis, while positive and negative ions move in opposite directions along 
this axis (see the third line of the table). The distinction between C? 
and C® can only be made for waves whose wavelength is small compared 
with the length of the crystal itself; the effective elastic constant for 
longer waves depends on the electrical boundary conditions. It follows 
that the crystal must become unstable against a mode whose wavelength 
is smaller than the crystal, which will therefore develop domains of alternate 
direction of shear and of polarization. These conclusions, and the features 
of fig. 2 shown by the full lines, follow from the theory of lattice dynamics, 
and one assumption, eqn. (3.14). 

One further assumption is necessary if the crystal is not to become 
completely unstable, and the same assumption determines the properties 
of the polarized phase. It is that the short-range force per atom between 
Bravais arrays displaced a comparatively large distance wu is anharmonic, 


F(u)=R,u+ Bui + Bu +.... a ae ee a 


The relation between applied field Z, and polarization P (tor the cubically 
clamped crystal) can then be shown to be 
By =P(0—T")|6 +P BPF[(Zy")8 + OB P5|(Zy")+..., (3.18) 


2G2 
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where the apparent charge Z,” is defined by 
AU CEA hl Sear ca 


(see Part I). The assumption (3.17) enables the curves of fig. 2 to be 
continued as shown qualitatively by the dotted lines. 

A crystal having the zinc-blende structure may in principle become 
unstable against other modes of vibration than the particular transverse 
optic mode considered above. Suppose an instability occurs, as a result 
of a linear variation of force constants with temperature, for a transverse 
mode for which q is at the boundary of the first Brillouin zone in the [001] 
direction, that is, q=(0,0,7/7)). In this mode the atoms move as shown 
in fig. 3 and such an instability would result in an antiferroelectric transition. 


Fig. 3 


Direction of displacement of the atoms in the zine blende structure in a trans- 
verse mode for which q=(0, 0, 7/79). The unit cell dimension perpen- 
dicular to the plane of the diagram is 279, and atoms at heights 0, +, 4and ? 
are shown by different amounts of shading. Atoms of type 1 are at 
heights 0 and 4, of type 2 at heights } and 8. 


Calling the frequency of this mode w,, the transition will take place at a 
temperature 7', for which w,~0. If 7, >7'.", in other words if the force 
constants are such that wy>w,, the crystal will exhibit no large anomaly 
in ¢,” or in C44". If however wy and wy, are simultaneously low, but the 
latter is the lower, then the variation of es” and of C,,® will not be very 
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different from that shown in fig. 2, except that C,," will not become zero 
unless 7’, <T'.’. This particular mode has been discussed because it appears 
to be closely analogous to that against which NH,H,PO, becomes unstable. 
As shown in fig. 3, the zinc-blende structure can be referred to a tetragonal 
unit cell with a= 4/2r,, c= 2ry and containing two each of atoms of type 1 
and type 2. This unit cell is preserved after the transition, but is not then 
body-centred. 

In the following section we consider how some of these results may be 
extended to crystals of.general symmetry. 


§ 4. GENERAL THEORY OF DIELECTRIC AND ELASTIC PROPERTIES 


The theory of the dielectric and elastic properties of crystals has been 
worked out by Born and Huang (1954) in terms of lattice dynamics for 
the rigid ion model. Their theory has been extended by Cowley (1961 b) 
to include the effects of polarizability; the results are then rather 
complicated and their significance for ferroelectric crystals can be more 
readily appreciated by discussing the rigid ion model, where the same 
principles apply but the equations simplify considerably. For the rigid 
ion model eqn. (2.11) becomes simply 


w?mgU,=(R+ ZaCZa)Ue. Pah oe (42) 


It follows from the work of Born and Huang that for the discussion of 
modes of small wave-number it is important to separate the matrix 
R+ ZqCZq into a part which is regular in q and a part which is not. 
Equation (4.1) may then be written 


w?mqU.= MU,— ZaEc, . . . e (4.2) 


where E, is a column matrix whose elements give the amplitude of the 
macroscopic electric field at each atomic site. (By restricting attention 
to transverse modes in Part I this term was excluded.) The matrix M 
may be expanded in powers of q: 


M=M+i > M0, +4 > M,,24,9- Do ORS 
y ya 


(In order to facilitate comparison with the results of Born and Huang, 
aBy will be used to denote direction instead of xyz of earlier sections and 
Part I.) 

The matrix M© is of order 3n x 3n and is singular. .An effective inverse 
of M® may be defined by constructing a matrix H® whose elements satisfy 
the relations 

H,,(KK’)=0 for K=1 or -K’=1, 
> A,,O(KK')M,,(K'K")=5,,d¢5 for K, Kk’, K"=2,3...n. (4.4) 


The following: results may now be obtained by specializing Cowley’s 
(1961 b) results to the rigid ion model, or by expressing those of Born and 
Huang (1954) in matrix notation. 
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The elastic constants are given by 
Coypan = [af, yA] ci [By; ar] as [BA, ay | oe (ay, BA) 


where 
(9B, Al= 5 5 (Mas) 
and 
(ay, BA) = — = > (M,°HOM,).p ee Ley 


the piezoelectric constants are given by 


if ~~ 
€,07= — 5 ©, (MOH Za)ag yal eenee 


and the elements of the (clamped) static susceptibility tensor by 


; | 
Xoe= 7% D,(ZaHZa)u: 9» + + + + (47) 


Each of the above summations is over all matrix elements with certain 
subscripts held constant, e.g. 


> (X)zg denotes > X,,(KK’). 
KB’ KK’ 


Next consider the solutions w,; of the equations 
w2mqUe = MU. Sh, i 


Three of these frequencies (j7=1, 2,3) are zero. The remaining (3n—3) 
frequencies are the dielectric dispersion frequencies. These frequencies 
are also the frequencies of modes of vibration for which the wavevector 
is perpendicular to an axis of the crystal and for which the electrical 
polarization is transverse, so that the macroscopic field is zero. They are 
therefore the infra-red absorption frequencies for special directions of the 
wavevector (Cochran and Cowley, to be published). Finally it may be 
shown that if A(j) is the amplitude of the polarization wave associated 
with the mode whose frequency is w;, and U(Kj) the corresponding 
amplitude of displacement of the Kth atom, the susceptibility at any 
frequency w can be written 


hy P(j)P pl J) 
Xap(@) =) > o;(w;* im w*) (4.9) 
where 
oj= > m,|U,(Kj)/?. Ayreon Mayet hiki) 
Ky 


The absolute magnitudes of the U,(Kj) are of course arbitrary; it is 
convenient to introduce normalized values U (Aj) such that 


> mglO(KHP= 1. at, ieee te 
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To appreciate the significance of these known results, for ferroelectric 
crystals one further result is required. It may be shown (see, for example, 
Born 1942) that if e,(Kj) are the eigenvectors and w;" the eigenvalues of 
a matrix D appearing in the equations 


este, ) => DiCKK jen K'), 300 a) 7.9.5 (4512) 
then ‘ne 
Do te, (Kije(Kj)= (DP) (KK), 1, 94,18) 


where the right-hand side of eqn. (4.13) is an element of the matrix 
(D raised to the power p). Applying this result to eqn. (4.8) with p= —1, 
and remembering that H® is the effective inverse of M®, one finds that 


: 3n A a 
| = ope CM AG)U (Kye eOUK KR \ ee 5.) eo 8 lo (4s14) 
I= 


This result enables the elements of the matrix H® to be expressed in terms 
of the dielectric dispersion frequencies. Suppose we now have to deal 
with a crystal which is nearly unstable against an optic mode of vibration 
for which g+0, and for which the electrical polarization is transverse. 
We shall refer to this as the ferroelectric mode and denote it by j=F. 
Furthermore let us assume that the frequency w, varies with temperature 
according to the law . . 

w pC (T—T.”), ters vy, 4, LD) 


which, by analogy with the simple structures already considered, is to 
be expected if the force constants and unit cell volume vary linearly 
with temperature. If these quantities have a comparatively small 
temperature coefficient (~10~ per °K, say) there will be no rapid variation 
of individual elements of the matrix M®, and the values of the U (Kj) 
will remain nearly constant, even for j=. From eqn. (4.14) however we 
see that the term w,* will cause certain elements of the matrix H® to 
vary rapidly with temperature, in fact to follow a Curie-Weiss law when 
eqn. (4.15) is satisfied. Since the matrix H® appears in eqns. (4.5), (4.6) 
and (4.7) it follows that certain elements of the (ay, BA), eg ,, and y,, tensors 
will follow a Curie-Weiss law. This will not be the case for all elements of 
these tensors, for the symmetry of the crystal makes the contribution of wp, 
to certain elements identically zero. As an example, consider a crystal 
having the symmetry of KH,PO, (in the paraelectric phase). If the 
ferroelectric mode has #,(F)=P,(F)=0, A,(F) 49, so that q is perpen- 
dicular to the c-axis, then the macroscopic quantities which exhibit an 
anomalous temperature dependence are (a8, 08)=C yg" —C og"; Cyan = ese 
and y,,,= X3 in the usual notation. The crystal will not make a ferroelectric 
transition at T=” however, nor will the frequency w, ever become 
precisely zero. This is because the ‘internal strain’ term (af, of) 
contributes negatively to C,,” (eqn. (4.5)) and causes the latter to vanish 
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at a temperature 7'.’>7'e"._ The wavelength of the acoustic mode against 
which the crystal becomes unstable must be small compared with the 
length of the crystal in the direction of q, so that the pattern of spontaneous 
polarization and shear will be that which is found experimentally for 
KH,PO, (Kanzig 1957). If the crystal could be ‘tetragonally clamped ’ 
the transition would take place at 7'=7'.” by the crystal becoming unstable 
against the ferroelectric mode. 

These conclusions have been reached by considering the rigid ion model 
of a crystal, which is not physically realistic. When electronic polariza- 
bility is taken into account, to give eqns. (2.11), the conclusions are 
unchanged. In eqn. (4.2) the matrix M must be redefined in a way 
analogous to the definition of M’ in Part I, and Zgq replaced by a matrix 
which is not diagonal. Additional terms then appear on the right-hand 
sides of eqns. (4.5) to (4.7) (Cowley 1961b). These additional terms are 

not however related to dispersion frequencies of the lattice, although they 

can in principle be related to electronic dispersion frequencies. They would 
not therefore be expected to show more than a weak temperature 
dependence. Equation (4.9) remains unchanged in form, but now gives 
the ionic susceptibility, that is, the difference between the total and 
high frequency (or electronic) susceptibilities. 


§ 5. DiscussIon 

The first point requiring discussion is the extent to which the foregoing 
constitutes a general ‘theory’ of ferroelectricity. The concept of crystal 
stability enables the salient features of a ferroelectric (or an antiferro- 
electric) transition to be described and correlated in a very natural way. 
It leaves open the explanation of the temperature variation of the 
frequency w,, except to show that this can be brought about in a variety 
of ways such as the temperature dependence of the volume of the unit 
cell or high frequency dielectric constant, or a change in the effective values 
of the parameters specifying the short-range interaction, either because 
of the anharmonicity of the lattice vibrations or because of a change of 
bond: character with temperature, as suggested by Megaw (1952). The 
theory is based essentially on the harmonic approximation, in which for 
every mode of vibration the atoms make small oscillations about one 
equilibrium position in the paraelectric phase. As such it is probably 
adequate for ferroelectrics of the perovskite type for example. It fails 
to account for certain features of the transitions in other crystals, and the 
direction in which it must be specialized in particular cases can be 
illustrated by discussing qualitatively the special features of KH,POQ,. 
The description of the transition which the pseudo-harmonic approxi- 
mation gives for this crystal is that in the ferroelectric mode the atoms 
move as shown in fig. 4 with the centre of the hydrogen bond as the 
equilibrium position of the proton. The pattern of first-order displace- 
ments in the acoustic mode controlled by the elastic constant CO,” is the 
same (see below). Approximately this pattern of displacement is ‘frozen 
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in’ below the temperature 7'c’ at which the acoustic mode becomes unstable. 
The hydrogen atoms are not ordered from a disordered configuration. This 
simple picture conflicts with the spectroscopic evidence that the length of 
the OH bond remains nearly constant. Blinc (1960) has summarized, and 
used as a basis of a theory of ferroelectricity in certain crystals, the 
evidence that the proton in a short hydrogen bond ‘tunnels’ in a potential 
field which has a double minimum. The potential ‘hump’ however is not 


Fig. 4 


Ok 


Displacement of the atoms in the ferroelectric mode of vibration of KH,PO,, 
i.e. values of U(K, j=), deduced from the structure analysis of Bacon 
and Pease. The hydrogen atoms are shown in the centres of the hydrogen 
bonds, but it is suggested that U(K =H, j=) is nearly independent of 
the actual positions of hydrogen atoms within the bonds, as shown in 


fig. 6. 


something which is independent of the position and state of polarization 
of the other atoms. Reid (1959) has emphasized that if, as appears very 
probable, configurations with other than two hydrogens immediately 
adjacent to a phosphate group are energetically very unfavourable, then 
there is no double minimum for the motion of a single proton, but only for 
the collective movement of a group. In the polarized phase the arrange- 
ment of the hydrogen atoms is as shown in fig. 5(a) (Bacon and Pease 
1955); in the antiferroelectric structure of NH,H,PO,, according to 
Nagamiya (1952), it is as shown in fig. 5(b). These are characterized by 
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(a) Positions of the hydrogen atoms in the ferroelectric structure, showing 
spirals of hydrogen bonds all of one kind. The heights of the centres 
of the PO, tetrahedra are shown. (b) In the antiferroelectric structure 


R and L type spirals of hydrogen bonds occur alternately. Arrows 
indicate the dipole moments of H,PO, groups. 


Crystal Stability and the Theory of Ferroelectricity 417 


the arrangement of hydrogen atoms in the spirals of hydrogen bonds 
denoted by the letters R and L. Interconversion of R and L would reverse 
the polarization in both cases. We suggest that in the tetragonal phase 
the spirals R and -L, and similar arrangements producing horizontally 
directed dipoles, occur at random as illustrated in fig. 6. The dipole 
moments shown will be much smaller in this phase because of the absence 


Segment of a possible structure of the paraelectric phase. A small square in the 
centre of a PO, tetrahedron indicates a dipole moment down into the 
plane of the diagram, a small circle indicates an upward dipole moment. 
Arrows attached to hydrogen atoms indicate their direction of displace- 
ment at a particular phase of the ferroelectric mode of vibration. 


of any cooperative effect. It appears plausible that the displacements of 
the hydrogen atoms in the ferroelectric mode are small displacements 
about their various equilibrium positions, as shown in fig. 6, while the other 
atoms move as in fig. 4. The ferroelectric mode then appears, in Reid’s 
terminology, as a collective movement for which the potential hump is 
greatly reduced, so that in addition to their small movements, groups of 
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hydrogens will jump, predominantly in a direction which enhances the 
local polarization set up by the ferroelectric mode. This will occur to an 
increasing extent as the instability is approached, since this is the point 
at which the potential wells become flat-bottomed. The antiferroelectric 
transition now appears as an instability against another mode which in 
NH,H,PO, happens to have a lower frequency than the ferroelectric 
mode, although the latter is also low. Inspection shows it to be a 
transverse mode with q=(0, 0, 27/c); compare the last paragraph of § 3. 
Few details of the structures of the phases of Rochelle Salt have yet been 
published, but the indications are that the transition cannot be described 
as a simple displacive one (Okaya et al. 1960). Nevertheless there is 
experimental evidence for the occurrence of a frequency w ,/(27) ~ 101° ¢.p.s. 
in this crystal, and this frequency is lowest at the Curie points (Akao and 
Sasaki 1955). If we assume that results which have been derived using the 
harmonic approximation apply to KH,PO, at least as a first approximation, 
then the corresponding frequency in this crystal can be estimated. From 
the fact that all involve the same (temperature-dependent) elements of 
the matrix H® it can be deduced that the pattern of displacement of the 
atoms is the same in all of the following circumstances: (a) in the ferro- 
electric mode of vibration, (b) when a field is applied parallel to the 
ferroelectric axis and (c) when the crystal is sheared about this axis. It 
follows that the value of F,?(F)/c, in eqn. (4.9) can be estimated from the 
spontaneous polarization and the measured displacements of the atoms 
in the polarized phase. Neglecting the contribution of modes other than 
j=F, eqn. (4.9) leads to 
ENA sities) 
Cpuz 


Saal 
Wn = 


(5.1) 


where P is the spontaneous polarization, @ the Curie constant, p the density 
and u? the weighted mean square displacement of the atoms, that is 


wy m= = mate (K), Pe a5 as eel Oe 


where u(K) is an individual displacement. Experimental values 
are approximately @=260°, T’—T7,"=175° at room temperature, 
P=4-8yc/cm? and p=2-35g/ec. The oxygen tetrahedron is approxi- 
mately undistorted after the transition, w(pot.)—w(oxy.)= +0-04A, 
u(phos.)—w(oxy.)= —0-08A parallel to the c-axis and u(hyd.)=0-2A in 
a perpendicular direction (Bacon and Pease 1955). The estimated value 
of w,/(27) is therefore 2-5 x 101? ¢.p.s. or about 85 cm-! at room tempera- 
ture, decreasing to a minimum of about 14cm! at the transition 
temperature. These cannot be expected to be correct in more than 
order of magnitude. 

The first dispersion frequency of ys is certainly higher than 2-5 x 101° ¢.p.s. 
(work by W. A. Yager, quoted by Kinzig 1957). Very recently the infra-red 
absorption in the range 500-100 cm of solid potassium, rubidium, and 
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ammonium dyhydrogen phosphates and of potassium dihydrogen arsenate 
has been studied by Hadzi (1961) who remarks that the group of ferro- 
electrics exhibits strong absorption, extending from above 500 em-! to 
below 200 cm~. In certain materials including KH,PO, an absorption 
band was found at about 120 cm~ and did not appear in the spectrum of 
the deuterated analogues, for which 7’— 7.” is smaller. (This band has 
been attributed by Hadzi to direct transitions between the split zero 
level of the OH oscillator.) The temperature dependence and polarization 
of these bands do not appear to have been measured. 

The Curie constants of ferroelectrics cover a wide range. In a crystal 
which is diagonally cubic, @ is equal in order of magnitude to the inverse 
of the temperature coefficient 8 appearing in eqn. (3.14), irrespective of 
the number of atoms per unit cell. However this relation does not hold 
for crystals of lower symmetry. In crystals having about the same values 
of wy, eqn. (4.9) indicates that different values of @ mainly reflect different 
values of the oscillator strengths W?(F)/c,, of their ferroelectric modes. 
It is therefore possible in principle for a crystal to make a transition of the 
displacive type discussed here without exhibiting any large anomaly in 
its dielectric properties, provided that the transverse optic mode with 
g->0 against which it becomes unstable is symmetric. A comparison of the 
structures of a and f quartz strongly suggests that the transition in this 
material provides an example. 

Broadly similar conclusions on the order of magnitude and temperature 
dependence of the frequencies of vibration modes in ferroelectrics and in 
quartz have been reached by Ginzburg (1960) and by Ginzburg and 
Levanyuk (1960), whose treatment is essentially a phenomenological one. 
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‘ ERRATUM 


‘In Part I it was deduced that in the cubic perovskite structure there 
are, for values of q close to zero, three transverse optic and three longi- 
tudinal optic modes of vibration. There is in fact a fourth T.O. mode and 
a fourth L.O. mode, they have the same frequency and are symmetric. 
Having no polarization associated with them, they do not affect the 
dielectric properties of the crystal, and their omission from formulae 
given in Part Lis correct, although their existence was in fact overlooked. 
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